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PREFACE 

This  roport  was  prepared  by  Dr.  John  W.  Perry,  Archie  W,  Straitori  and  Boh  M, 
Fannin  ol  ttie  Millimeter  Wave  Gionp,  Eleclrioal  Ent;ineerinf;  Research  Laboratory, 

Tlie  University  of  Texas  at  Austin,  under  Grant  DA-ENG-27021-73-G39.  The  work  was 
perlormed  as  part  of  the  Advanced  Research  Projects  A{;ency's  Arctic  Surface  Fffect 
Veliicle  program  under  ARPA  Older  leif),  Dr.  Pieter  Hoekstra  and  Dr.  Ricliard  Munis 
monitored  the  project  for  USA  CRREL. 

Technical  review  of  the  report  was  performed  by  Dr.  Rof;er  Beif-er  of  USA  CRREL. 
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complex  index  of  refraction  at  97  ^Hz  I 

jSLr?al  !5"lo"l8orin°f"  was  carried  out  over  the  timperaturj 

C In  a specially  constructed  environmental  ^ 
chamber  of  approximately  70  cubic  foot  volume.  A Fabry-Perot 
cavity,  consisting  of  a 6-in.-diam  aluminum 
mirror  of  a 72-in.  radius  of  curvature  and  a copper 

ti  i"*  fJ’om  the  spherical  mirror,  was^used 

to  make  the  electromagnetic  measurements.  Waveguide  feed*? 
brought  through  the  back  of  the  flat  plate  co^d  Lerg?  into 
and  out  of  the  cavity.  The  complex  index  of  refractlon^of  the 
ce  fog  medium  in  the  fields  of  the  Fabry-Perot  resonant  cavltv 
were  found  from  the  change  in  the  resonant  frequency  and  ampll- 
tude  response  of  the  cavity  when  the  fog  was  present  Thr^og 
could  be  controlled  and  Introduced  to  the  cold  chamber  frL  an 

® blower  located  inside  the  chamber 
was  used  to  mix  and  dissipate  the  fog  as  desired.  A phase 

ovee^Th  ^^ed  to  sweep  the  klystron ’signal  source 

over  the  response  range  of  the  Fabry-Perot  cavity.  The  swe^^ 
frequency  was  beat  with  known  marker  frequency  oLlllaL^rto 
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13.  Abstract  (Cont'd) 

provide  reference  markers,  and  this, 
plus  the  output  of  the  loosely 
coupled  square  law  detector  connected 
to  the  Pabry-Perot  cavity,  was  put  on 
an  oscilloscope  and  photographed  for 
later  study.  The  cavity  frequency 
change  was  measured  by  positioning  the' 
marker  frequency  beats  at  the  peak  of  i 
the  cavity  response  before  and  after 
the  Ice  fog  was  admitted.  Ice  fog  ' 
particle  number-density  measurements,  i 
together  with  the  measured  extinction  ' 
of  a helium-neon  laser  (6328a),  were  ' 
used  to  help  characterize  the  fog 
properties  at  the  time  of  the  electro-j 
magnetic  measurements.  For  the  fogs 
studied  In  this  experiment,  correspond- 
ing to  1 dB/m  of  laser  extinction,  | 
measured  refractive  Index  changes  of  IJ 
to  2 N-unlts  occurred,  and  the  extlnc^ 
tlon  was  about  5 x 10“'  cm“^  at  the  ' 
lower  temperatures,  and  reached  3 or 
4 X 10  cm  at  the  -30°C  upper  tem- 
perature limit,  A literature  search 
Indicated  there  were  no  reported  com- 
plex dielectric  constant  measurements 
of  Ice  at  97  GHz.  The  complex 
dielectric  constant  of  Ice  Is  Importer 
both  for  understanding  3~mm  radio  wave 
propagation  In  an  arctic  environment 
and  for  some  of  the  theoretical  develcp 
ment  of  this  Ice  fog  research. 
Accordingly,  these  measurements  were 
made  for  laboratory  prepared  slab-lce 
samples  and  the  results  are  reported 


herein. 
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COMPLEX  REFRACTIVE  INDEX  OF  ICE  FOG  AT 
A RADIO  WAVELENGTH  OF  3 MM 

by 

John  W.  Perry,  Archie  W.  Straiten  and  Bob  M.  Fannin 


I.  INTRODUCTION 


General 

The  last  decade  has  seen  a strong  trend  toward  expanding  commercial  communications  from  the 
highly  congested  low-frequency  channels  into  the  UHF  microwave  and  millimeter  wave  spectrum. 

This  move  toward  using  higher  and  higher  frequencies  has  been  helped  along  considirably  by  the 
increasing  use  of  satellit  j-borne  millimeter  wave  transpondevs  and  the  increasing  pressure  for  higher 

data  bit  rates. 

The  millimeter  wave  region  is  a transition  region  between  the  radio  (long-wave)  domain  and  the 
infrared  region.  Many  of  the  problems  and  peculiarities  of  this  domain  are  unique.  Molecular 
absorption  becomes  important,  and  absorption  and  scattering  by  raindrops  represents  a severe  atmo- 
spheric limitation  for  communication  links.  Thermal  (blackbody)  radiation  becomes  increasingly 
more  important  as  a source  of  background  noise,  and  at  the  same  lime  millimeter  wave  hardware  such 
as  harmonic  generators,  mixers  and  sources  are  at  about  their  maximum  usable  frequencies,  particu- 
larly when  Iterating  at  90  GHz  and  above. 

Since  atmospheric  behavior  represent;  one  of  the  most  important  factors  in  millimeter  wave 
propagation,  much  effort  has  been  made  to  u.  . 'stand  and  study  its  effects.  However,  most  of  this 
effort  has  been  devoted  to  the  warm  atmosphere. 


Problem  desaiption 

Very  little  is  known  of  the  effects  of  arctic  and  subarctic  atmospheric  propagation  on  milli- 
meter waves.  It  is  important  to  understand  and  study  these  effects  because  in  recent  years  exten- 
sive efforts  have  been  underway  to  link  arctic  communities  with  the  rest  of  the  world  by  ground- 
based  and  sateUite-based  links.  Recently  a 95-GHz  pulse-type  raoar  was  planned  to  be  used  on  a 
surface  effect  vehicle  for  navigation/guidance  purposes.  This  proposed  vehicle  would  trave  a 
high  speeds  (up  to  100  knots  or  about  115  mph)  over  water,  ice  and  open  tundra  surfaces  and  smee 
it  has  limited  deceleration  characteristics  it  must  be  capable  of  detecting  ice  objects  and  other 
obstacles  far  enough  in  advance  to  be  able  to  change  course  or  slow  down.  This  3-mm  radar  looks 
horizontally  at  low  grazing  angles  and  must  be  capable  of  penetrating  fog  and  snow  (at  least  light. 

dry  snow). 

In  arctic  areas  ice  fog  is  a common  occurrence.  Ice  fog  is  a low  temperature  phenomenon,  and 
almost  never  forms  at  temperatures  above  about  -30°C.  Ice  fogs  form  where  open  sea  water  leads 
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form  in  cracks  in  sea  ice  packs:  where  warm,  moist  air  is  advected  over  cold,  dry  areas;  or  wliere 
moist  air  from  natural  liot  springs,  cooliiifr  tower  stacks,  or  industrial  or  automobile  exliausts  is  in- 
jected into  cold,  dry  air.  The  transition  from  a liquid  state  to  an  all  solid  state  takes  place  be- 
tween -30"^  iuid  about  -36‘^C,  and  it  is  in  this  repon  tliat  tlie  niaxinmm  radio  siftnal  alrsorptioii 
effects  are  expected  to  occur.  Extremely  dense  ice  fofts  liave  Ireen  olrserved  in  arctic  areas,  and 
their  extent  mid  duration  liave  proven  to  fie  very  disruptive  to  airport  and  otlier  commercial  and 
private  activities. 

Very  little  is  known  about  ttie  effects  of  ice  fojts  on  3-nini  radio  wave  propaftation,  particularly 
in  tlie  transition  temperature  interval  -30'’  to  -36°C.  It  lias  been  tlie  ftoal  of  tliis  researcli  to 
measure  the  complex  index  of  refraction,  phase  constant,  and  signal  loss  associated  with  ice  fof;s 
of  uniform  density  formed  over  the  temperature  interval  -30''  to  -48  C. 

Since  this  was  to  be  a laboratory  experiment,  it  was  necessary  to  construct  an  environmental 
cliamber  especially  suited  for  making  these  low  temperature  fog  nieasnrements.  Tlie  chamlier  so 
constructed  has  a volume  of  70  ft  * and  a temperature  control  capability  down  to  50'  C. 

Fog  patticlfc  size  iiieasuiemeiils  were  made  with  an  impact  or  to  help  characterize  the  tog  but, 
for  reasons  snown  later,  a far  more  penrate  indication  of  fog  uniformity  was  given  by  continuously 
monitoring  the  enJiiHd.luii  ^ r a ^ hsui'  bfew,.  sflM  lUhAig!!  a preBcribfctl  path  in  the 

ice  fog.  In  this  manner,  a “reference  fog  level’’  can  be  readily  attained  for  each  operating  tempera- 
ture. 

A high-Q  Fabry-Perot  cavity  was  designed  mid  constructed  for  making  the  electromagnetic 
measurements  to  determine  the  complex  index  of  refraction  of  the  fog.  The  resonant  frequency 
change  and  signal  amplitude  chmige  of  the  cavity  due  to  the  introduction  of  the  fog  were  measured 
by  a nietliod  devised  for  pluise  locking  and  sweeping  the  klystron  signal  source  over  the  cavity 
response  frequency  range.  A high  degree  of  frequency  stability  and  precision  frequency  measure- 
ment was  attained  in  tliis  manner. 

As  a suppleiiieiit  to  these  electi oiiiagiiet ic  ice  log  measurements,  the  complex  dielectric 
constant  of  slab  ice  was  measured.  The  authors  could  not  find  any  repoited  dielectric  constant 
measurements  of  ice  at  Jl'  mm  Tbnia  ire  »i  ulwtuus  assol  for  tin 

lation  of  ice  reflections,  and  tlie  backscattering  and  extinction  cross  sections  of  hail  at  3 mm. 

Tliese  ice  measurements  were  made  separately  using  a free-space  quasi-optical  technique  commonly 
described  in  the  literature. 

Method  of  presentation 

As  far  as  the  presentation  of  the  material  is  concerned,  it  was  deemed  appropriate  to  present 
first  the  results  for  the  slab  ic3  complex  dielectric  constant  measurements.  The  literature  review 
of  solid  ice  permittivity  measurements,  instrumentation,  and  calibration  techniques  was  incorporated 
in  the  same  chapter  (Chapter  II). 

Chapter  III  is  concerned  witli  tlie  meteorology  and  the  appropriate  literature  research  of  arctic 
ice  fog  formation. 

Chapter  IV  is  devoted  to  a theoretical  development  of  the  dielectric  behavior  of  a gas  ’vith 
homogeneously  distributed  small  dielectric  spheres. 

Chapl«fs  V and  VI  ate  devoted  to  describing  thbotelically  the  Fabry-Petol  cavity  useJ  di  these 
measurements  and  the  means  of  electronically  generating  and  sweeping  the  signal  over  the  response 
frequency  of  the  cavity.  It  is  in  this  section  m iiutmui'LiiitJitl.ii  Utat  Ihu  Mftifh ftu iiX'  rlmiiibur  and 
laser  extinction  measurement  methods  are  described.  Literature  surveys  are  included  where 
appropriate. 

ir 
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The  ineasuienients  and  results  ate  treated  fully  in  Chapter  Vll.  The  early  sections  in  this 
chapter  are  devoted  to  aiialyziiit^  empirically  the  effects  of  ice  accumulation  on  the  mirrors.  Data 
analysis  and  application  to  one-way  sif^nal  propagation  are  treated  here,  and  the  results  appropriate 
to  propagation  through  fogs  of  a given  reference  laser  extinction  are  plotted  and  tabulated. 


II.  COMPLEX  PERMITTIVITY  MEASUREMENTS  OF  SOLID  ICE 


Introduction 

In  microwave  remote  sensing  applications,  the  dielectric  constant  of  the  material  being  sensed 
is  a very  impoiiaiit  parameter.  The  backscattenng  cross  section  of  an  ice  spbere  in  a thunderstorm, 
when  struck  by  a radar  wave,  is  a function  of  the  sphere  circumference  in  wavelengths  and  the 
dielectric  constant  of  the  material.  Mie's  (1908)  treatment  of  tlie  scattering  and  extinction  cross 
sections  of  dielectric  spheres  was  restated  by  Stratton  (1941),  and  numerical  methods  of  machine 
computing  the  n^''  electric  and  magnetic  modes  involved  in  these  calculations  have  been  devised 
(Aden  1950). 

In  order  to  model  the  complex  permittivity  of  a low  temperature  ice  fog,  it  is  necessary  to  have 
a knowledge  of  the  complex  permittivity  of  ice.  No  published  data  appropriate  to  the  operating 
frequency  of  96  GHz  exist.  Inference  from  published  data  at  other  frequencies  leads  to  inconsistent 
conclusions.  It  is  therefore  appropriate  to  make  dielectric  constant  measurements  of  ice  at  96  GHz. 
This  chapter  contains  a literature  survey  of  the  available  publislied  data  for  ice,  and  details  the 
theory  of  dielectric  measurements  and  the  instrumentation,  calibrations  and  results  for  ice  at  this 
frequency. 

Review  of  the  literature  on  dielectric  behavior  of  ice 

Ice  can  occur  in  six  allotropic  forms  for  pressures  in  excess  of  2000  atmosplieres  (Bertie  and 
Wlialley  19b4),  and  for  pressures  ol  luOO  atmospheres  or  less,  ice  can  occur  in  the  noniial  hexagonal 
(1),),  cubic  (/p),  or  amorphous  state,  depending  upon  temperature  and  the  method  of  preparation 
(Bertie  and  Whalley  1964,  Lonsdale  19.58).  Between  about  0°  and  -80'’C  the  equilibrium  form  of  ic:e 
is  the  normal  hexagonal  form  /,j.  For  temperatures  below  -80”C  and  aliove  -130°C  ice  may  exist  in 
the  cubic  form  1^,  and  below  about  -130'^C  ice  may  be  in  the  amorphous  form.  These  temperatures 
are  opprjxtinalt  thsagre  .-metit  iX-Hif*  h.  tlri  Ihitalwc  aa  fc  the  exitet  butH.dtwy  vaiueh.  Huw 
evei,  as  it  occurs  naturally  in  the  terrestrial  environment,  the  normal  form  of  crystalline  ice  is  1^^. 
While  it  is  important  to  recognize  what  possible  crystalline  forms  ice  can  take,  it  is  not  important 
to  the  purpose  of  this  research  to  further  discuss  these  forms  since  for  all  practical  temperatures 
and  forms  (slab  ice,  snow,  hail  and  ice  fog)  only  /,|  exists. 

Irvine  and  Pollack  (1968)  reviewed  the  literature  on  the  absorption  coefficient  and  reflectivity 
of  both  water  and  ice,  and  for  each  phase  they  tabulated  both  the  real  and  iinaglnafy  parts  of  tlie 
index  of  refraction  (ra  = n'  - /n")  for  wavelengths  between  0.95  and  152  microns  Their  values  of 
the  complex  index  of  refraction  are  well  accepted  in  this  wavelength  interval. 

The  microwave  values  for  the  complex  refractive  index  of  ice  given  by  Gunn  and  East  (1954) 
are  for  n'  = 1.78  at  all  temperatures  and  n"  -24  x lO”**  at  0°C,  8 x 10“'*  at  -10°C,  and  6 x 10““*  at 
-20'^ C.  These  values  ate  based  upon  the  work  of  Dunsmuir  and  Lamb  (1945),  Lamb  (1946),  Lamb 
and  Turney  (1949),  and  Gumming  (1952)  over  ine  centimeter  band.  Lamb  (1946)  used  a cavity 
operating  at  3 cm  (10  GHz)  and  found  a relative  permittivity  of  3.05.  However,  he  later  found  a 4% 
cTFuT  h.  his  njeMHtdfiietrtH  tiue  lo  from  cup  in  c- wity.  wal  Lw.ib  md  Tiwuey 

(1949)  repeated  these  cavity  dielectric  measurements  at  3-cm  and  1.25-cm  wavelengths  (10  and  24 
GHz).  They  reported  a relative  permittivity  near  3.17  at  both  frequencies  near  0^'G,  and  this  value 
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changed  only  slightly  for  temperatures  down  to  -195'^C,  Dunsniuir  and  Lamb  (1945)  and  Lamb 
(1946)  reported  that  the  loss  tangent  of  their  ice  samples  varied  from  12  x 10"'*  at  0°C  to  1 9 x 
10  ^ at  -20°C,  and  to  nearly  1.3  x 10“'*  at  -40°C. 

Gumming  (1952)  later  used  a measuring  technique  employing  a section  of  waveguide  filled  with 
ice  to  make  permittivity  and  loss  measurements  at  3.2-cm  wavelength.  Gumming  found  a relative 
permittivity  of  3.15  independent  of  temperature  (over  a temperature  range  of  0°  to  -20°G)  at  tliis 
wavelength.  In  terms  of  tlie  imaginary  part  of  the  complex  index  of  refraction,  Gumming  found  loss 
tangent  values  about  twice  those  of  Lamb,  with  measurements  of  23  x 10"'*  at  0"G,  7,9  x 10"“*  at 
-10  C,  and  5.5  x 10  "*  at  -20"G.  His  measurements  have  been  considered  to  be  the  liest.  even 
though  made  only  over  a 0°  to  -20°G  temperature  range  and  at  3.2-cm  wavelength. 

Perry  and  Sti^iton  (1972)  reported  making  complex  permittivity  measurements  on  slab  ice  at 
35  and  94  GHz.  These  measurements  were  made  using  a qnasi-optical  free  space  technique.  A 
iscussion  of  tliis  technique  is  deferred  to  a later  section  dealing  witli  tlie  theory  of  dielectric 
measurements  as  made  on  slab  ice  at  96  GHz.  The  results  presented  by  Perry  and  Straiten  for  the 
relative  dielectric  constant.  1.91.  are  better  than  30%  below  the  low  frequency  value  given  by  Lamb 
etc.  Since  the  temperature  control  exerted  over  these  ice  samples  was  poor,  and  since  an  error  in 
the  phase  shift  interpretation  was  found,  these  measurements  were  repeatr'ci  in  an  environmental 
chamber.  The  results  of  these  recent  measurements  are  deferred  to  a later  section. 

Von  Hippel  (1954)  reports  measurements  made  by  Westphal  on  ice  at  -12°G  over  the  frequency 
range  0.1  MHz  to  10  GHz.  His  results.  3.20  at  3 GHz  and  3.17  at  10  GHz,  are  in  accord  with  the 
above  results  of  Gumming.  Between  1 MHz  and  1 GHz  his  measured  results  were  surprisingly  high, 
being  4.8  at  0.1  MHz.  4.15  at  1 MHz,  and  3.7  at  10  MHz.  At  these  frequencies  tan  B was  found  to 
be  9 / 10  . 7 X 10  '*,  8000  x lO""*.  1200  x 10"'*  and  180  x 10  "*.  respectively. 

Yoshino  (1961)  made  dielectric  constant  measurements  between  1 MHz  and  3 GHz  on  Antarctic 
ice  core  samples  over  the  temperature  range  -18°  to  -36°C.  His  results  are  in  accord  with 
Von  Hippel  s measurements  over  this  frequency  range,  being  nearl,’  4.0  at  5 MHz,  3.8  at  10  MHz 
3.45  at  100  MHz  and  about  3.3  at  3 GHz. 

These  results  in  the  intermediate  (HF-UHF)  frequency  range  of  both  Von  Hippel  (1954)  and 
Yoshino  (1961)  are  surprisingly  high  and  in  good  agreement  with  each  other  in  this  interval. 

More  recently.  Ragle.  Blair  and  Persson  (1964)  reported  on  Westphal’s  dielectric  measurements 
of  glacier  ice  made  from  150  MHz  to  iGHz  and  from  1°  to  -60°G  (Table  1).  These  measurements  are 
significantly  below  what  Von  Hippel  (1954)  and  Yoshino  (1961)  had  reported  over  this  frequency 
range,  and  aside  from  the  decrease  in  the  second  significant  figure  of  the  dielectric  constant  with 
increasing  frequency,  the  dielectric  constant  is  approximately  constant  at  2.9.  both  with  temperature 
and  frequency.  The  loss  tangent,  however,  decreases  both  with  decreasing  temperature  and  increas- 
ing frequency,  except  at  1 GHz,  where  it  is  reported  to  be  constant  with  decreasing  temperature. 

Ray  (1971)  made  an  empirical  model  of  the  complex  refractive  indices  for  ice  and  liquid  water 
which  is  applicable  from  -20°  to  0°G  for  ice  and  from  -20°  to  50°G  for  water.  His  model  extends 
from  2 /i  to  ! everal  thousand  kilometers  in  wavelength  for  ice.  His  expressions  are  of  value  in  that 
they  fit  quite  well  the  majority  of  the  absorption  bands  in  the  infrared  and  optical  parts  of  the 
spectrum.  In  the  radio  frequency  part  of  the  spectrum  he  adjusted  his  expressions  to  best  fit  tlie 
data  of  Gumming  (1952),  Lamli  (1946)  and  Lamb  and  Turney  (1949),  and  for  a single  relaxation  time 
effect  and  no  losses  {o  - 0)  his  model  reduces  to  the  classical  Debye  expression. 

Worz  and  Gole  (1969)  made  static  dielectric  constant  measurements  of  redistilled  water  samples 
further^ purified  by  zone  refining  techniques  to  yield  dc  conductivities  as  low  as  4 x 10"**  mho/cm 
at  -10  G.  Their  measurements  were  made  over  the  temperature  interval  0°  to  -80°G.  From  their 
measured  data,  the  equation  best  describing  the  static  dielectric  constant  is  of  the  form 
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Table  I.  Complex  permittivity  of  glacier  ice  core  samples 
from  150  MHz  to  1 GHz. 

(After  nieasureitienis  made  by  ‘(Vesiphal  and  reiioned  by  Ragle, 
Blair  and  Persson  1964.) 

Temp 


'O 


rc) 

150  MHz 

300  MHz 

500  MHz 

1000  MHz 

Dielectric  constant 

1 

2.95 

2.94 

2.93 

2.93 

10 

2.92 

2.92 

2.915 

2.915 

-20 

2.91 

2.91 

2.91 

2.91 

-00 

2.90 

2.90 

2.90 

2.90 

Loss  tangent 

- 1 

0.00200 

0.00108 

0.000.52 

0.0004 

-10 

0.00110 

0.00056 

0.00028 

0.0004 

-20 

0.00068 

0.000,34 

0.00019 

0.0004 

-60 

0,00005 

0.00003 

0.00003 

0.0004 

20,715 

7’  -38 


and  is  taken  to  be  3.1.  The  relaxation  time  of  the  ice  thus  prepared  was  foimd  to  obey  an 
Arrhenius  rate  eqitation  of  the  form 

where  A ,'3.3  x 10' sec 

C 13.2  kcal^mole  r,.r,2657(;  ' lO"*  J/niole 
R 8.3143  J/K  - mole 

and  is  in  good  agreement  with  the  relaxation  limes  measured  by  Auty  and  Cole  (1952)  in  the  tem- 
perature range  above  -50°C.  Worz  and  Cole  (1969)  give  slightly  higher  values  for  the  static 
dielectric  constant  of  ice  at  very  low  temperatures,  and  their  results  are  probably  better  than  Auty 

and  Cole  s (1952)  findings  (0°  to  -66°C)  due  to  the  care  taken  to  reduce  the  dc  conductivity  of 
their  water  samples. 

n 'understand  what  is  meant  by  the  relaxation  time  is  to  visualize  the 

molecular  dipoles  as  tiny  spheres,  as  Debye  (1929)  did  in  his  original  theory.  The  rotation  of  these 
spheres  in  a viscous  medium  is  opposed  by  forces  related  by  Stoke’s  law  to  the  macroscopic 
viscosity  Tj  of  the  surrounding  medium,  and  Debye  (1929)  deduced  that 

4nr]ii^ 

kT 


where  a molecular  radius. 
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It  is  interesting  to  put  some  niimi)ers.  appropriate  for  water,  "ito  this  simple  expression  for 
relaxation  time.  Water,  at  room  temperature  (about  ^20''C),  has  a viscosity  of  n 0.01,  and  the 
’■radius”  of  the  water  molecule  is  about  2 x 10“^*^  ni,  giving  a r of  0.25  x 10~^“  sec.  The  relaxa- 
tion time  of  water  is  in  fact  associated  w th  an  absorption  peak  in  tlie  centimeter  baud.  This 
"spliere”  model  of  Debye’s  is  really  only  important  in  that  it  postulates  that  tlie  orientation  of 
polar  molecules  in  liquids  and  solids  leads  spectroscopically  to  a simple  relaxation  spectrum 
whicli  is  indeed  wliat  is  observed. 

More  recent  approaclies  to  tlie  problem  of  dipole  relaxation  leave  involved  the  idea  that  the 
dipole  activity  rotates  between  two  positions  of  equilibrium  sepaiated  by  a potential  energy  barrier. 
If  the  height  of  this  energy  barrier  is  written  as  D,  then  the  relaxation  time  - which  is  a measure 
of  the  transition  probability  - is  given  by  a relation  of  the  form 

pD/kT  _ E/RT 
f - e Tq  e 

Ui 

where  R - kyM,C\^  only  slightly  t»u*.i..'jature  dei;r.!,denl,  and  rr/w  is  th/e  average  time  requited  by 
an  excited  molecule  to  turn  from  one  equilibrium  direction  to  the  other  (Saxton  1952). 

Ihe  relaxation  time  of  ice  increases  with  decreasing  temperature.  Auty  and  Cole  (1952)  and 
'iote  and  Cult  ( iyt! it)  tevjtl  t to  ratige  Irom  abetrt  2 x w'*  at  u C lo  2.5  x lu‘  sec  at  -Hi>  C. 

It  is  interesting  to  compare  these  relaxation  times  of  ice  with  the  relaxation  times  of  water.  Saxton 
(1952)  has  measured  r of  water  to  lie  1.87  x 10'^’  at  0"C  to  about  6 x 10  at  40^C.  Thus  water 

has  a PblaxatitMi  Hint,  abuut  lu  that  of  ice. 

At  least  in  a first  approximation,  both  water  and  ice  traverse  a simple  Debye  relaxation 
spectrum  since  only  one  relaxation  time  has  been  found  for  eatdi  phase  over  the  known  radio  wave 
( tru.fl  iitid  fC..,uEri,attii  191*9).  Thus  the  rtlutlvt  Jiclectii^  eonsl.din  i /iq  and  dielectric 

loss  ("/fp  are  given  by  tlie  expressions 


'o  'o  ,y(l  , oi^r^) 
'O  (q(1  t w^r^) 


where  is  the  limiting  value  of  the  high  frequency  dielectric  constant  and  f'  /fn  is  the  static 

(de)  dielectric  constant. 

The  principal  dielectric  loss  of  water  at  low  frequencies  (frequencies  below  100  MKz)  is  due 
primarily  to  free  cliarge  carriers,  such  as  ions  and  other  free  radicals.  This  dielectric  loss  mech- 
anism, caused  by  tlie  dc  conductivity  of  the  water,  decreases  with  increasing  frequency  according 
to  the  relationship 


f"(oj)  /; 


'0  'o'' 


where  o is  measured  in  mhos/nieter. 

At  high  frequencies,  the  dielectric  loss  of  water  is  due  primarily  to  viscous  friction  losses 
associated  with  the  rotation  of  the  polarized  water  molecules  as  they  follow  the  incident  electro- 
magnetic radiation.  Ions  and  other  radicals  cannot  contribute  significantly  to  the  dielectric  loss 
since  they  aie  just  too  massive  to  frV.ow  the  incident  electric  field  variations. 
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Figure  1 shows  the  dielectric  spectrum  of  ice  and  water  at  various  temperatures  us  calculated 
from  the  Debye  equations  with  the  appropriate  relaxation  times  and  static  dielectric  constants.  For 
ice,  the  Worz-Cole  expression  for  the  static  dielectric  constant  was  used,  aud  tlie  relaxation  times 
for  each  temperature  were  taken  directly  from  Auty  and  Cole’s  paper.  In  the  case  of  ice,  was 
taken  to  he  .3.1.  For  water,  tlie  relaxation  times  and  static  dielectric  constants  appropriate  for 
eacli  temperature  were  obtained  from  a tabulation  in  Kerr’s  hook  (1951).  The  value  for  for  water 
appropriate  for  all  temperatures  was  taken  to  lie  5.5.  In  plotting  these  curves,  the  effect  of  the  dc 
conductivity  on  dielectric  loss  was  assumed  to  be  negligible.  The  dielcciiic  spectrum  of  water  at 
0°.  10'’  and  20°C  and  ice  at  -0.1°,  -20.9°  and  32'’C  was  plotted. 

Table  11  suiumarizes  tiie  reported  values  of  the  complex  permittivity  for  ice  over  the  electro- 
magnetic spectrum. 

Theory  of  dielectric  measurements 

For  general  nonmagnetic  materials  where  /<g.  two  independent  experimental  quantities 
must  be  found  to  allow  the  analytic  calculation  for  the  simultaneous  solution  of  r'  and  f". 

The  propagation  constant,  detined  (or  the  field  quantities  E or  H,  can  he  expressed  as 
1 ’ = a + 

where  a is  tlie  attenuation  in  nepers/m  and  fi  is  the  phase  constant  in  radians/m.  The  propagation 
constant  of  a material  is  also  written 


V = = j(o 

where  is  the  complex  permittivity  and  is  .vritten 
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Table  II.  Reported  complex  permittivity  of  ice. 


k' 

k" 

t(or  X.) 

Remarks 

1.09 

0.2  X 10*^’ 

1 micron  A 

Irvine  and  Pollack  (1968) 

1.3^7 

9.2  X 10 

10  micron  A 

Irvine  and  Pollack  (1968) 

2.12 

1 

o 

> 

71 

100  micron  A 

Irvine  and  Pollack  (1968) 

3.00 

' 3.00  ■ 10  ^ 

3 Gllz 

Dunsinnir  aiid  Lamb  (1945)  iind  Lamb  (1940^);  mea- 
sured lull  S over  0°C  to  40^C  range 

3.17 

— 

10.  24  Gllz 

Lamb  and  Turney  (1949) 

3.10 

'7.0  < lO”-’ 

9.375  Gllz 

Cumming  (1952) 

1.91 

■ 0.0  X 10  ^ 

30.  94  Gllz 

Perry  iind  Straiton  (1972  ’ >)or  control  over  temper- 
ature and  an  error  in  phase  sliilt  calculation 

3.20 

2.88  X 10”^ 

3 Gllz 

Von  llippel  (1954)  reporting  .some  of  WestpliaPf; 
early  measurements  of  ice  at  -12°C 

3.17 

2.24  -■  10”^ 

10  Gllz 

Von  llippel  (1954)  reporting  sonir  of  Westplial's 
early  measurements  of  ice  at  -l2'’C 

4.80 

3.84 

0.1  Mllz 

Von  llippel  (1954)  reporting  some  of  Westplial's 
early  measurements  of  ice  at  -12  C 

4.10 

0.498 

1.0  Mllz 

Von  llippel  (1954)  reiiorting  some  of  Westphal’s 
early  measurements  of  ice  at  -12'^C 

3.70 

0.0007 

10  Mil/ 

Von  llippel  (1954)  reporting  some  of  WestphaPs 
early  measurements  of  ice  at  -12°C 

M.O 

5 Mllz 

Yoshino  (1961);  Antarctic  ice  core  samples 

3.3 

- 

3 Gllz 

Yoshino  (1961);  Antarctic  ice  core  samples 

2.9  to 

2 X 10  X lO"'"’ 

0. 10  10 

Ragle,  Blair,  and  Persson  (1964)  reporUng  measure. 

2.95 

to  4 > 10^ 

1.0  Gllz 

meats  of  West))lial  on  ice  core  samples  over  tlie 
temperature  range  -I'^C  to  -60°C 

91.0 

- 

dc 

Worz  and  Cole  (1969);  Anty  and  Cole  (19.52),  T - 

0.1  c 

< 


c 


I I << 


and  thus 


‘o  ^0  ^0 


k'  - jk" 


and  ..  is  the  permittivity  of  free  space,  8.854  x 10“^*^  farads/m.  The  quantity  k'  is  referred  to  as 
the  relative  permittivity  or  the  dielectric  constant  of  the  material,  and  k"  is  associated  with  the 
dielectric  loss  of  the  material  and  is  known  as  the  relative  dielectric  loss  factor. 

By  equating  these  expressions  for  the  propagation  constant  I , 

a ^ ifi  j(o  - jk") 


then 


+ iZdfi  -o/pfo(k'  - ;k") 
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an  1 equating  tho  rsal  and  imaginary  parts  ^ives 


k" 


2afi 


Tims  tlie  complex  dielectric  constant  may  be  computed  from  a knowledge  of  tlic  measured  attenuation 
and  pliase  sliift  of  the  material. 

A wide  variety  of  metliods  are  used  for  obtaining  an  independent  pair  of  tliesc  quantities.  Fol- 
low frequencies  one  teclinique  places  the  dielectric  material  between  the  parallel  plates  of  a 
capacitor  in  a resonant  circuit.  From  a measure  of  the  Q of  tlie  resonant  circuit  with  and  witlioiit 
the  d.elecfic  material  between  the  capacitor  plates,  and  from  a measure  of  the  frequency  shift  of 
the  circuit,  one  can  calculate  k'  and  fc".  As  the  frequency  of  measurement  increases  tlie  cavity 
replaces  the  resonant  lumped-parameter  circuit.  Other  techniques,  such  as  standi.ig  wave  measure- 
ments made  in  conjunction  with  the  null  shift,  have  lieen  successfully  employed  (see,  for  example. 
Von  Hippel  1954). 

One  approach  successfully  used  at  X-band  with  long  (in  terms  of  wavelength)  samples  is  to  fill 
a section  of  waveguide  with  the  dielectric  material  and  to  compaiu  the  phnse  shift  and  attenuation 
through  the  sample  with  that  through  an  empty  section  of  guide  (Straiten  and  Tolbert  1948.  Wiebe 
1971).  For  yet  higher  frequencies,  where  the  dimensions  of  the  guide  prove  too  small  for  practical 
applications  of  this  technique,  tlie  wave  is  propagated  throiigli  a free-space  path  and  througli  a 
section  of  waveguide  to  a detector  where  the  signals  from  each  path  are  combined  (Hertel,  Straiten 
and  Tolbert  1953.  Montgomery  1947,  Saxton  et  al.  1962).  Tlie  amplitude  change  of  the  signal 
propagated  through  the  free-space  and  material  sample  gives  information  about  the  attenuation 
constant,  and  the  null  shift  at  the  detector  gives  a measure  of  the  phase  shift  due  to  the  sample. 

This  particular  technique  was  used  to  make  the  dielectric  measurements  of  ice  at  96.5  GHz,  and  a 
more  thorough  discussion  of  this  appro.'  'h  is  deferred  to  a following  section. 

In  practice,  a number  of  material  samp...  thicknesses  are  used  since  in  general  a plot  of  attenua- 
tion versus  sample  length  will  not  be  a linear  function  due  to  variations  in  the  wave  impedance  as 
the  electrical  length  of  the  sample  is  changed.  The  attenuation  and  phase  shift  factor  are  obtained 
by  plotting  a large  number  of  points  and  taking  the  slope  of  a straight  line  drawn  as  the  average  of 
file  poinlH  Tluit  this  Ls  irun  ran  bi  .seen  from  lhn  following  argument  In  Appendix  A it  is  shown 
that  for  a plane  electromagnetic  wave  normally  incident  upon  a three-layer  dielectric,  with  tlie 
bounding  dielectric  materials  having  a complex  propagation  constant  1'^,  and  with  the  inner,  sand- 
WK!l«eJ  layer  having  a oomplix  prupagatkxi  ttetmtatit  I'g.  Hhn*  rtutio  ^ llw  wnplUtnirt  of 
advancing  components  of  tlie  electromagnetic  field  incident  upon  and  leaving  the  layered  dielectric 
are  given  by 
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Figure  2.  Geometry  of  dielectric  muterial  holder. 


For  ice  no  "holdinn”  dielectric  slabs  are 
neeesstiry.  Therefore,  r,  r/^  at.d  rf  • 0 and 
the  expression  collapses  to 


- (i/o-V„)”e"'  "^1 


where  1’^  piopagali on  constant  of  the 
dielectric 

r/,j  --  intrinsic  wave  impedance  of  the 
dielectric. 


Since  for  a low  loss  dielectric  1'^  then 


cos  (/^/)  , /2(,/‘5  r ril;  sin  . 
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Thus 


tan  0 = 


tan  (/a^  L) 


and  so 


^*3  "'-'/o'/d 

and  the  advancing  component  in  medium  3 ca 


can  be  written  as 

E = /I,  q-!P{L  + z) 

li 

a 

itc!:  fp's  r:  .tirr,: . I'j  ^ 'n 

plot  of  a and  </>,  m nepers/m  and  radians/m  respective^lv”'wili°'^  f k ^ measur 

a slowly  varying  sinusoidal  function  about  a straight  r ^ function,  but  rath( 
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figure  4.  Qunsi-opticnl  ftee-space  method  of  measuring  complex  permittivity,  p and  r . 

wavelength.  This  ,s  in  agreement  with  the  measured  data  for  ice,  whieli  vary  from  0.02  to  about 
0.06  wavelength  above  and  Ixilow  the  straight  line  fit. 

Instrumentation  for  ice  dielectric  measurements 

at  fo''  >»^king  complex  dielectric  constant  measurements 

at  96.5  GHz  involves  measuring  the  amplitude  and  phase  change  of  a signal  propagated  throuch  a 
ree  .space  path  containing  varying  thicknesses  of  the  dielectric  material.  With  this  technique  a 
.meter  wave  signal  from  a klystron  source  is  split  by  a -3  dB  coupler  and  travSs  ove^?::  ...e- 

fain  n I”"  attenuator,  but  only  the  attenuator  in  the  path  con- 

tain.iig  the  dielectric  sample  between  the  two  horns  must  be  a precision  attenuator.  After  pa  s^ 

i:  iTed  line"  ^ ^ section  of  p"i^ 

n . ■ indicates  this  experiment;  1 measurement  technique.  The  isolator  is  to 

c rnres''?n  0^1^^  '''"'"'I"" 

renecL  effect^  ' ^ 

The  horns  are  made  of  electroformed  copper,  with  an  opening  of  about  V.  in  The  minimum  horn 
separation  distance  is  set  by  the  Fraunhofer  far  field  criterion,  d 2AD2  ^ 

penin^  in  wavcleiij^tlis  and  A is^Uie  operating  wavelength,  ^ 

for  balanced  by  removing  the  slotted  line  section  and  adjusting  the  attenua- 

o s for  equal  signal  levels  over  each  path.  For  phase  measurements,  the  attenuators  me  not 
touched  again,  since  it  was  discovered  that  changing  the  attenuator  setting  by  several  dB  would 

h » ^„r " •>“-11 

11.  the  otanding  wave  pattern.  However,  a high  degree  of  balance  is  not  really  required,  since  the 
detection  system  was  sensitive  enough  to  follow  a null  shift  even  when  more  that.  6 dB  of  imbalance 

thickest  ice  samples  used,  was  under  2 dB.  The  measurements  with  water  presented  the  greatest 
problem.  The  details  of  how  this  measurement  was  accomplished  are  delayed  to  a later  section,  but 
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for  water  samples  less  than  about  0.010  in,  thick  satisfactory 
attenuation  and  phase  measurements  were  made. 

F’or  attc  luation  measurements,  propapation  er.iy  over  the  free 
space  patli  is  important.  A power  detector  wa  j mounted  on  the 
receivinp  horn,  and  the  attenuation  measurements  were  accom- 
plished in  either  one  of  two  ways.  For  moderate  and  high  loss 
dielectrics,  tlie  dB  change  in  the  precision  attenuator  setting 
necessary  to  restore  the  same  power  level  at  tlie  det'jclor  with  the 
dielectric  in  the  path  as  was  present  without  the  dielectric  is  a 
measure  of  the  dielectric  loss.  For  low  loss  dielectrics,  where  a 
cliange  in  the  precision  attenuator  setting  is  iurd  to  read,  the 
dielectric  loss  was  calculated  from  the  change  in  the  output  meter 
reading  on  the  high-gain  tuned  amplifier.  This  amplifier  was 
checked  against  various  dB  settings  of  the  precision  attenuator,  and  its  output  meter  reading  was 
very  linear  in  power,  thus  indicating  that  the  crystal  detector  was  operating  in  its  square-law  region. 
With  several  hours  warm-up,  both  the  klystron  and  this  amplifier  proved  to  be  quite  stable,  and  over 
tlifc  conrse  of  the  n’leaKurenieRtB  the  output  luvul  wturnod  to  its  initial  value  on  the  output  na  tei  woh 
the  initial  (reference)  precision  attenuator  setting  and  no  dielectric  in  the  path. 

To  make  phase  shift  measurements,  each  propagation  path  is  connected  up  to  one  side  of  the 
slotted  line  section  after  having  been  balanced  as  well  as  possible,  as  described  earlier.  The  probe 
in  the  slotted  line  section  is  adjusted  for  minimum  penetration  depth  consistent  with  a workal)le 
signal  level.  This  assures  a minimum  perturbation  of  the  fields  within  the  waveguide.  With  no 
dielectric  material  in  the  path,  the  probe  position  is  adjusted  for  a null  in  the  standing  wave  pattern. 

A dielectric  sheet  placed  in  this  free  space  path  causes  this  null  to  shift  in  position,  and  the 
amount  of  shift,  in  hundredths  of  a millimeter,  can  be  read  by  the  gauge  attached  to  the  carriage  of 
the  slotted  line.  This  null  shift,  in  millimeters  or  fractions  thereof,  can  then  be  put  in  terms  of 
fractional  parts  of  a wavelength  shift  and  plotted  as  a function  of  dielectric  thickness,  as  is  usually 
done. 

With  a slotted  line  technique  of  making  phase  measurements,  the  measured  phase  shift  is  actually 
half  the  true  phase  shift.  Let  E^  and  E^  be  the  reference  and  the  transmitted  E-field  vectors  in  the 
slotted  line  at  the  probe  for  no  dielectric  in  the  free-space  path  and  for  the  probe  positioned  at  a 
null  in  the  field.  With  the  introduction  of  the  dielectric  slab  in  the  free-space  path,  the  phase  of  the 
transmitted  wave  shifts  E^  to  E|  (see  Fig.  5)  by  an  amount  0.  To  reestablish  a null,  the  probe  is  re- 
adjusted to  shift  the  phase  by  6/2.  since  in  moving  E|  to  E^'  through  an  angle  0/2,  E^  is  moved  6/2 
in  the  opposite  direction  to  E^,  thus  making  E’  and  E|^  again  180  degrees  out  of  phase  with  each 
other. 

For  all  dielectric  measurements  made,  the  klystrom  repeller  was  modulated  at  2 kHz,  and  the 
high-gain  amplifier  used  after  the  detector  was  sharply  tuned  to  2 kHz.  This  procedure  is  usually 
followed  to  avoid  the  problems  of  flicker  noise  (If  noise)  and  drift  associated  with  a high  gain  dc 
amplifier.  There  is  no  special  reason  to  choose  a 2-kHz  modulation  frequency,  but  in  general  the 
modulation  frequency  should  be  chosen  to  avoid  the  60-Hz  line  frequency  and  its  harmonics. 

Calibrations  with  polystyrene  and  water  samples 

In  order  to  establish  the  accuracy  and  to  becon»  familiar  with  this  dielectric  measurement 
technique,  dielectric  measureiiieiits  were  uiade  Ofi  polystitrciic  tOid  wattf  prk*  to  the  iec  tiieasiiri!  ■ 
ments.  The  measurements  were  all  made  at  room  temperature. 
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Figure  ,5.  Phase  relationship 
of  transmitted  and  reference 
fields  in  the  slotted  line  sec- 
tion for  the  null  condition. 
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Figure  6.  Polystyrene  phase  shift  curve. 

Polystyrene.  Slabs  of  polystyrene  were  made  by  niaebining  disks  from  a 13-iii.-diam  polystyrene 
rod.  l-ileven  sucb  disks  were  macliined  in  thicknesses  I’roin  0.025  in.  to  over  0.5  in.  Figure  6 is  a 
plot  of  the  phase  shift  curve  for  the  polystyrene  samples.  The  error  bars  were  set  by  the  uncertainty 
in  measuritin  the  null  shift  with  the  slotted  line  and  in  the  variations  of  the  observed  data  points  from 
repeated  measurumeiits.  Finure  7 is  the  ailenualioii  curve  measured  for  this  polystyrene  sample  set. 
The  error  in  these  data  is  about  * 0.1  dB,  This  inaccuracy  arises  principally  due  to  cliaiit’es  in  the 


tilt  angle  of  the  dielectric  slabs.  Since  it  was  impossible  to  orient  the  samples  exactly  the  same 


each  time,  many  repeated  measurements  were  made  to  average  out  this  effect. 

Polystyrene  was  chosen  for  calibration  purposes  since  its  complex  dielectric  constiuit  is  well 
known  and.  except  for  its  dielectric  loss,  essentially  constant  (ITT  Reference  Data  for  Radio 
Engineers  Handbook,  1970).  The  measured  and  published  complex  dielectric  constants  for  this 
material  are  given  in  Table  III. 

Water.  Water  was  chosen  for  measurement  purposes  because  its  behavior  is  well  described  by 
the  Deliye  expressions.  In  addition,  measurements  liave  been  made  of  the  complex  dielectric 
constant  at  millimeter  wave  frequencies  :uid  reported  in  the  literature  (Zolotarev  et  al.  1969, 
Cliamlierlaiii  et  al.  1966.  Collie  et  al,  1948  and  Lane  and  Saxton  1952). 

Water  is  a pmlicularly  difficult  material  to  work  with.  Because  of  its  great  dielectric  loss  at 
this  frequency  (96.5  GHz)  it  is  practical  to  make  measurements  only  on  very  thin  layers.  It  was 
deterniined  tliat  the  upper  limit  to  the  water  thickness  was  about  0.010  in.  This  limit  is  clearly  set 
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Figure  7.  Polystyrene  nttenuation  curve. 


Figure  8.  //older  arrangement  for  water  samples. 


by  the  available  klystron  power  which  in  this  case  is  about  100  mw,  and  which  suffers  about  a 
4-dB  loss  in  just  getting  to  the  horns.  Varying  thicknesses  of  water  were  obtained  by  using  two 
pieces  of  very  thin  slide  glass  separated  by  pieces  of  shim  stock  of  a known  thickness.  A thin 
pipette  was  used  to  place  the  liquid  between  the  bounding  glass  layers.  Capillary  action  then  drew 
the  water  into  the  enclosed  space.  Figure  8 diagrams  this  arrangement.  With  this  arrangement,  the 
slotted  line  was  adjusted  for  a null  with  the  holder  present  in  the  propagation  path  but  with  no  water 
present.  The  null  shift  caused  by  the  addition  of  water  to  the  holder  was  then  plotted  as  a function 
of  water  thickness. 
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Figure  9.  Phase  sliid  curve  (or  distilled  water. 


Attenuation  measurements  were  made  in  essentially  the  sante  manner,  that  is,  the  sipial  level 
was  recorded  witli  and  witliout  water  in  the  holder. 

Fifruies  9 and  10  are  two  phase  sliilt  curves  ol  distilled  water  and  lake  water,  respectively. 
Tlie  measured  dc  conductivity  ot'these  two  water  samples  was  1.6  pmlio/cm  at  25"C  and  565  /mdio/ 
cm  at  25  C,  respectively.  Fjf;ures  11-13  give  the  attenuation  curves  lor  distilled  water,  rainwater 
and  lake  water.  Tlie  rainwater  sample  had  a measured  conductivity  of  33.6  /imlio/cni  at  25"C,  but 
tliis  sample  was  lost  iKifore  a pliase  shift  curve  for  it  could  be  obtained. 


} 


Considering  tlie  practical  problems  involved  in  working  with  thin  water  samples,  the  agreement 
is  surprisingly  good  lietween  tliese  measurements  and  values  publislied  in  the  literature.  Tlie  major 
sources  of  erre.-  in  this  work  were  in  determining  the  water  thickness  accurately  and  in  reading  the 
"bottom"  of  the  null  shift  when  water  was  introduced  in  the  holder.  The  nulls  tended  to  he  some- 
what flat,  especially  for  thicknesses  greater  than  about  0.006  in.  In  this  case,  the  average  value 
of  two  points  of  ecpial  response  on  either  side  of  the  null  was  taken  to  he  the  true  null  point.  The 
error  bars  shown  in  the  following  figures  were  foiiiid  by  assnniiiig  the  measurement  errors  were 
essentially  random,  and  thus  the  error  shown  is  n^/\N , where  o.  refers  to  the  error  associated  with 
an  individual  iiieasuremeiit  mid  N is  the  nniuber  of  individual  meastirements  made  for  a given  sample 
thickness. 

lahle  111  gives  the  complex  permittivity  values  for  polystyreni"  and  water  as  derived  from  these  | 

measun'iiieiits.  The  values  of  the  appiopriate  inihlislu  d data  mi*  also  listed  for  comparison  purposes.  j 
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Figure  12.  Attenuation  curve  for  r.'iin  water. 
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Pigarc  13.  Attenuation  curve  [or  lake  water. 


3 


i 

I 


1 


COMPLEX  REFRACTIVE  INDEX  OF  ICE  FOG  AT  A RADIO  WAVELENGTH  OF  3 MM 


19 


Table  III.  Complex  permittivity  of  polystyrene  and  water  samples  at  96.5  GHz. 


Material 

f^rad/m) 

a(  neper  s/m) 

k' 

k"(10^) 

Remarks 

Polystyrene 

3275  + 55 

7,25  + 1.35 

2.63  + .09 

11.74  + 2,46 

Perry,  U.T. , 1973 

Di.stilled  de- 
ionized water 

7390  + 1000 

4900  ' 900 

7.5  + 1.54 

17.9  t 6.2 

Perry,  U.T.,  1973: 
o 1.6  /miho  cm  at 
2.5"  C 

Lake  water 

7740  t 1000 

5100  t 900 

8.45  * 1.5 

19.8  t 6 

Perry,  U.T.,  1973: 
o .565  /iinho  'cm  at 
25"C 

Rain  water 

4810  t ,300 

~ 

— 

Perry,  U.T.,  1973, 

<J  33.6  /iinlio  cm  at 
25°C 

Poly.styrene 

- 

- 

2.. 54 

3.05 

Ref  Data  Radio  Engrs 
(ITdtT),  for  25  GIIz 

Polystyrene 

— 

- 

2.55 

0.842 

Ref  Data  Radio  Engrs 
(I’f&T),  for  3 GHz 

Polystyrene 

- 

- 

3.4 

5.78 

Pate  (1967),  94  GIIz 

''Pure"  water 

— 

- 

9.33 

16.3 

Calculated  from  Debye 
etjuations 

Water 

— 

- 

19.2 

29.1 

Collie  et  al.  at 
1.25  cm 

Water 

- 

- 

- 4.44  + .08 

- 

Chamberlain  et  al., 
A .3  mm 

Water 

- 

- 

n^  = 4.30  i .08 

- 

Chamberlain  et  al., 
A 0.118  inm 

Water 

— 

- 

n = 3.55 

- 

Zolotarev  et  al., 
A = 3. 1 mm 

Water 

*- 

- 

13 

23 

Lane  and  Saxton, 

A 6.2  mm  at  20°C 

Water 

— 

- 

7.84 

14.1 

Lane  and  Saxton, 
A = 6.2  mm  at  0°C 

Since  the  measured  ^ is  the  difference  between  the  phase  shift  in  1 m of  air  and  1 m of  the 
sample,  then  the  true  is  obtained  by  adding  2>r/.\^,  or  i?020  rad/m,  to  twice  the  measured  B.  It  " 

is  this  value  that  is  given  in  Table  111.  : 

The  values  foi  /.'  and  k"  measured  here  for  polystyrene  agree  well  with  the  published  values 
for  this  material.  Water  is  well  described  by  a Debye  expression  for  the  complex  permittivity,  at 
least  in  this  part  of  the  millimeter  wave  spectrum.  The  values  of  k'  and  k"  measured  for  distilled, 
deionized  water  agree  well  with  the  Debye  predictions,  and  are  close  to  other  reported  measurements 
in  the  millimeter  wave  spectrum.  From  these  measurements  it  appears  that  the  water  samples  with 
the  higher  amounts  of  impurity  had  a slightly  higher  dielectric  loss  and  dielectric  constant.  How- 
ever, the  error  bars  are  overlapping,  so  one  cannot  draw  any  strong  conclusions. 
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Results  for  ice 

The  icc  samples  used  in  making  these  measurements  were  prepared  from  two  different  water 
samples,  lake  water  and  distilled,  deionized  water.  These  water  samples  had  the  measured  dc 
conductivities  given  in  the  previous  section.  All  ice  samples  were  prepared  in  specially  machined 
aluminum  holvlers  of  prescribed  thickness.  Once  the  samples  were  frozen,  the  holder  sides  could 
be  removed  to  release  them.  Ice  samples  were  made  in  thicknesses  from  0.125  to  0.404  in.  The 
advantage  of  preparing  ice  in  this  manner  is  that  the  frozen  samples  had  flat,  parallel  sides.  In 
order  to  make  ice  samples  as  free  from  internal  stress  cracks  as  possible,  freezing  was  accom- 
plished in  two  stages.  The  water  was  slowly  frozen  over  a 12-hour  period  at  about  -1°C,  then 
transferred  to  the  environmental  chamber  where  it  was  chilled  down  to  about  -28°C  before  the 
dielectric  measurements  were  made. 

Waveguide  was  run  through  the  walls  of  the  environmental  chamber  to  the  interior  horn  antennas 
where  all  measurements  were  made. 

Several  individual  sets  of  measurements  were  made  with  ice  prepared  from  each  type  o'"  water . 
Figures  14-17  are  representative  of  the  phase  shift  curves  obtained  in  this  manner.  The  measured 
phase  shifts  varied  between  about  118  and  128.5  wavelengths/m.  This  phase  shift  is  just  half  the 
true  phase  shift,  as  explained  earlier.  The  ice  made  from  the  lower  impurity  level  water  showed  a 
lower  average  phase  shift  (about  121.3  wavelengths/m),  compared  with  about  129  wavelengths/m 
for  the  other  ice  samples.  Figures  18  and  19  are  attenuation  curves  for  these  ice  types.  All  of 
these  ice  samples  proved  to  be  rather  sensitive  to  the  orientation  angle,  and  it  was  necessary  to 
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Figure  18.  Attenuation  curve,  ice  made  (rom  distilled  water,  T 
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Pignre  19.  Attenuation  curve,  ice  wade  from  lake  water,  T 
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Table  IV.  Measured  complex  permittivity  of  ice  at  96.5  GHz. 

a(  nepers/m ) 


2.27  + 0.8 

3545  t 120 

3.08  + 0.2 

3.99  + 1..50 

nemarKS 

Distilled,  deionized 
water  sample 
-28°C  <T<  -29'^C 

2.32  + 0.8 

3635  + 120 

3.24  + 0.2 

4.17  + 1.69 

Lake  water  sample 
-28..5°C  < r<-3i°c 

average  over  many  data  point,  to  obtain  consistent  results.  The  error  bars  shown  apply  to  all 
points,  and  while  the  points  plotted  are  the  average  of  a'  le.ast  ten  different  nieasurernts  the 
Dos°sihr^  indicative  of  the  uncertainty  in  the  individual  ice  sample  measurements  The 
possible  error  associated  with  the  indh-idual  loss  measurements  was  less  than  .Q  1 dB  The 

attenuation,  in  dB.  was  calculated  from  relative  changes  in  output  meter  level,  and  not  changes  in 
the  precision  attenuator  setting.  i-uanges  in 

in  n-Msuremeiits.  The  fl  appearine 

table  ts  tmce  the  measured  phase  shift.  In  rad.'n  plus  the  pltase.shm  of  free  space.  2n/X.. 
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III.  THE  METEOROLOGY  OF  ICE  FOG 


Introduction 

f OK  is  a common  arctic  siKlit.  Polar  ocean  meteoroloKical  stations  report  a frequency  of 
100  days  per  year  of  foR  occiirrences.  and  Greenland  coastal  reportinR  stations  observe  50  to  60 
days  a year  of  foR  occurrence  (Orvip  1970).  Summer  foRS  in  tlie  Arctic  are  Renerally  caused  by  ad- 
vection  of  relatively  v.arm  and  moist  air  over  meltinR  ice  or  cold  water.  This  foR  is  patchy  and  of 
short  duration,  and  does  not  occur  for  wind  speeds  over  10  m/sec.  DuriiiR  winter,  steam  foR  forms 
over  open  water  leads  in  the  ice  pack  where  the  sea  water  may  be  near  0°C  and  the  air  temperature 
may  l)e  at  -40‘  C,  thus  settir.R  up  a stronp  thermal  pradient  over  a relatively  short  distance.  Such 
a foR  is  Renerally  known  as  an  ice  foR  since  it  is  composed  larpely  of  minute  ice  crystals.  This 
foR,  also  known  as  'arctic  sea  smoke,”  develops  when  very  cold  air  flows  over  open  ocean. 
causinR  a rapid  discharpe  of  moisture  and  heat  to  the  air. 

One  of  the  principal  requirements  for  fop  formation  is  tlie  existence  of  a stable  atmospheric 
cuiditititi  It  is  well  ktijwi.  that  the  mocl  stable  atniosplieTic  coiiJitiuns  occur  unJeT  lUiJ  within 
temperature  inversions.  There  are  three  main  causes  of  temperature  inversions. 

1.  When  a warm  air  mass  flows  over  a cooler  surface,  a temperature  inversion  may  form  in  the 
’jotSAtay  layer  belweei,  tint  iw. . Her.,  the  ttmipt^tatuit  is  relatively  couslam  fiom  dm  Riuunci  up 
to  the  boundary  layer,  then  increases  sharply  upon  enterinp  the  upper  air  mass.  This  type  of 
inversion  is  responsible  for  entrapment  and  intensification  of  the  infamous  Los  Anpeles  smop. 

2.  Wheti  there  i.s  a net  loss  of  enerpy  from  the  eartii  s suriace  (by  outpoinp  lonp  wave  radiation) 
a temperature  inversion  forms  at  pround  level.  These  inversions  become  especially  well  developed 
over  snow  and  ice  surfaces  and  are  quite  common  botli  day  and  nipht  in  the  Arctic  and  Antarctic. 

3.  Subsidence  is  generally  associated  with  hiph  pressure  areas  where  there  is  a descent  and 
warming  of  deep  layers  of  air.  This  effect  causes  a decreasing  of  the  temperature  gradient  and  is 
a frequent  cause  of  temperature  inversions  of  the  lower  troposphere. 

The  type  of  inversion  described  by  2.  is  expecially  common  in  the  Fairbanks  area,  and  has 
been  responsible  for  producing  some  of  the  strongest  temperature  gradients  in  the  world,  with 
values  of  10°  to  30°/100  m (Benson  1965.  1970).  These  inversions  rarely  extend  over  50  or  100  m 
in  height,  hut  within  this  ground-based  inversion  layer  dense  ice  fop  can  form.  Not  only  are  these 
some  of  the  severest  inversions  known,  but  they  are  also  the  most  persistent,  sometimes  lasting  for 
weeks  at  a time  (Ball  1959,  Ber.son  1965.  1970).  Surface  winds  within  these  inversions  are  ideal 
nurseries  for  the  development  and  intensification  of  ice  fogs  where  conditions  are  right  for  their 
production. 

Ice  fog  is  not  restricted  to  coastal  and  ice  pack  areas,  and  commonly  occurs  ss  an  inland  form 
of  low-temperature  air  pollution.  Aside  from  local  sources  of  water  vapor  in  an  inland  area  such  as 
hot  springs  and  caribou  herds,  ice  fog  is  restricted  to  populated  areas.  It  is  produced  by  the  water 
vapor  output  from  automobile  exhausts,  power  plant  stacks,  household  chimneys,  cooling  water  dis- 
charge areas  of  electric  generating  plants,  and  other  sources  associated  with  urban  development. 
Since  air  pollution,  i.e.  the  presence  of  unwanted  material  in  the  air,  varies  in  specific  detail  from 
lafc  luiAtily  to  U)n  ttf^nx  irtt-  fur  ran  qualify  as  si  air  pollullon  n.aterial  ;-r  il.e  Fairbauks  Ft.  I^aln- 
wright  area  of  interior  Alaska.  For  tliis  reason  this  literature,  which  treo  s ice  fog  as  a low  tem- 
perature form  of  air  pollution,  has  been  consulted  for  details  concerning  ice  fog  crystal  measure- 
ments and  the  mechanisms  for  ice  fog  formation. 
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In  arctic  and  subarctic  regions  where  sufficient  volumes  of  moisture  or  open  water  exist,  and 
where  temperatures  go  down  to  about  -30°C,  a fog  composed  of  minute  ice  crystals  appears.  Be- 
cause of  its  close  association  to  human  activities  (such  as  automobiles,  home  heating,  power  plant 
exliausts)  as  well  as  its  obvious  adverse  effects  on  traffic  and  airport  visibility,  Oliver  and  Oliver 
(1949)  suggested  that  ice  fog  l)e  studied  as  an  air  pollution  problem.  Accordingly,  Robinson  et  al. 
(1955),  Benson  (1965,  1970)  and  Ohtake  (1970)  have  shown  that  ice  fog  does  qualify  as  a low 
temperature  form  of  air  pollution  and  have  published  detailed  characteristics  of  the  ice  fog  particles 
and  of  the  meteorological  conditions  occurring  during  ice  fog  events. 


Physical  distinguishing  characteristics  and  optical  properties 

The  optical  properties  and  ice  crystal  types  occurring  in  an  ice  fog  are  different  from  those  of 
other  concentrations  of  ice  crystals  in  the  atmosphere  such  as  cirrus  clouds  or  ’“diamond  dust” 
displays.  Ice  crystals  which  have  well  formed  hexagonal  plates  and  columns  are  called  "diamond 
dust”  because  of  the  way  they  twinkle  under  a low  sun  or  nighttime  illumination.  “Diamond  dust” 
ice  crystals  in  the  atmosphere  are  responsible  for  the  interesting  halo  effects  sometimes  seen  when 
cirrosiratus  clouds  come  between  the  observer  and  the  sun  or  moon.  This  halo  takes  the  form  of 
a luight  ring  having  a radius  of  about  22  degrees,  and  besides  this  halo  there  may  be  numbers  of 
light-bows  and  concentrations  of  light  in  spots,  each  with  a name  of  its  own,  which  combine  to 
form  the  halo  phenomena.  Minnaert  (1954)  describes  in  detail  each  of  these  halo  phenomena,  as 
well  as  the  physical  cause  underlying  the  display.  All  clouds  giving  rise  to  halos  are  composed  of 
small  ice  crystals,  and  the  regularity  of  the  shapes  of  these  crystals  is  responsible  for  the  beautiful 
symmetry  of  the  hght  phenomena.  The  small  halo  phenomenon  is  caused  by  ice  crystals  in  the 
shape  of  a pribfii,  wiTli  the  crystals  generally  between  du  and  300  n in  diameter.  Vertical  pillars  of 
liglit,  appearing  at  night  above  automobile  headlights,  and  above  and  below  the  rising  or  setting  sun, 
are  caused  by  ice  flakes  inclined  nearly  to  the  horizontal  in  the  air.  Crosses,  halos,  light  pillars, 
double  suns,  and  parhelia  or  mock  suns  at  the  same  altitude  as  the  sun  are  all  optical  phenomena 
associated  with  the  passage  of  light  through  ensembles  of  geometrically  configured  euhedral  crys- 
tals. However,  halo  effects  are  entirely  absent  in  ice  fogs.  The  optical  properties  of  an  ice  fog 
are  similar  to  those  of  a water  droplet  fog,  yet  water  droplets  are  not  present  in  ice  fogs  (Robinson, 
Thuman  and  Wiggins  1957).  In  addition,  visibility  is  much  lower  in  an  ice  fog  than  in  a "diamond 
dust'  shower.  These  differences  are  explainable  in  terms  of  the  crystal  sizes  involved.  Smaller 
particles  tend  to  have  a more  pronounced  scattering  effect  than  larger  particles.  The  similarity  in 
the  optical  properties  of  ice  fogs  and  liquid  water  droplet  fogs  is  shown  by  Kumai  and  Russell 
(1969).  This  optical  similarity  between  water  and  ice  fogs  is  due  to  their  optical  values  of  the 
complex  index  of  refraction,  wliich  are  very  nearly  the  same.  They  calculated  the  extinction 
coefficients  of  both  water  and  ice  spheres  (for  an  assumed  particle  density  distribution)  at  wave- 
lengths of  2.2,  2.7,  4.5,  5.75,  9.7  and  10.9  microns.  The  attenuation  coefficients  and  backscattering 
functions  of  ice  fog  and  water  fog  were  within  the  same  order  of  magnitude  for  the  same  fog  con- 
centrations and  wavelengths.  Recent  measurements  of  laser  extinction  in  an  ice  fog  were  made  by 
Munis  and  Delaney  (1972).  They  made  measurements  of  laser  extinction  at  0.6382,  1.15  and  3.39 
microns,  and  compared  these  measurements  with  the  theoretically  calculated  values  given  by  the 
Mie  extinction  coefficients  for  the  measured  particle  distributions  involved.  In  general,  the 
measured  extinction  coefficients  were  within  about  a factor  of  5 of  the  Mie  extinction  coefficients, 
with  the  best  agreement  between  theory  and  experimental  measurements  being  at  1,15  microns. 

Their  experimental  finds  at  0.6328  fi  are  of  special  interest  to  this  research  work  in  that  a helium- 
neon  laser,  operating  at  this  same  wavelength,  was  used  to  obtain  additional  insight  into  the  fog 
statistics. 
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Figure  20.  Small  (J  mm)  snow  crystals. 


Ice  fog  crystal  size  and  formation 

Ice  fog  crystals  range  bciween  2 and  20  /r  in  diameter  (Robinson,  Thuman  and  Wiggins  1957, 
Kumai  1964).  Ice  crystals  larger  than  30  /i  can  occur,  and  in  general  they  are  more  euliedral.  Many 
of  the  small  crystals  tend  to  lack  a definite  crystalline  form.  Because  of  this  Thuman  and  Robinson 
(1954)  coined  the  term  droxtal  in  reference  to  these  ice  fog  particles.  According  to  them,  a droxtal 
is  a small,  poorly  formed  ice  crystal.  It  is  an  equant  solid  particle  with  rudimentary  crystal  faces 
which  seems  to  have  the  characteristics  of  both  droplets  and  crystals.  The  term  droxtal  appearing 
in  the  literature  does  not  mean  that  ice  fog  particles  are  not  crystals.  When  a supercooled  droplet 
freezes,  it  becomes  crystalline  ice  regardless  of  its  shape  and  the  presence  cr  absence  of  any 
crystalline  faces.  There  has  been  no  evidence  of  amorphous  ice  existing  in  ice  fog  particles.  This 
is  not  to  say  that  only  poorly  defined  geometric  forms  exist  for  ice  fog  particles.  Ohtake  (1970)  has 
found  14-  and  20-face  polyhedral  crystals  in  ice  fog  at  -47°C,  although  it  is  true  that  the  vast 
majority  of  the  droxtals  are  “spherical”  ice  crystals. 

Figures  20-22,  from  Benson’s  (1965)  report,  show  the  difference  between  small  1-mm  snow 
crystals,  diamond  dust  crystals  (from  1 to  about  300  g in  diameter),  and  ice  fog  crystals  with 
diameters  of  25  p or  less.  In  general,  ice  fog  particles  are  an  order  of  magnitude  smaller  than 
diamond  dust  crystals  which  in  turn  are  an  order  of  magnitude  smaller  than  snow  crystals.  The 
differences  in  size  and  shape  of  the  ice  fog  droxtals  and  the  cirrus  cloud  or  diamond  dust  crystals 
are  due  to  the  different  rates  of  cooling  and  growth.  Gotaas  and  Benson  (1965)  measured  cooling 
rates  as  high  as  4°C  per  12-hour  period  for  the  arctic  free  atmosphere,  with  typical  cooling  rates  on 
the  order  of  one-half  to  one-third  this  value.  However,  such  low  cooling  rates  in  a very  low  tempera- 
ture environment  cause  water  vapor  to  precipitate  out  of  clear  air  (the  typical  diamond  dust  shower) 
in  the  form  of  well  developed  large  crystals.  Where  the  atmosphere  and  the  water  vapor  contained 
in  it  cool  simultaneously  the  ice  crystals  formed  grow  in  a saturated  air  volume  as  long  as  the  cool- 
ing continues.  The  size  of  the  crystals  increases  with  both  the  original  moisture  content  of  the  air 
and  the  length  of  the  cooling  period  involved,  the  ultimate  result  being  a snowstorm  containing  large 
and  well  formed  ice  crystals.  ■ 


Figure  22.  Ice  fog  crystals  collected  from  a dense  ice  fog  in  a Fairbanks,  Alaska,  supermarket 
parking  lot.  Average  droxtal  diameter  is  in  the  10-20  u range. 


oQ 
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These  conditions  are  in  sharp  contrast  to  the  formation  of  ice  fog  crystals.  When  droxtals  form 
from  coolinR  exhaust  {rases.  ori{rinally  at  temperatures  exceedinp  100°C,  the  cooling  takes  place  in 
less  than  10  seconds  when  the  gases  are  ejected  into  air  of  temperature  -35°C  or  less  (Benson 
1965).  These  cooling  rates  are  on  the  order  of  10®  times  greater  than  those  in  the  free  atmosphere. 
Cooling  tower  fogs  and  steam  fog  from  open  areas  of  water  also  have  rates  of  cooling  that  are 
orders  of  magnitude  greater  than  in  the  free  atmosphere.  Such  a rapid  cooling  condenses  the  water 
vapor  into  many  small  droplets  which  become  supercooled,  and  at  temperatures  below  about  -SO^C 
they  begin  to  freeze.  This  process  also  differs  from  that  in  the  free  atmosphere  in  that  the  crystals 
do  not  grow  in  a continuously  saturated  environment.  When  studied  under  the  microscope,  many  ice 
fog  particles  appear  to  be  twinned  (two  crystals  joined  along  a common  plane).  Collections  of  ice 
fog  crystals  show  many  twinned  crystals,  whereas  twinning  of  cirrus  cloud  or  diamond  dust  crystals 
and  snow  crystals  is  not  prevalent. 

Ice  crystal  density  distribution 

Thuman  and  Robinson  (1954)  made  some  of  the  first  studies  of  ice  fog  size  density  distributions 
at  Eielson  Air  Force  Base  (Alaska).  They  found  that  the  average  particle  size  in  ice  fog  at  -40°C 
was  13  n diameter,  and  the  size  increased  to  16  g at  -30°C.  Kumai  (1964)  showed  some  size  dis- 
tributions for  ice  fog  particles  with  sizes  from  4 to  12  g for  -30°  and  -41°C  in  downtown  Fairbanks 
and  at  the  airport,  respectively.  Kumai's  distributions  were  based  upon  a count  of  the  crystals 
precipitated  on  a collection  surface  (a  glass  slide)  and  are  subject  to  some  doubt  since  the  fall 
speed  of  the  ice  fog  particles  may  be  less  than  local  turbulent  air  updrafts.  Ohtake  (1970)  made 
calculations  of  the  terminal  velocity  versus  particle  size  using  Stokes’  law  for  spherical  particles. 
For  33-g  particles,  the  terminal  velocity  is  3.47  cm/sec,  while  for  9-/i  particles  the  terminal  velocity 
is  only  2.7  mm/sec.  Thus  only  slight  air  motion  ciui  cause  considerable  uncertainty  in  any  statis- 
tical properties  deduced  from  ice  fog  precipitation  studies. 

For  this  reason,  Ohtake  made  studies  of  ice  fog  density  distributions  using  an  impaction  device. 
This  device  consisted  of  a cylinder  with  a nozzle  at  one  end  facing  an  interior  glass  slide  coated 
with  a thin  film  such  as  silicone  oil.  A measured  volume  of  air  was  drawn  through  the  nozzle,  and 
the  ice  crystals  collided  with  and  were  captured  on  the  silicone  oil  film.  The  slide  could  then  be 
withdrawn  for  microscopic  study  and  photographing.  Using  this  technique,  samples  of  ice  fog  were 
taken  and  studied  over  tU-  li-mperaturw  rang«  30°  to  - 41°C.  AlUrough  his  data  were  swutwhat 
scattered,  it  was  clear  tlsat  the  mean  diameter  of  the  ice  fog  particles  was  a linearly  increasing 
function  of  temperature.  Figure  23,  taken  from  Ohtake,  illustrates  this  relationship  between  mean 
diameter  and  temperature.  From  this  figure  it  is  clear  that  30-g  size  particles  are  the  rule  for 
-31°C  ice  fogs,  with  particle  sizes  of  10  g or  less  dominating  in  -40°C  ice  fogs.  Ohtake  found  that 
the  lower  the  temperature  became,  the  narrower  was  the  breadth  of  the  distribution  in  crystal  sizes. 

In  additiv/Ti,  fot  teiiiptratufeS  of  about  -3T  0 and  below,  the  spherical  particles  tlhuman  and 
Robinson’s  ’’droxtals”)  make  up  the  majority  of  the  fog  particles,  whereas  at  higher  temperatures 
other  crystalline  forms  are  more  common. 

The  role  of  nuclei  and  low  temperatures  in  crystal  formation 

Ice  fog  is  a phenomenon  associated  only  with  subfreezing  temperatures.  Naturally  occurring 
ice  fog  has  always  been  observed  to  disappear  when  the  temperature  rose  above  -30°C  (Benson 
1965).  The  occurrence  of  liquid  water  droplets  in  subfreezing  temperatures  is  the  commonly  ob- 
served phenomenon  of  supercooling  which  occurs  in  fogs  above  about  -30°C.  The  only  condition 
necessary  for  producing  supercooled  water  is  that  no  ice  be  present  vben  the  water  is  cooled  down 
to  and  below  0°C.  That  is.  the  melt  must  not  be  "seeded”  with  the  solid  phase.  Furthermore, 
supercooling  always  takes  place  in  unseeded  melts  before  freezing  occurs.  This  is  also  true  for 
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Figure  23.  Relationship  between  mean  diameters  o( 
total  ice-fog  crystals  and  temperatures.  (From  Ohtake 
1970.) 


supersaturated  solutions  prior  to  crystallization  and  precipitation  of  one  of  the  components  of  the 
solution.  In  a cloud,  freezing  starts  either  when  the  cloud  is  cooled  below  a certain  temperature  or 
when  the  cloud  is  seeded  or  nucleated  with  some  suitable  foreign  substance. 


Many  observers  have  reported  that  at  a temperature  of  -40  i 1.5°C  clouds  are  always  composed 
of  ice  crystals  (Civilong  1945,  1947a,  b,  Schaefer  1946,  1948,  D’Albe  1948.  Oliver  and  Oliver  1949, 
Mason  1950.  Civilong  used  a Wilson  cloud  chamber  in  his  study,  and  he  found  that  when  the  tem- 
perature was  greater  than  -35°C  at  the  end  of  expansion  progress,  when  condensation  had  just 
started  to  occur,  no  ice  crystals  were  to  be  observed.  However,  when  the  temperature  was  below 
-35°C,  then  a few  ice  crystals  formed  in  the  fog.  Larger  numbers  of  ice  crystals  began  to  appear 
as  the  temperature  was  lowered,  and  at  about  -41°C  the  fog  was  composed  entirely  of  ice  crystals. 
Civilong  found  that  below  this  temperature,  -41°C,  sublimation  was  involved  where  ice  crystals 
and  no  water  droplets  formed.  This  “threshold  temperature"  was  found  to  increase  with  an  in- 
creasing nucleating  impurity  added  to  the  air. 


Schaefer  (1946,  1948)  made  supercooled  fog  by  blowing  moist  air  into  a cold  chamber.  He  found 
that  ice  crystals  did  not  form  until  the  temperature  was  lowered  below  about  -39°  C.  or  upon  the 
introduction  of  an  appropriate  nucleating  substance.  Vonnegut,  in  a series  of  papers  (1947,  1949a,  b), 
investigated  the  effects  of  silver  iodide,  potassium  iodide,  and  several  lead  oxide  compounds  as 
nucleating  agents  for  ice  crystal  formation  in  clouds.  He  found  silver  iodide  (Agl)  to  be  the  most  ef- 
fective nucleating  agent,  and  attributed  its  effect  to  its  close  crystalline  similarity  with  ice  crystals. 
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thus  allowing  water  vapor  in  the  air  to  sublimate  directly  onto  tlie  silver  iodide  crystals  In  urban 
areas,  combustion  byproducts  of  sizes  generally  less  than  about  0.1  p diameter  can  alL'be  e ei 
tive  ice  fog  nucleating  agents  (Otitake  1967,  Kumai  and  O’Brien  1965)  The  role  tint  foreign 
tnipunt.es.  of  all  types,  play  as  nucleating  agents  in  ’’clean”  air  is  well  rocunmuted  eTneciallv 

TraTth^’""'.”*^  the  electron  microscope  (Ohtec  1949,  Mason  1950.  Kuroiwa  1957.  KunL  1950 
tliese  condensation  nuclei”  play  an  important  role  in  ice  fog  formation  is  undoubted  since’ 

area  o Al  itt,  a ’ .-.■porta  that  In  tlte  Faltl,anks-Ft.  Wainwrlpl.t 

area  of  Alaska  ice  fog  begins  to  set  in  at  aliout  -35  C and  is  always  well  developed  at  -40°C  Mnd 

30  , and  the  pollutants  need  be  of  no  special  type  in  order  to  be  effective. 

fott„au‘l„° ' m,ll!."r'  “I  ..“1...  a.-e  todiapenaable  to  Icc  fop 

““  '■'>«  l'.«l  .'itclel.  01, take 

me  temperature  is  -37  C or  lower.  Ut.uke  (1970)  also  reported  that  71  out  of  713  crystals  from 
hena  and  Manley  Hot  Springs  inspected  under  the  electron  microscope  did  not  have  any  nuclei 

t.rks  1 temperatures  below  about  -40»C.  Furthermore,  at  the  center  of  Fair- 

banks  1.7.0  of  236  crystals  examined  from  an  ice  fog  at  -39'^C  had  no  nuclei. 

Phase  change  effects  on  development  of  ice  fog 

fiir  particles  is  a linearly  increasing  function  of  atmospheric  tempera- 

ruo  r re'r?““T  ‘t  P“'  '■>'  molehtre  vX  InTcterl 

a atmosphere.  Typical  coricentralions  for  Fairbanks,  as  reported  hv  Ohtake  (1970)  are 

rToTl,  ' ■ s o.'ly  abon,  300™,densa. 

that  concp*  ‘'f  C'>ena  Hot  Springs  under  ice  fog  conditions,  he  found  about  twice 

at  concentration  of  ice  crystals.  Other  observers  (Kumai  1964,  Benson  1965)  are  siibst  iiitiallv  in 
agreement  with  these  measured  concentrations.  simstaiit tally  in 

and  tlreffecrof  th?!"'^^  ‘■^"'“spl.eric  temperature 

en  observed  that  air  containing  high  concentrations  of  combustion  byproducts  can  furnish 

-SrorS'c ' tTc"?  ''"'f temperatures  as  high  as 
sfl  r~r  k ■ to  ice  crystals  at  this  temperature  has  a 

droplet  fogs  are  forming,  such  as  over  open  water  areas,  then  tlie  low  level  air  will  be 
saturated  or  close  to  saturation  with  respect  to  water. 

For  example.  -34‘-C.  the  capacity  of  the  air  to  hold  water  vapor  decreases  with  decreasing 
temperature  at  about^-0.0279  g/m’*’  °C  (over  water).  Thus,  if  the  air  is  saturated,  cooling  the  L 

sunnose  tf  ^ r”  ^«"densed  moisture.  Now 

trhoirn  r additional  reduction  in  the  .air’s  capability 

o hold  moisture  because  tlie  saturation  vapor  pressure  must  be  found  with  respect  to  ice  The  ' 

o^  ro  ri38?/m3  tT"  T."  density  is 

n n^4  e/n?3^  ^'i  ''"P”  P'^d^^d  simply  by  freezing  droplets  at  -34°C  is 

h k o no,n  -34^C  to  -35"C.  whereas  only 

about  0.0279  g/m  water  was  condensed  in  going  from  -33°  to  -34°C.  Tims  the  excess  moisture 

TooHng  at  ^3? c"' N 3.16  times  as  great  as  the  moisture  condensed  by  one  degree  of 

cooling  at  34  C.  Now  suppose  that  freezing  is  postponed  (for  whatever  reason)  until  -44°C  Tlie 
excess  moisture  produced  by  freezing  at  -44°C  is  now  3.36  times  as  much  as  is  produced  by  ’one 

gree  of  coohng  at  -44  C.  Therefore  the  relative  contribution  of  water  vapor  due  to  freezing 
increases  with  ii  decreasing  freezing  point. 
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Uiis  .iiguiiieiit  is  mount  to  (‘tiiphasizo  and  siibstantiato  two  ohservational  tacts  about  ice  I'oks’ 
tliey  liavo  been  observed  to  appear  suddenly  over  wid(‘spread  areas,  and  denser  logs  liave  been 
observed  at  low  temperatitres  thati  at  bigber  temperatures.  The  above  calculations  were  based  up- 
on the  densities  of  water  vapor  at  saturation  over  water  and  over  ice  (List  1966  Tables  108  and 
109). 


IV.  DIELECTRIC  BEHAVIOR  OF  A MEDIUM  WITH 
HOMOGENEOUSLY  DISTRIBUTED  PARTICLES 

Dielectric  behavior  of  a medium  of  small  particles 

Van  (le  lliilst  (1957)  derived  an  expression  for  tbe  complex  refractive  index  of  a inedinm  coii- 
taimtig  many  stnall  scatteriMs.  His  expression  for  m n'  In"  is 

in  ' 


n'  1 . Ini  1S(0)| 


n"  RelS(0)| 


u,  \ 

seatterer  tumiliet  detisity 

cotnplex  atnplitnde  coelficient  for  forward  scattering 

1 , ixa;  . !,„») 

and  where  a|;  and  are  tbe  complex  atnplitnde  coefficients  for  the  nth  magnetic  and  electric  tnodes 
respectively. 

This  expression  was  derived  to  appi  to  the  case  wbcre  tbe  distance  between  particles  is 
large  compared  with  both  the  wavelength  and  tlie  size  of  the  particles.  Dense  packing  of  particles 
is  thus  excluded.  Scattering  eoutritiutions  from  the  small  particles  which  affect  the  phase  of  the 
transmitted  wave  come  not  trom  particles  on  all  sides,  but  only  from  particles  that  are  nearly 
ditectly  behind  the  point  ot  observation.  This  expression  is  accurate  for  the  complex  refractive 
index  of  the  medium  close  to  1, 

An  alternate  expression  for  the  opposite  case  is  the  e.\  *ssion  derived  by  Lorentz-Lorenz 
(Jackson  1963): 


where  k 
N 
S(0) 


4rraN 


I 2 

m ^ 


where  a - polarization  of  the  particle 

p oEq,  the  induced  dipole  moment  due  to  the  incident  electric  field  E^. 

This  expression,  derived  as  a part  of  the  theory  of  molecular  optics,  is  clearly  not  the  same  as 
tbe  Van  de  Hulst  expression,  yet  both  expressions  give  the  compl  r refractive  index  of  a medium  as 
a whole  in  terms  of  the  scattering  properties  of  the  individual  particles. 
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The  Lorentz-Lorenz  expression  is  derived  for  the  conditions  tliat  tlie  mutual  distances  of  the 
particles  are  small  compared  witli  the  wavelength,  which  of  course  means  that  the  particles  are 
miicli  smaller  than  a wavelengtli.  The  field  at  a given  point  of  observation  is  due  to  the  scattered 
fields  of  a fairly  large  number  of  nearby  particles  on  all  sides  of  the  point  of  interest.  These  fields 
are  nearly  in  phase  with  the  incident  field,  and  it  is  the  statistical  effect  of  these  fields  which 
gives  rise  to  tlie  Lorentz-Lorenz  formula. 

However,  tlie  Lorentz-Lorenz  formula  can  give  the  Van  de  Hulst  expression  for  tlie  case  of  a 
dilute  gas  of  small  scatlerers  where  . 1.  Here, 


^nuN 


- 1) 
ml  f 2 


3l2(/«,„  1)1 

2(m,„  - 1) 


• 1 ^ 2nuN 

and  refers  to  the  complex  dielectric  constant  of  the  medium  as  a whole.  For  scattereis  small 
in  terms  of  a wavelength 

S(0)  ~ jk^u 


since  for  spheres  of  radius  a 


where  m is  the  complex  refractive  index  of  the  sphere  material.  So 


is  the  coefficient  for  the  first  electric  mode  of  a dielectric  sphere,  which  is  just  the  familiar 
expression  for  the  case  of  Rayleigh  scattering.  Thus. 

tn„,  l+2nNS(0) 

which  is  just  Van  de  Hulst ’s  expression. 

This  expression  is  appropriate  to  use  to  model  the  complex  refractive  index  of  an  ice  fog 
since  an  ice  fog  is  comiiosed  of  a dilute  assemblage  of  ice  particles  that  are  small  in  terms  of  the 
incident  electromagnetic  wave. 

Sample  calculations  using  the  Van  de  Hulst  expression  for  tlie  index  of  refraction  of  the  ice 
fog  medium  are  delayed  to  a later  section  pending  the  fog  impactor  measurements  and  a discussion 
of  the  results.  However,  it  is  clear  that  this  expression  predicts  a refractive  index  which  is 
sensitive  to  the  total  volume  of  ice  per  unit  volume.  Thus,  on  a per  ice  sphere  basis,  an  icte 
particle  30  g in  diameter  has  nine  times  the  effect  of  one  10  g in  diameter,  A relatively  small 
number  of  large  particles  can  strongly  influence  the  measured  refractive  index.  Thus,  in  order  to 
accurately  apply  Van  de  Hulst 's  formula  you  must  know  not  only  the  numbers  involved  hut  also  the 
particle  sizes. 
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Radio  wave  energy  loss  mechanisms 

When  a plane  electroiiiaKnetic  wave  is  incident  upon  a dilute  cloud  of  small  particles,  energy 
is  extracted  from  the  incident  wave  in  proportion  to  the  absorption  and  scattering  cross  sections 
presented  by  the  small  particles.  For  the  case  where  the  incident  electromagnetic  wave  has  a 
wavelengtli  much  larger  than  the  circumferential  size  of  the  aerosol  particles,  the  Rayleigh  approxi- 
mations for  the  scattering  and  absorption  cross  sections  (Stratton  1941)  are  valid  and  are  given  by 
tlie  expressions 

Im  ( X) 

where  Im  imaginary  part  of  argument 

0^  particle  diameter  in  wavelengths 

f 2 

and  m is  the  complex  refractive  index  of  the  aerosol  dielectric  material. 

A computer  program  was  used  to  determine  the  extinction,  scattering  and  absorption  cross 
sections  for  ice  spheres  from  1 to  35  fx  in  size.  This  computer  routine  was  used  since  it  had 
already  been  written  for  use  at  optical  frequencies  to  calculate  6328  A helium-neon  laser  extinction 
due  to  a fog.  This  routine  calculates  the  complex  amplitudes  of  the  nth  electric  and  magnetic  modes 
set  up  in  eacli  dielectric  sphere  by  the  incident  wave.  These  amplitude  coefficients,  a and  b 
(Stratton  1941),  were  in  turn  used  to  calculate  the  following  parameters  appropriate  for  each  fog 
particle  size  from  the  following  expressions: 

Scattering  cross  section 

Q _ total  power  scattered  over  An  .steradians 
® incident  power  density 


n~l 

Extinction  cross  section 

total  power  extracted  from  the  wave 


Qv 


incident  power  density 


Re 

2rr 


2 


n=l 


(1) 


(2) 


Absorption  cross  section 


(3) 
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All  equivalent  expression  for  tlie  extinction  cross  section  is  f;iven  liy 
Qg  2kn" 

whore  n"  is  the  imaginary  part  of  the  Van  de  Huls.  expression  for  the  complex  index  of  refraction 
given  earlier,  and  k 2n  'A. 

While  the  Van  de  Hulst  expression  gives  an  extinction  coefficient,  it  does  not  give  separate 
expressions  for  scattering  and  absorption.  Since  these  quantities  may  be  of  interest,  expressions 
1-3  were  used  to  deterniine  these  parameters  with  the  aid  of  a previously  written  computer  program. 
In  computing  the  scattering,  absorption  and  extinction  cross  sections  in  the  manner  defined  above 
tive  terms  were  taken  in  the  summation  indicated.  However,  this  is  not  necessary  as  all  terms 
greater  than  the  tirst  are  orders  of  magnitude  smaller. 

Table  V lists  the  dominant  Rayleigh  parameters  and  the  attenuation  and  scattering  cross  sec- 
tions tor  ice  fog  particles  at  96.5  GHz  for  a complex  dielectric  consiani  of  3.08  - ;(0.004)  which  is 
appropriate  for  the  actual  measured  complex  dielectric  eoiistani  at  iliis  frequency.  It  is  clear  that 
at  ,3  mm  wavelengtli  absorption  is  the  dominant  extinction  meclianism. 

For  a wave,  passing  through  such  a dissipative  medium,  Hie  extinction  coefficient  per  meter 
thickness  of  the  fog,  r,  is  found  from 


E 

all 


nD"^ 


!^(D)Q^(m,D)  meters 


-1 


particle.s 


where  D particle  diameter,  meters 

N(D)  number  of  particles  of  diameter  D per  cubic  meter 
D)  = normalized  extinction  coefficient 

It  is  clear  that  this  implies  that 


P . P,  e- 

where  P received  power  density 
Pg  initial  power  density 
X absorption  path  length. 

The  application  of  these  expressions  to  the  electromagnetic  measurements  is  deferred  to  a 
later  chapter  pending  the  actual  fog  impactor  and  electromagnetic  measurements.  However,  it  can 
be  said  here  that  this  model  is  in  reasonably  good  accord  with  the  EM  measurements  considering 
the  uncertainty  which  is  shown  to  exist  in  the  fog  impactor  measurements. 


COMPLEX  REFRACTIVE  INDEX  OF  ICE  FOG  AT  A RADIO  WAVELENGTH  OF  3 MM  35 


Table  V.  Complex  scattering  coefficients  and  attenuation  cross 
sections  for  ice  particles  at  96.5  GHz. 

Particle  Complex  first  electric  Absorption  Scattering  Normalized  extinction 


diani 

nmde^aett 

cross  section 

cross  section 

cross  section 

(P) 

“ '>1 

re 

'V 

im 

(m‘^1 

lQp/lrrl}~''l)] 

1 

3.25  x 

lO'*'* 

2.87  < 

10  »° 

1.48  X 10  *° 

• >r. 

3.75  X 10 

1.88  X 10  ® 

o 

2.61  V 

10- 

2.30  X 

10-° 

1.19  X 10  17 

2.40  > 10  °° 

3.78  • 10  ® 

3 

8.80  X 

10 

7.75  X 

10-° 

4.01  V 10-^7 

2.73  .<  10 

5.67  X iQ-** 

■1 

2.09  :< 

10' 

1.84  X 

10"° 

9.. 50  10  *7 

1.54  >.  10-^* 

7.56  / 10-® 

r> 

4.07  X 

10-11 

3.59  V 

10  ° 

1.86  :•  10  *® 

5.86  > 10-^1 

9.46  X 10'® 

() 

7.04  X 

10  11 

6.20  X 

10-° 

3.21  X 10-^® 

1.75  - 10-°° 

1.13  X 1Q-® 

7 

1.12  X 

10- in 

9.84  X 

10-° 

5.09  < IQ-*® 

4.41  X 10 

1.32  X 10~® 

8 

1.67  X 

iQ-in 

1.47  X 

10  ’ 

7.60  X 10"*® 

9.83  ■ 10-°° 

1.51  X 10"® 

il 

2.38  X 

lo-in 

2.09  V 

10-' 

1.08  X 10  *® 

1.99  X 10-1° 

1.70  X 10"® 

10 

3.26  X 

lo-in 

2.87  X 

10-7 

1.49  X 10-1® 

3.75  x 10-1° 

1.89  X 10  ® 

11 

1.34  X 

lo-in 

3.82  X 

10-7 

1.98  X IQ-^® 

6.65  X 10-1° 

2.08  X 10-® 

12 

5.63  X 

10-  in 

4.96  X 

10-7 

2.. 57  X IQ-*® 

1. 12  > 10-*° 

2,27  X 10“® 

13 

7.16  X 

10-1° 

6.30  X 

10-7 

3.26  X 10-^® 

1.81  --  10-*° 

2.46  10-® 

H 

8.95  -• 

10-1° 

7.87  X 

10-7 

4.08  10  *® 

2.82  X 10-*° 

2.65  X 10"® 

ir> 

1.10  > 

10-° 

9.68  X 

10-7 

.5.01  X 10  '® 

4.27  10  *° 

2.84  -■  10"® 

16 

1.34  < 

10-° 

1.18  X 

10  ° 

6.08  X 10"^® 

6.29  • 10-*° 

3.03  >.  10"® 

17 

1.60  X 

10-° 

1.4  1 X 

10  ° 

7.30  X 10'^® 

9.06  < 10"*° 

3.22  X 10-® 

18 

1.90  X 

10  ° 

1.67  X 

10*° 

8.66  X 10"'® 

1.28  X U)-*7 

3.41  X 10'® 

19 

2.24  • 

10  ° 

1.97  X 

lQ-° 

1.02  X 10-'“’ 

1.77  X 10-*7 

3.60  X 10"® 

20 

2.61  X 

10“° 

2.. 30  V 

10  ° 

1.19  X IQ-^'l 

2.40  X 10  *7 

3.79  X 10"® 

21 

3.02  X 

10-° 

2.66  < 

10  ° 

1.38  X lO-*"* 

3.22  X l0-*7 

3.98  X 10"® 

22 

3.48  X 

lQ-° 

3.05  < 

10  ° 

1.58  > 10- 1"! 

4.25  X 10-*7 

4 17  X 10-® 

23 

3.98  X 

10-° 

3.19  • 

10  n 

1.81  X IQ-*'’ 

5.55  < 10- *7 

4.36  >.  10-® 

24 

4.52  > 

10° 

3.97  •- 

10  ° 

2.05  X lO-^"! 

7.17  X 10-*7 

4.56  X 10~® 

2.5 

5.11  X 

10"° 

4.48  X 

10"° 

2.. 32  < lO"''’ 

9.16  X in-*7 

4.75  X 1Q-® 

26 

5.75  X 

U)-° 

5.04  X 

10-® 

2.61  X IQ-l'* 

’.16  ■ U)-*® 

4.94  X 10-® 

27 

6.45  X 

lQ-° 

5.65  X 

lO'® 

2.93  X U ’ ' 

1.45  •:  10  *® 

5.13  X 10-® 

28 

7.19  X 

10  ° 

6.30  X 

1Q-® 

.3.26  X IQ-*'! 

1.81  X 10  *® 

5.33  X 1Q-® 

29 

8.00  y 

10-° 

7.00  X 

1Q-® 

3.62  X 10 

...23  X 10-*® 

5.  .52  X 1Q-® 

30 

8.86  X 

lQ-° 

7.75  X 

IQ-® 

4.01  X 10-*“* 

2.74  X 10-*® 

5.72  X 10“® 

31 

9.78  X 

10-° 

8.55  X 

10“® 

4.43  X 10"^® 

3.33  X 10-*® 

5.91  X 1Q-® 

32 

1.08  X 

10-° 

9.40  X 

IQ-® 

4.87  X 10-^“* 

4.03  X 10" *® 

6.11  X 10~® 

33 

1.18  X 

10-° 

1.03  X 

IQ-® 

5.. 34  X IQ-^'^ 

4.85  X 10-*® 

6.30  X 1Q-® 

34 

1.29  X 

10-° 

1.13  X 

IQ-® 

5.84  X IQ-*'* 

5.80  X IQ-*® 

6.50  X 1Q-® 

35 

1.41  X 

10-° 

1.23  X 

10"® 

6.37  X IQ-*'* 

6.90  X IQ-*® 

6.70  X 10  ® 
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V.  METHOD  FOR  MAKING  COMPLEX  PERMITTIVITY  MEASUREMENTS  OF  A FOG 

Literature  review  of  millimeter  wave  Fabry-Perot  interferometry 

If  a cavity  is  used  to  measure  low  losses  and  refractive  index  cliaiif’es  of  ten  A/-units  or  less 
(1  A/ -unit  (n  - 1)  ■ 10®,  where  n is  the  index  of  refraction),  it  must  liave  a high  Q.  Using  a cavity 
with  a Q of  10®  o.ie  can  easily  detect  a 1 N-unit  change  (97  kHz  at  97  GHz),  One  type  of  high  Q 
cavity  successfully  used  at  millimeter  wave  frequencies  is  the  spherical  mirror  Fabry-Perot 
interferometer  (FPS),  The  FPS  is  widely  used  in  maser  and  laser  teclinology  because  of  its  inlier- 
eiitly  high  resolving  power  and  its  relative  ease  of  adjustment.  This  cavity  has  a much  higher  Q 
than  the  quarter-wave  coaxial  cavity  more  commonly  used  in  the  microwave  spectrum  (which  may 
have  unloaded  Q’s  of  up  to  10,000).  An  extensive  body  of  literature  exists  on  the  Fabry-Perot 
cavity,  and  some  good  tutorial  papers  are  those  of  Fox  and  Li  (1961),  Schawlow  and  Townes  (1958), 
Boyd  and  Gordon  (1961),  and  Boyd  and  Kogelnik  (1962),  The  paper  liy  Kogelnik  and  Li  (1066)  is 
an  excellent  review  of  open-sid^-d  tcsoruUots  and  contains  niurdi  maiiirial  of  a tutorial  natu'-p. 
Culshaw,  in  a series  of  papers  (1959,  1960,  1961,  1962),  discussed  the  application  and  design  of 
Fabry-Perot  interferometers  for  the  millimeter  wave  spectrum.  Mucli  of  Culshaw's  emphasis  was 
placed  upon  the  tlieory  and  design  of  the  parallel-plate  Fabry-Perot,  and  methods  of  coupling  the 
electromagnetic  energy  into  and  oirt  of  this  type  of  structure.  His  1960  paper  describes  one  such 
plane-parallel  resonator  built  for  operation  at  a wavelength  of  6 mm  which  had  a Q of  about  100,000. 
Zimmerer,  Anderson,  Strine  and  Beers  (1963)  discussed  the  construction  of  both  plane-parallel  and 
splierical  Fabry-Perot  resonators  and  considered  the  problem  of  obtaining  efficient  power  transfer 
into  tlie  resonant  structures.  A separate  paper  by  Zimmerer  (1963)  considered  measurements  of  the 
transverse  E-field  variations  of  a spherical  mirror  Fabry-Perot  resonator  as  a function  of  reflector 
size  and  spacing.  Valkenburg  and  Derr  (1966),  Culshaw  and  Anderson  (1962).  Frenkel  and  Woods 
(1966),  Degenford  and  Coleman  (1966),  and  French  and  Arnold  (1967)  are  among  some  who  have 
built  and  used  the  FPS  to  make  measurements  of  gas  absorption  and  permittivity  in  the  millimeter 
wave  spectrum.  Dees  and  Slieppard  (1965)  reported  on  the  construction  of  a semi-confocal  FPS  with 
a Q of  88,000  at  168  GHz.  Zimmerer  (1962)  reports  the  construction  of  a novel  high  precision  wave- 
meter  covering  the  50  to  75  GHz  band.  Zimmerer  mounted  his  conlocal  resonator  on  a precision 
screw  thread,  and  took  advantage  of  the  fact  that  cavity  resonances  occur  at  intervals  of  A/2  to 
measure  the  operating  frequency  to  a high  degree  of  accuracy. 

There  are  two  basic  geometric  forms  for  the  Fabry-Perot  interferometer;  the  parallel  flat  plate 

Srid  tlw  Ttiii-  pntallrd-plab  is  tlie  form  of  the  original  oj  tical  resonator 

as  constructed  by  Fabry  and  Perot,  and  is  perhaps  the  most  widely  known  version  of  the  two.  Tliis 
type  uses  very  highly  reflective  plane  parallel  plates  whose  surfaces  iiiusl  both  he  flat  to  wilWii 
about  1/lOOOA  and  parallel  to  a very  high  degree  to  minimize  diffraction  losses  in  order  to  attain  a 
high  operating  Q.  The  millimeter  analog  of  this  type  has  been  built  (Culshaw  1961,  for  example), 
but  plate  alignment  is  a nontrivial  problem  and  in  general  relatively  large  plates  are  required  to  re- 
duce diffraction  losses,  even  for  relatively  close  plate  separations.  Both  Culshaw  in  his  series  of 
papers  (op.  cit.)  and  Zimmerer  et  al.  (1963)  have  treated  the  problem  of  coupling  energy  into  and  out 
of  these  flat-plate  resonators. 

The  confocal  type  of  the  Fabry-Perot  interferometer  has  all  of  the  advantages  and  few  of  the 
disadvantages  of  the  flat  plate  type  of  interferometer.  The  spherical  mirrors  eliminate  the  critical 
alignment  problem,  and  the  focusing  action  of  these  mirrors  decreases  the  diffraction  loss  signifi- 
cantly. The  confocal  FPS  has  two  mirrors  of  equal  radius  of  curvature,  with  the  radius  of  curvature 
of  ^dCii  TuiTiui  coi.icitfetit  With  Hit  smface  jf  the  dhct  (set  Fig.  21a).  i<oolo«’al  FPS  la  a 

version  of  the  confocal  type,  with  one  spherical  mirror  imaged  by  the  flat  plate.  The  semi-confocal 
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Figure  24.  Confocal  and  semi-confocal  Fabry-Perol  inlerferomeiers. 


FPS  is  quite  practical  and  easily  built,  and  the  coupling  problem  is  simplified  in  that  waveguide 
feeds  may  be  brouglit  tlirough  tlie  back  of  the  flat  and  located  symmetrically  about  the  center  line 
of  the  cavity.  Tlie  separation  of  the  flat  and  the  spherical  mirror  is  at  a distance  of  lialf  the  radius 
of  curvature  of  tlie  spherical  mirror  (Fig.  24b).  The  field  and  mode  distribution  for  semi-confocal 
and  confocal  FPS  are  identical  for  equal  mirror  sizes  and  equal  radii  of  curvatu-'e.  However,  the 
Q of  a semi-confocal  FPS  is  lialf  that  of  the  confocal  FPS  since  energy  has  to  travel  only  lialf  as 
far  in  the  semi-confocal  case. 

Theoretical  considerations 

In  order  to  discuss  the  relationship  between  cavity  Q,  mode  structure  and  mirror  size,  it  is 
necessary  to  discuss  briefly  the  basic  theoretical  development  behind  these  resonators.  In  this 
sense,  the  approach  of  Fox  and  Li  (1961)  will  be  followed.  Their  approach  to  the  open  resonator 
problem  was  to  consider  a propagating  wave  which  was  reflected  back  and  forth  between  two  plane 
mirrors.  By  assuming  an  arbitrary  initial  field  disUibution  at  the  first  mirror,  the  field  at  the 
second  mirror  was  calculated  using  the  scalar  formulation  of  Huygen’s  principle.  This  field  was 
used  to  calculate  the  field  back  at  the  fL-st  mirror  on  the  second  transit,  etc.  Only  after  several 
hundred  iterations  did  the  field  approach  a steady  state,  differing  only  by  a complex  constant  from 
mirror  to  mirror.  Without  further  derivation,  their  iterative  equation  for  confocal  spherical'mirrors  of 
circular  cross  section  is 


a 2jt 


rjdi^jdrj 
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with 


Q - iirjrg  cos  (</.! 

I 1.2 

_ — , I 1,  2 for  larue  h/n. 

2b 

UsiriK  this  simplification,  ^o)  simplifies  to 

0 0 

which  is  strictly  valid  for  (a^/b,\)  « (b/a)^. 

FiRure  2r>  defines  the  terms  appropriate  for  the  above  terminoloRy.  In  this  equation,  the  single 
most  important  parameter  is  (a~/bA).  This  parameter,  called  the  Fresnel  number,  is  approximately 
the  nuniher  of  Fresnel  zones  seen  in  one  mirror  from  the  center  of  the  other,  and  determines  the 
numbers  of  ripples  in  the  amplitude  luid  phase  distribution  of  the  electric  field  of  eacli  mirror. 

Calculations  oy  Fox  and  Li  made  using  the  iterative  procedure  and  the  above  equations  indicate 
that  of  all  th“  TEM|,ju,|  modes,  those  with  higher  order  of  radial  variation  (m  1)  and  angular  varia- 

tion (n  _ 1)  have  higher  losses  and  phase  shift  than  does  the  TEMqqjj  mode.  Higher  order  mode 
losses  decrease  sharply  with  increasing  Fresnel  number  {N  ■ To  decrease  both  the  level 

and  the  number  of  undesired  higher-order  mode  frequency  responses,  the  Fresnel  number  of  the  milli- 
meter wave  FPS  is  made  as  small  as  possible,  with  the  mininmn  N set  by  the  mirror  diffraction 
losses.  This  effect  is  well  shown  by  Figure  15  of  Fox  and  Li,  reproduced  here  as  Figure  26.  From 
this  figure  it  is  clear  that  tlie  power  loss  is  20  times  greater  for  the  TEMqj^  mode  than  for  the 
TEM  OOq  mode.  From  this  figure  it  is  also  clear  liow  the  power  loss  of  the  mirror  (diffraction  loss) 
varies  as  a function  of  the  Fresnel  number  N.  At  N 1.0,  this  loss  (in  percent)  is  about  0,05,  and 
for  N 0.5,  the  loss  is  about  10.  Thus,  for  a factor  of  two  change  in  N , the  loss  increases  200 
times.  It  is  clear  that  the  optimum  range  for  N is  closely  restricted  to  a narrow  range  about  1.0  for 
there  to  be  a minimum  diffraction  loss  togetlier  with  good  higher  order  mode  suppression.  Kogelnik 
and  Li  (1966)  calculated  the  phase  shift  per  transit  for  both  TEM^j^^  and  TEMqq^i  modes  in  a 
confocal  resonator.  For  N 1.0,  the  phase  shift  for  the  TEMpp^  mode  is  about  95°,  while  for  the 
TEMgjq  mode  it  is  190°.  This  phase  shift  difference  of  95”  implies  that  the  two  modes  will 
resonate  at  different  mirror  separations.  This  separation  should  be  about  A/4,  and  in  fact  the 
higher  order  mode  for  the  cavity  built  for  this  experiment  had  this  higher  order  mode  displaced 
0.26A  from  the  TEM^^^^  mode,  in  good  agreement  with  the  above  prediction. 

Both  Boyd  and  Gordon  and  Fox  and  Li  give  as  the  Q-factor  for  a multimode  confocal  resonator: 

1^  2nb 
uA 


Q 


u\  \ A ^ 


and 


R - lb-;  I r; 


b,  b-Aj- 


A.  b - (b- 


A 


P^igtire  25.  Terminology  nppropriute  to  theoretical  treatment 

of  FPS. 


P'igtire  26.  Power  loss  per  Uansit  vs  N a'/b\ 
for  confocal  spherical  mirrors.  Dashed  curves 
for  circular  plane  are  shown  for  comparison. 
(After  Fox  and  Li  1961.) 


If 
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siiu’e  A[n  A/2?r/)l  ^ A 

b mirror-to-niinor  separation 
u >■  losses  ol  cavity,  per  transit 


A wavelentrth. 

This  relatiottslnp  would  seem  to  indicate  that  the  cavity  Q-factor  increases  directly  with  mirror 
separation.  This  is  true  up  to  the  point  wliere  mirror  dil'l'raclioii  effects  increase  the  losses  as 
fast  as  the  separatioti  iticreases,  even  allowiiift  for  the  iiicreasiiif;  tiiirror  size  (as  dictated  by  the 
fresiiel  eriterioti)  with  separation.  Cavity  Q-factors  greater  than  10*’  at  150  GHz  have  been  re- 
ported by  fretich  atid  Arnold  (1967).  Dependiiif;  upon  the  cotistriietioii.  coiipliiif;  technique,  atid 
the  precision  with  which  the  FPS  is  made  and  aligned.  Q's  of  50.000  and  greater  me  not  difficult 
to  achieve  at  millimeti'r  wavelengths. 

In  order  to  calculate,  at  least  to  first  order,  •'  requited  separation  of  tlie  mirrors  for  a given 
some  estimate  of  tlie  losses  mtist  lie  obtained,  ibis  loss  estimate  will  be  made  from  summing 
the  e.xpected  reflection,  diffraction  and  coupling  losses. 

Mirror  reflection  loss  is  given  by 


where  I'  }—LL  _ i 
S 

reflection  coefficient  of  the  metal  surface 

and  where 


and  all  other  symbols  have  their  standard  meaning.  For  tlie  aluminum  spherical  mirror 
3.72  ■ 10'  mho^m 


and  the  copper  flat  plate 

5.80  • 10^  mho/m 


and  since  ( 

9.7  10‘°c/s.  then 

1.857  X 10*^ 

^Cll 

2,32  X lO'*. 

Thus 

1-t 

2 

M - - = 1.08  X 10“*^,  for  aluminum 

S 

8.62  X 10  for  copper. 
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Figure  27.  Amplitude  distribution  of  TEM^^^^ 
mode  (at  any  beam  cross  section). 


The  splierical  mirror  was  made  of  aluminum 
since  a suitable  piece  of  round  stock  was  avail- 
able for  cutting,  grinding  and  polishing.  Time 
and  the  lack  of  a large  enough  piece  of  copper 
did  not  permit  this  mirror  to  be  made  of  copper. 
Tlie  initial  flat  was  turned  from  aliiniiiium,  but 
tlie  later  availability  of  a piece  of  copper  plate 
allowed  this  piece  to  be  made  from  copper.  This 
operation  was  not  time  consuming,  since  the  face 
of  the  flat  was  latlie  turned  within  half  a day’s 
time. 


The  coupling  loss  into  and  out  of  the  cavity 
is  the  ratio  of  the  total  waveguide  cross  sectional 
area  to  the  beam  waist  area  on  tlie  flat  mirror.  The  amplitude  distribution  for  the  fundamental  mode 
(TEMqqjj)  across  the  waist  (in  a radial  direction)  is  Gaussian  in  form.  The  parameter  of  the 
beam  waist  is  the  radial  distance  from  the  center  one  must  move  for  the  amplitude  to  fall  to  1/e  of 
its  maximum  value,  ’■’igure  27  is  a sketch  of  the  radial  beam  amplitude  distribution.  uJq,  often 
called  tlie  beam  radiu.s  or  spot  size,  is  given  by  the  relationship 


O 


bA 
2n  ' 


where  b is  Mie  radius  of  curvature  of  the  spherical  mirror  and  A is  the  operating  wavelength.  Wq 
designates  beam  radius  at  the  team  waist  where  d b/2. 

For  a mirror  of  72-in.  radius  of  curvature  and  jui  operating  wavelength  of  3.09  mm,  the  spot 
size  is  of  30-mm  radius.  In  assuming  a spot  size  of  30-mni  radius,  this  means  that  the  majority  of 
the  power  is  concentrated  in  a spot  with  this  radius  on  the  flat  mirror.  This  power  within  the  spot 
size,  as  a fraction  of  the  total  power  on  the  flat  mirror,  may  easily  be  found  as  follows. 

, For  F(x)  - cumulative  distribution  function  of  standardized  normal  random  variable 


then  2(1  - F(x))  area  in  wings  of  the  Gaussian  function  fur  ^ \ 2.  From  a table  of  the  normal 
distribution  function  1 - F(x)  0.0793,  and  hence  the  fractional  power  enclosed  by  the  spot  size  is 
0.84  of  the  total  power  on  the  flat  mirror. 

For  a type  RG-cg  waveguide,  the  area  enclosed  by  the  guide  v/alls  is  4.802  mm".  Thus,  the 
coupling  loss  is  approximately 


9.604  mm" 
2.82  X lO'"^  mm^ 


= 3.4  y 10'^. 


Mirror  diffraction  losses  are  not  so  easily  obtained.  However,  Fox  and  Li  indicate  that  a 
diffraction  loss  on  the  order  of  5 > 10““*  may  be  expected  foi'  a Fresnel  number  N i.o. 

Thus,  the  total  losses  (per  transit)  for  a Fabry-Perot  cavity  operating  in  the  dominant  mode 
with  the  mirrors  constructed  of  copper  and  aluminum  and  fed  with  two  type  99  waveguides  located 
symmetrically  about  the  center  of  the  flat  plate  are  the  sum  of  the  I'^R  mirtor  losses,  the  coupling 
losses,  and  the  diffraction  loss.  This  comes  to  = 5.84  x 10"^,  but  since  further  losses  must 
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arise  in  practice  due  to  mirror  alignment  errors  and  surface  irregularities  in  the  metal  optics  due  to 
machining  and  grinding  a more  reasonable  total  loss  is  taken  to  he  1,0  x 10“^. 

For  tliis  experiment,  a "practical”  FPS  size  would  be  a mirror-to-mirror  separation  of  72  in. 
(corresponding  to  a semi-confocal  mirror-to-flat  separation  of  36  in.).  For  these  data,  tlie  design  Q 
can  be  calculated.  Thus, 

n 2/rb 

Q = ^ 186,000. 


A Q of  186,000  implies  an  effective  path  length  of 


= 91.5  m. 


It  is  obvious  from  this  relationship  why  a high  Q is  desirable  for  increased  measurement  sensitivity. 

With  a mirror-to-mirror  separation  of  72  in.  (confocal),  tlie  optimum  mirror  radius  for  N 1,0  is 
found  from  the  relationsliip 


Thus  a 75.15  mm 
or  a diameter  of  5.92  in. 

This  completes  tlie  tlieoretical  considerations  appropriate  to  tlie  design  of  the  semi-confocal 
FPS  for  an  operating  wavelength  of  3.09  mm. 

The  spherical  mirror,  as  previously  indicated,  was  made  from  a 6-in.-diam  section  of  aluminum 
rolled  stn^'k.  Initial  forming  and  shaping  of  the  mirror  was  done  on  a lathe,  with  final  grinding  and 
polishing  done  in  an  optical  laboratory.  Since  fticilit  es  and  equipment  existed  for  grinding  and 
polishing  telescope  mirrors  for  a 72-in.  radius  of  curvature  (36-in.  focal  length),  this  was  set  as  the 
design  parameter. 

Deviation  from  a true  spherical  surface  is  measured  by  an  optical  fringe  pattern  test  known  as 
a Ronchi  test  (Ingalls  1970).  A true  spherical  surface  gives  a finge  pattern  consisting  of  alternating 
straight  light  and  dark  bands  across  the  visual  surface  of  the  mirror.  Bowed  bands  indicate  varying 
degrees  of  parabolism  and  hyperbolism  of  the  surface.  This  mirror  was  found  to  give  a nearly 
perfect  straight  fringe  pattern,  with  only  a slight  tendency  toward  parabolism. 

The  flat  pla.e  was  machined  on  the  lathe  from  a '/j-in. -thick  copper  sheet  to  a tolerance  of 
* 0.005  in.  The  final  diameter  of  the  flat  was  about  8.5  in.,  which  is  certainly  mucl!  larger  than  re- 
quired by  the  beam  waist  diameter.  RG-99  waveguide  coupling  feeds  were  brought  through  the  back 
of  the  flat  mirror  and  located  symmetrically  about  the  center  line  of  tlie  cavity,  witli  a guide-to- 
guide  separation  of  1.5  in. 

Figure  28  shows  the  completed  cavity  mounted  on  an  optical  bencli  and  installed  in  tlie  environ- 
mental chamber. 

Details  of  the  cavity  Q measurements  and  response  are  delayed  to  Chapter  VII,  pending  details 
of  the  electronic  phase  locking  and  frequency  sweeping  teciinique  described  in  the  following 
chapter. 
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Rigure  2S.  Two  views  of  cold  clumber  interior,  showing  detnils  of 
evaporator  mounting  and  FPS  placement. 


Electronic  frequency  locking  system 

The  task  of  making  frequency  shift  measurements  as  small  as  a few  kilohertz  '.t  07  nw.  , 


44  COMPLEX  REFRACTIVE  INDEX  OF  ICE  FOG  AT  A RADIO  WAVELENGTH  OF  3 MM 


To  FPS 


Figure  29.  Block  dhignim  of  clusaicul  millimeler-wiive  frequency  measuring  system 

using  variable-frequency  markers. 


fiom  a directional  coupler  and  routed  to  a harmonic  mixer.  The  other  input  to  the  harmonic  mixer 
is  typically  an  accurately  known  suhmultiple  of  a frequency  in  the  vant'e  ol  the  klystron.  This 
lower  frequency  signal  is  usually  established  hy  an  X-hand  phase-locked  source  refereiu;ed  to  a 
free  running  oscillator  near  100  MHz.  Since  the  100-MHz  oscillator  frotiuency  can  lie  measured 
with  a high  degree  of  precision,  the  frequency  of  the  X-hand  source  can  Ire  established  with  the 
same  degree  of  precision.  Beats  between  the  swept  klystrom  and  multiples  of  the  X-hand  source 
establish  frequency  “heiicli  marks”  appropriate  for  calibration  and  measurement  ol  the  high  Ire- 
quency  klystron  source.  Varying  the  frequency  of  the  100-MHz  oscillator  "moves”  the  marker  on 
tlie  oscilloscope,  tints  establishing  a calibration  of  the  frequency  scale.  This  system  is 
diagrammed  in  Figure  29.  The  principal  disadvantage  of  this  system  lies  with  the  stability  of  the 
essentially  free-running  high  frequency  klystron.  Modulator  amplitude  stability,  klystron  filament 
temperature,  and  klystron  beam  and  reflector  voltage  changes  all  contribute  to  drift  in  the  P-curve 
displayed  on  the  oscilloscope. 

Another  technique  for  phase-locking  klystrons  employs  a fixed  low  frequency  oscillator  whose 
signal  is  harmonically  multiplied  up  to  the  klystron  frequency.  This  harmonic  multiple  is  then 
mixed  with  the  klystron  signal  to  produce  an  IF,  which  is  in  turn  used  to  phase  lock  to  a reference 
signal  at  the  same  intermediate  frequem  y.  The  reference  oscillator  itself  may  be  linearly  swept 
F \f  about  /g.  The  function  of  hanuonically  multiplying  the  fixed  low  frequency  oscillator  signal 
up  to  the  klystron  operating  frequency  ;uid  simultaneously  mixing  this  signal  with  a portion  of  the 
klystron  output  can  be  accomplished  by  one  diode.  This  is  no  trivial  feat,  since  the  low  frequency 
fixed  oscillator  may  run  near  2 GHz,  and  thus  this  harmonic  generator/mixer  crystal  must  generate 
a significant  amount  of  harmonic  energy  at  the  50th  or  so  multiple  and  mix  this  with  the  klystron 
output  to  provide  an  IF  signal.  Recent  advances  in  high-performance  harmonic  mixer  diodes  have 
made  possible  the  phase-locking  of  klystrons  to  well  above  100  GHz.  The  principal  application  of 
this  phase-locking  system  to  data  has  been  in  the  area  of  millimeter  wave  radio  astronomy  for  the 
study  of  extraterrestrial  molecular  line  spectra  (Weinreb  1970).  A phase-locking  system  using  this 
special  type  of  mixer  diode  was  built  to  take  advantage  of  the  extreme  stability,  simplicity,  and 
precision  frequency  measurement  offered  by  this  technique.  The  following  discussion  refers  to  the 
block  diagram  of  this  instrumentation  shown  in  Figure  30. 
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\w  vhe  10-clB  directional  coupler  and 
The  output  from  the  Varian  VC-7 13H  klystron  ‘^^^f^tween  the  directional  coupler  and  the 
fed.  via  an  attenuator,  to  the  . j.^om  entering  the  harmonic  nnzer  and  causing  the 

FPS  cavity  prevents  reflections  from  tie  couoled  bv  a triplexer  assembly  to  the 

Z::,o\Lom  «nl0Ced.  a.  2 

IF  sysiom  and  a 2-OHz  treqnency  exln.  whlcl,  Inck"  .0  a l.isl.  m.lhple 

Is  no  problem  since  fairly  liir.b  P0W«  <-  ^ 100-MHz  osr-.ialor  can  oc  diroclly 

rrsrr:.:- 

''tlllrple  Of  Ibis  kno«n  2-»a  -^^cy  m ^ 

I by  an  absorption  ly2=  „ic  mTx«t  l“  IF.  Pbas.  locking  is  done  a.  the  IF  fre- 

the  fteqiwncy  penerated  within  freoiiencv  relationships  wllhout  dislortion  "> 

qiiency  since  mixers  Iranslale  amplitode  and  ’ • phase  lookinB  system  proved 

frequencies  for  ease  in  further  pocessrnR.  A 60  MHx 

convenient  for  use  with  this  system.  oscillator  input  fte- 

It  is  the  purpose  of  the  triplexe.  mixer  crystal  bias  curreni.  Since 

queue,  horn  the  OO-MHa  b„  S;  “ o^snainp  this  componen,  tor  a eWHx  IF,  its  croni. 

Zitre=Xre  me  Lear,  rn  F.me  31.  _ 

The  60.MHX  IF  »T' i*JXLe'L,ir.tom  ,h.;  IF  system,  .hen  amplified  prio.  ,o  inirc 
isolate  the  high  power  2-GHz  reference  -.cu 

duction  into  the  phase  locking  system.  ^ ^ 

The  second  reference  input  frequency  o ^ Af  MHz,  and  is  linearly  swept  by  an 

sweep  oscillator.  The  sweep  oscillator  frequen  y 
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Port  (2)  - (3)  0.40  dB@  60  MHz 
Port  (1)  - (3)  3.  3 dB  @ 2.  17  GHz 
measured  with  HP432A  power  meter. 


Components  mounted  in  2- chamber 
milled  aluminum  box  with  interior 
compartment  sizes  1.2"  wide  x 
1 . 55"  long  X 0.  68"  deep. 


AT  A RADIO  WAVELENGTH  OF  3 MM 

Insertion  Loss 

Port  (1)  - (3)  3.  35  dB  @ 2.  17  GHz 

Port  (2)  - (3)  0.  514  dB  @ 60  MHz 

(measured  with  HP85 1 B /85 1 B 

spectrum  analyzer) 

Cj  =2.7  p£  (low  ind.  cap.  ) 

(JFD  2R7C)  in  parallel 
with  3.  3 pf  ceramic 

==  43  pf/lOOOv  mica 

CFjiCF^  = .001  mfdfeed 
through  capacitors 

Lj  - 10  turns  (/24  nyclad  wire  on 
15000,  1/2  watt  carbon  re- 
sistor 

= 189  ph  peaking  coil 

: 1 mh  coil 

Rj  = meter  Shunt  resistor,  lOf) 
to  give  0-lOma  full  scale 
reading 

■ TNC  female  co.’inectors 
= Americon  #26805 


Figure  31.  Triplexer  'assembly, 


external  sawtooth  function  generator.  This  sawtooth  voltage  is  also  used  to  control  tlie  sweep  of- 
the  honzontal  time  base  of  the  oscilloscope.  The  sweep  oscillator  has  an  internal  mixer  and 
provisions  for  external  input  sources  to  be  added  to  the  mixer.  Tims  marker  pips  are  generated  by 
the  sweep  oscillator  frequency  at  frequencies  corresponding  to  the  external  oscillator  frequencies. 
These  pips  are  displayed  on  one  of  the  vertical  channels  of  the  oscilloscope  and,  by  adjusting 
these  pips  to  correspond  to  the  half-power  points  of  the  cavity  response  and  to  tlie  peak  of  the  Q- 
curve,  tlie  cavity  Q and  frequency  shift  can  be  accurately  measured  by  the  difference  in  frequency 
of  these  60-MHz  reference  oscillators.  In  practice,  frequency  measurement  to  within  a few  hundred 
cycles  precision  is  possible  at  an  operating  frequency  of  97  GHz. 

Since  the  output  of  the  phase-sensitive  detector  is  a dc  error  voltage  which  is  added  to  the 
klystron  reflector  voltage,  the  klystron  frequency  is  electrically  tuned  in  a direction  to  minimize 
this  dc  error  signal.  The  effect  of  this  is  to  phase  lock  the  klystron  in  such  a way  that  it  is  swept 
Af  in  frequency  by  the  60-MHz  sweep  oscillator. 

Whether  the  klystron  frequency  is  above  or  below  the  reference  frequency  by  the  IF  frequency 
can  be  resolved  by  observing  the  lock  points  as  the  klystron  reflector  voltage  is  increased  and/or 
decreased. 
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\ 


tigiiro  32.  View  of  electronic  burdwure  ufisociuted  with  klystron  phase  locking  system. 


A wave^Miide  feed  couples  the  cavity  to  an  ordinary  inillimeler  wave  crystal  video  detector  to 
sample  tlie  eavity  response.  This  siptial  is  coupled  to  one  of  tlie  vertical  ehannels  for  visual 
presentation  and  measurement  of  eavity  response.  Between  tlie  crystal  detector  and  tlie  flat  plate 
a precision  BF  power  attenmtor  is  located  for  calibration  and  test  purposes.  l''ifinre  32  is  a view 
of  the  microwave  plnnibiiift,  iiielnding  the  triplexer  and  liannonic  mixer  assembly. 

Environmental  chamber  and  refrigeration  system 

The  enviroiniiental  .system  consists  of  a Styrofoam  cliamber  of  interior  dimensions  6 ft  S in.  ■ 

4 ft  '21.8  in.  The  wall  thickness  of  this  cliamber  is  8 in.,  and  access  to  the  interior  is  pained 
tlironpli  the  front  side  wliicli  is  held  to  the  chamber  by  six  screw-type  extension  clamps.  This 
method  of  seenrinp  the  door  provides  a maximum  access  area,  rpiick  entry,  and  a mininuini  amount 
of  heat  leakage. 

The  refrigerant  evaporator  system  consists  of  two  parallel  evaporators,  each  with  its  own 
independent  evaporator  control  valve.  Each  evaporator  consists  of  approximately  72  ft  of  %-iii.- 
diam  copper  pipe  with  alnminnm  "disk”  fins  in.  in  diameter  and  a fin  radiator  density  of  S/in. 
One  evaporator  is  6 ft  .'5  in.  long  by  3 ft  wide  and  is  mounted  on  Styrofoam  support  blocks  which 
hold  it  near  the  chamber  roof.  The  other  evaporator  is  more  compact,  consisting  of  three  stacked 
coils  surrounded  by  a fan  shroud,  and  is  mounted  at  one  end  of  the  chamber.  The  fan,  a 600-cfni 
sriuirrel  cage  type,  draws  air  over  this  set  of  coils  and  discharges  this  air  upwards  and  over  the 
roof  evaporator  coils.  This  fan  povides  a considerable  amount  of  air  circulation  within  the  cham- 
ber. 

The  optical  bench,  containing  the  Fabry-Perot  cavity,  sits  on  the  floor  of  the  chamber  with  the 
flat  mirror  moniited  near  one  end  wall.  'Vavcguide  runs  tlirougli  the  chamber  walls  to  connect  the 


. vicv/  ui  cnamoer  ana  electronics. 


Figure  33.  Further  details  of  environmental  chamber  and  peripheral  equipment. 

millimeter  wave  electronic  components  to  the  Fabry-Perot.  Figure  28  is  a front  view  of  the  cliambe 
looking  in  on  the  system.  Fi^i.re  33  is  an  overall  view  of  the  chamber  and  a partial  view  of  the 
electronics.  The  access  port,  visible  in  the  door,  is  shown  in  this  photo.  Figure  34  is  a sketch  of 
the  interior  floor  plan  of  the  chamber.  Note  that  in  this  sketch  the  laser  beam  is  shot  through  the 

bv  UiflaseT  ^^^france  to  and  from  the  cLmber 

by  the  User  beam  is  aained  by  V/z-m.-diam  Lucite  tubes  sealed  at  each  end  with  clear  slide  cover 

K cSs  and  containing  several  ounces  of  calcium  sulfate,  a good  dessicant.  Thus  no  moisture  will 
condense  on  the  interior  glass  covers  as  tlie  temperature  is  loweied.-and  there  is  no  scattering  or 
dif.raction  of  the  beam  except  by  the  ice  fog  aerosols.  A photocell  detector  and  dc  amplifier  drive 
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ciB  loss  = -10  loK 


L 


Iq  = no-fog  quiescent  output  current  as  sliown  on  chart  paper 


reduced  output  ciurent  due  to  fog  scattering  and  absorption  effects. 


ni7e  fnfr'  ^ absorption.  From  this 

density  opt.;!  ftt  bl‘ ^ "--1 


94  " single-stage  3-hp  Copeland  condensing  unit  (Model 

. ) Inch  drives  the  paralleled  evaporators  within  the  chamber.  Each  evaporator  is  fed  throneh 

.ts  own  separate  expansion  valve,  as  each  expansion  valve  has  its  own  separate  tte^mls^a^^rbull 
f iermometer  attached  for  continuous  temperature  readout.  Thus  each  evaporator  may  be  -balanced" 
a justing  the  valves  to  give  the  same  or  nearly  the  same  refrigerant  evaporation  temperature. 

untiHhe  n"'!'"'  (empirically)  by  adjusting  these  valves  in  combination 

ntil  t'>e  suction  line  input  pressure  is  -5  in.  of  mercury  (vacuum).  Refrigerant  502  was  chosen  as 

wor  mg  uid  for  this  system  as  the  best  compromise  between  low-temperature  caoabilitv  sv^tP 
co„pa..t,ll.„  (.l,e  Pinal  oompiesso,  was  <leai„ned  to  R-im  „ne,.  and  elow 
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Figure  35.  Laser  extinction  due  to  0.3  normal 
density  optical  filter. 


Temperature  control  of  the  chamber  was  achieved  by  adjustment  of  the  pressure-sensitive 
switclies  on  tlie  inlet  and  outlet  lines  of  the  compressor.  While  not  designed  to  function  as  a 
tliermostat,  these  limit  switches  in  fact  work  quite  well  for  this  since  there  is  an  explicit  relation- 
sliip  between  evaporator  temperature  and  evaporator  pressure  for  each  type  of  refrigerant  used,  This 
system  does  not  allow  for  rapid  changes  in  setting  the  chamber  temperature,  but  after  the  appropriate 
pressure  settings  for  eacli  desired  chamber  temperature  were  experimentally  found,  this  technique 
would  hold  tlie  temperature  constant  within  ± 1°C  for  all  settings  used. 

Warm  fog  was  admitted  in  controlled  amounts  to  the  cUambeJ'  through  a slim  glass  pipe  in  one 
wall.  This  fog  was  generated  in  a teapot  on  an  electric  heater  (Fig.  36),  and  seemed  to  produce  a 
much  steadier,  more  uniform  fog  tlian  earlier  techniques  which  were  tried.  Fog  was  added  to  the 
chamber  only  after  electrical  and  temperature  stabilization  was  achieved  (meaning,  generally,  only 
after  several  liours  of  eperation),  and  in  measured  amounts.  A crude  measurement  was  achieved  by 
timing  the  fog  admission  time  with  a stopwatch.  Five  to  ten  seconds  would  be  sufficient  to  produce 
a very  dense  fog,  wliich  was  mixed  and  dissipated  by  operating  the  internal  chamber  fan.  Generally 
this  mixing/dissipation  was  done  until  the  laser  level  returned  to  within  -3  dB  of  its  initial  level. 
Tills  level  of  operation  was  chosen  since  this  fog  density  provided  for  an  easily  observed  cavity 
frequency  shift  and  satisfactory  impactor  measurements,  and  compared  with  a known  calibrated 
extinction  point  (the  0.3  ND  optical  filter)  of  the  laser  on  the  chart  record. 

Once  tlie  desired  fog  density  was  reached  (by  observing  the  chart-recorded  laser  level)  both 
fan  and  compressor  were  shut  off  to  minimize  vibrational  effects  on  the  cavity  operation.  All 
mechanical  operations  were  stopped  whenever  a photographic  record  and/or  frequency  measurement 
was  made.  At  this  level,  an  oscilloscope  picture  was  made  of  the  cavity  response  to  compare  with 
the  no-fog  cavity  response.  Interestingly  enough,  for  the  small  amount  of  steam  fog  admitted,  the 
chamber  temperature  changed  very  little  from  its  no-fog  condition.  This  change  was  at  most  2°C 
and  typically  1°C  or  less,  which  is  an  adequate  testimonial  to  the  low  temperature  capacity  of  this 
system. 
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Figure  3(1.  Teupol  mode  of  aieum  fog  production. 


VII.  MEASUREMENTS 


FPS  operation 

Prior  to  all  operations  and  ineasiueinents  tlie  klystron  and  klystron  power  supply  were  tlioroiighly 
warmed  np.  All  low  fieqneney  oscillators,  ineliidinp  the  97-MHz  relerence  oscill.’tors  and  tlie  60- 
MHz  heat  (marker)  oscillators  were  run  contiinionsly  dnrinp  tlie  data  tiikin;;  operations.  Prior  to  data 
takinft  and  nieasinenients  in  tlie  cold  chamher  it  is  necessary  to  rim  for  several  lioiirs  at  tire  desired 
temperature  to  assure  tliat  adequate  tliermalizalion  in  the  nieehaiiieal  parts  of  tire  cavity  has  taken 
place. 

Longitudinal  mirror  iiioveinent  was  remotely  controlled  by  a motor  suitably  peared  down  so  as 
to  move  the  splierical  mirror  very  slowly.  Thus,  as  the  air  was  cooled  down  and  dried  out  tire  index 
of  refraction  decreased,  and  lienee  the  resonant  frequency  of  the  cavity  increased.  For  air  saturated 
with  water  vapor  at  room  temperature  there  is  a chanpe  of  over  100  N-uuits  as  the  air  is  cooled  to 
-35°C  and  becomes  saturated  with  lespect  to  ice  at  this  temperature.  Rather  than  disturb  the 
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reference  2-GHz  oscillator  setting  to  track  the  cavity  frequency,  th.i  spherical  mirror  was  moved  by 
the  remote  controlled  motor  to  keep  the  Q-ciirve  on  the  oscilloscope. 

The  reference  oscillator  frequency  controllinK  the  2-GHz  locked  oscillator  was  kept  nem’ 

97,9213  MHz.  and  the  locked  oscillator  output  was  1762.5834  MHz  (corresponding  to  the  18th  multi- 
ple of  the  97-MHz  input).  A spectrum  analyzer  check  indicated  1760  MHz.  For  a multiplicative 
factor  of  .55.  the  klystron  frequency  is  therefore  96.942.087  MHz  ^ IF.  An  absorption  wavemeter 

check  indicated  97  GHz. 

The  purpose  of  i.icludi.if;  the  RF  waveguide-type  calibrated  attenuator  in  the  crystal  detector 
waveguide  from  the  FP  cavity  was  to  check  the  square-law  characteristic  of  the  detector  diode. 

For  a square-law  detector,  tlie  output  voltage  of  the  detector,  y.  goes  as 

y k ■ 

where  k < onsta.it,  charactc.ized  by  the  amplitude  of  <»ie  sine  wave  signal  divided  by  the  RMS 
noise  amplitude 

V voltage  impressed  upon  the  crystal  by  the  input  waveform  prior  to  detection. 

Since  input  power  is  proportional  to  v=^.  the  linear  measure  of  y is  a measurement  of  the  input 
signal  power  to  the  crystal  detector. 

The  square-law  character  of  the  video  crystal  was  checked  by  changing  the  attenuator  sefiiig 
3 dB  and  observing  the  response  on  the  oscilloscope.  The  detector  law  is  in  fact  square-law  since 
the  peak  of  the  curve  dropped  by  one-half. 

Figure  37  shows  two  measurements  of  the  cavity  Q made  in  the  open 
During  this  measurement  the  relative  humidity  was  35%  and  the  temperature  74  F.  The  Q is  found 

from  the  relation 


Q ^ 

A( 

where  f klystron  center  frequency.  - 97  GHz 
\f  : half-povier  frequency  separation. 

The  beat  markers,  visible  at  the  bottom  of  the  photographs,  were  set  to 

left-half  power  points  of  the  Q-curve  by  adjusting  the  frequencies  of  the  two  60-MHz  marker  oscilla 
tors.  The  variation  in  the  Q between  these  two  curves  is  attributed  to  the  difficulty  of  setting  these 
marker  oscillators  exactly  at  the  half-power  points. 

Figure  38  shows  the  cavity  response  curves  in  the  cold  chamber  at  a temperature  of  -35°C 
The  video  amplifier  gain  and  the  horizontal  sweep  width  (the  Af  of  the  60-MHz  sweep  genera  or) 
were  reduced  to  bring  the  curve  down  to  full  scale  on  the  oscilloscope  and  spread  the  curve  out 
increased  Q measurement  accuracy,  respectively.  Still,  it  is  difficult  to  J 

cavity  Q because  of  this  difficulty  in  determining  the  marker  separation.  However,  the  cavity  Q 
a cold,  dry  environment  was  taken  to  be  140,000,  in  good  agreement  with  the 
186  000  Undoubtedly  this  measured  Q would  have  been  higher  but  the  mirror  vertical  travel  slide 
, “upper  prop  beL  the  amplUude  of  ,he  <!  curve  peaked.  No  etfor.  was 
mirror  so  as  to  Increase  the  availabie  travel  because  of  the  lime  limitations  placed  upon  some  ot 

the  borrowed  equipment. 
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0 ' 107 ,000  : 9081  MHz 

inure  37.  hibry-Peroi  response 
curves  in  luborutory  environment  (rel- 
ative humidity  35%.  temperature  74  F). 
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Covity  responi*  m dry  , cold  chomber 
(-35-C  ) , At  = 0.6619  MHz  0 = 146,500 


(-35*0  ) , At  '0.7205  Mhz  Q'135,000 

f'ignre  3S.  Cavity  response  curves 
in  dry,  cold  chamber  with  no  log. 


ij3 


or  ,h.  dielecric 


wliere  <?„  cavity  Q in  a cold,  dry,  fotrless  clianiber 
cavity  Q after  introduction  of  fop 
A operating  wavelength,  in  cm 
with  a loosely  coupled  square-law  detector. 


tlierefore 


i 

i 


X 
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, — plot  Mirror 


I- 

Figure  39.  Method  for  pliotogruphing  ice  particle 
accnmiilutions  on  mirrors. 


Changes  .n  amplitude  are  far  more  accurately  and  easily  measured  than  changes  in  Q.  Therefore, 
loss  coefficient  measurements  are  best  found  in  this  manner  where  tlie  amplitudes  are  directly 
scaled  from  the  oscilloscope  photographs  taken. 

Since  it  is  possible  to  read  the  amplitude  of  the  cavity  Q-response  curve  on  the  oscilloscope 
photographs  to  about  0.3%,  the  minimum  detectable  loss  will  be  about 


.V,nin  1.45.10-'‘ 


ler-l 


2 • 10  ^ cm  K 


FPS  mirror  effects 

The  effects  of  ice  accumulation  >n  the  cold  mirrors  was  investigated  in  an  empirical  manner. 
Several  schemes  were  attempted  and  discarded  in  an  effort  to  measure  mirror  ice  accumulation. 

Direct  photograptiy  of  the  mirror  surfaces  invariably  was  inconclusive  for  any  viewing  condition 
witti  white  light.  However,  a technique  whereby  red-sensitive  film  was  used  to  photograph  laser 
scattering  due  to  small  particles  on  the  mirror  surface  worked  surprisingly  well  for  accumulation 
and  documentation  studies. 

With  this  technique,  a laser  beam  is  directed  tangentially  across  the  flat  mirror  and  so  adjusted 
that  there  is  only  the  barest  grazing  contact  with  the  mirror.  When  so  adjusted  the  bea.m  itself 
leaves  no  visible  trace  across  the  mirror.  However,  dust  particles  on  the  mirror  surface  too  small  to 
see  directly  witli  ttie  unaided  eye  cause  a distinct  and  higlily  visible  amount  of  light  scattering. 
Photographic  exposures  of  this  light  scattering  effect  are  made  tlirough  a 2-in.-diam  Lucite  tube 
penetrating  through  the  Styrofoam  wall  and  sealed  at  each  end  witli  clear  slide  glass.  Frost  forma- 
tion on  the  interior  glass  cover  is  prevented  by  the  addition  of  a small  amount  of  drying  compound, 
Ca^SO^,  The  viewing  tube  is  so  positioned  that  about  a 3‘/a-in.  length  of  the  grazing  beam  on  the 
mirror  is  visible,  and  the  scattered  light  to  the  viewer  is  at  an  angle  of  about  40'  from  the  forward 
direction.  Figure  39  is  a sketch  of  this  arrangement. 

With  this  viewing  technique  it  is  easy  to  see  that  dust  in  the  air  is  amazingly  abundant.  A 
clean  mirror  surface  does  not  stay  that  way  for  long.  Any  passing  air  current  carries  with  it  dust 
which  builds  up  a rapid  accumulation,  and  in  a short  period  of  time  these  dust  particles  sticking  to 
the  mirror  “puint  in"  a solid  red  line. 

A relative  size  calibration  can  be  had  by  dusting  the  mirror  surface  with  a 3-p  polystyrene  dust 
and  photographing  the  scattered  light.  Figure  40a  shows  a clean  mirror  surface.  Several  dust  grains 
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C/ean,  dust-free  mirror. 


N 


3(1  polystyrene  ■‘sphere"  dust  on  mirror. 


Chomber  temp  -4b  C,  no  lop,  ahoiil  III  dust  p.irtlcles  visible. 


Chumher  (omp  -4b.5''C.  Picture  mode  4 m/niilcs  niter  strum  loi; 
introduced  In  c/iamPcr-iipoui  4:i  purllclrs  visible. 


Figure  40.  Phologriiphs  made  from  6328^ 
light  scattered  hy  particles  on  mirror  sur- 
face. 


are  visible  in  the  rittlit  center.  Figure  40b  is  a 
photograph  of  this  same  surface  very  lightly  dusted 
with  the  3-p  polystyrene  dust.  This  dusting  was 
accomplished  by  putting  a small  amount  of  the 
polystyrene  dust  on  a piece  of  clean  tissue  and 
shaking  the  tissue  about  6 in.  in  front  of  the  mirror. 
While  it  was  strikingly  visible  in  the  photograph, 
this  amount  of  dust  was  almost  undetectable  to  the 
eye  over  other  parts  of  the  mirror.  Some  clumping 
is  evident,  as  seen  by  the  larger  areas  of  scattered 
light,  but  it  is  clear  that  this  technique  is  quite 
sensitive  to  detecting  light  scattering  from  small 
particles.  The  mirror  was  then  cleaned  as  well  as 
possible,  and  the  chamber  temperature  lowered  to 
-47°C.  Figure  40c  is  a picture  of  the  mirror  taken 
at  this  temperature  with  no  fog  having  been  admitted 
to  the  chamber.  Dust,  from  the  interior  of  the  cham- 
ber, stirred  up  by  the  evaporator  fan,  is  clearly  seen 
Figure  40d  was  made  about  four  minutes  after  a 
dense  ice  fog  was  formed  in  the  chamber.  The 
accumulation  seen  here  formed  very  quickly  after 
the  admission  of  the  warm  steam  fog,  and  no  further 
accumulation  was  noted  as  time  went  on,  even 
though  the  fog  remained  optically  quite  dense  for 
better  than  eight  minutes.  It  is  clear  that  the  fog, 
while  still  in  the  liquid  state,  will  adhere  to  the 
mirror  and  freeze,  but  once  frozen  the  ice  fog  crys- 
tals bounce  off  the  cold  surface.  This  effect  was 
noted  as  the  fog  particles  became  visible  in  the 
beam,  recoiled  from  the  smface,  and  disappeared 
out  of  the  beam.  The  freezing  time  for  these  small 
droplets  at  this  temperature  is  apparently  quite 
rapid.  If  a measure  of  the  freezing  time  can  be  had 
by  measuring  the  accumulation  time  during  which 
the  particles  will  adhere  to  the  mirror,  then  the 
steam  fog  is  largely  crystallized  within  5 seconds 
of  admission. 


The  cavity  response  was  undisturbed  by  a thick  polystyrene  dust  coat  on  the  mirror.  Polysty- 
rene dust  was  applied  to  the  mirror  until  a thin  white  coat  covered  the  surface.  Figure  41  i , the 
FPS  response  curve  made  before  and  after  this  dust  was  applied.  There  is  no  visible  change  in  the 
Q-curve.  The  same  sort  of  effect  was  later  noted  with  a thin  film  of  frost  on  the  mirrors.  This 
frost,  brought  on  by  applying  warm,  moist  air  directly  to  the  cold  mirrors,  had  no  discernible  effect 
on  the  cavity  response.  However,  any  degree  of  liquification  of  the  frost  caused  the  cavity  ampli- 
tude response  to  die  out  dramatically.  This  same  effect  was  previously  demonstrated  during  some 
early  experimentation  with  this  cavity  (Fig.  42). 


Ice  fog  particle  measurements 

Physical  measurements  of  the  ice  fog  were  made  separately  from  the  electromagnetic  measure- 
ments. These  measurements  of  the  particle  sizes  and  densities  were  made  at  the  same  temperatures 
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Figure  41.  Liberally  dusting  mirror  with 
3~n  polystyrene  dust  produces  almost  no 
effect. 
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Figure  42.  Q-perturbation  due  to  a thin  film  of  water  on 
mirror  surface  (expanded  s'  eep). 


as  electromagnetic  measuiements  of  the  fog  were  made,  i.e.  -30°,  -35°,  -39°,  -43°  and  -47°C. 

These  fog  measurements  were  made  at  or  near  the  -3-dB  extinction  level  which  was  taken  to  be  the 
standiird  operating  point  for  each  fog  and  electromagnetic  measurement.  Thus  all  measurements  over 
the  temperature  interval  above  are  referenced  by  the  common  laser  extinction  level.  This  laser 
level,  as  mentioned  earlier,  was  chosen  because  of  the  calibration  established  by  the  known  trans- 
mission of  an  optical  filter,  as  well  as  its  effect  in  producing  a visible  change  of  the  cavity  response, 

A UNICO  four-stage  cascade  impact  or  was  used  to  collect  fog  samples  for  study.  This  impactor 
uses  two  glass  slides,  coated  with  magnesium  oxide,  as  impact  areas  for  the  ice  fog  particles.  The 
n ignesium  oxide  slide  coating  is  formed  by  holding  the  clean  slide  over  a burning  magnesium  ribbon 
(Rothman  and  Ledbetter  1968).  This  impactor  is  so  con.structed  that  it  contains  a sequence  of  four 
internal  knife  edges  over  which  the  sample  volume  of  air  containing  the  aerosol  is  drawn.  These 
knife  edges  are  seriesed,  and  each  edge  is  closer  to  the  slide  than  the  preceding  knife  edge.  In 
theory,  the  most  massive  aerosol  particles  have  too  gieat  a mouientmii  to  make  the  turn  with  the 
airstream  around  the  first  sh.arp  edge.  Their  momentum  sends  tliem  crashiiig  into  the  magnesium 
oxide  layer  to  leave  a crater,  the  size  of  which  is  a measure  of  the  size  of  the  ice  fog  particle. 

The  aerosol  particles  are  “sized”  stage  by  stage  in  this  manner,  v/ith  the  fourth  stage  collecting 
the  smallest  particles.  A rotometer  and  stop  watch  are  used  to  determine  the  volume  of  air  drawn 
through  the  impactor,  and  from  a count  of  the  number  and  sizes  of  the  impact  craters  a number-density 
distribution  can  be  found. 

While  in  theory  this  measurement  technique  sounds  simple,  in  practice  it  is  subject  to  many 
errors  which  make  these  counts  unreliable.  Since  this  device  is  hand-controlled  and  -operated,  it  is 
clear  that  each  individual’s  technique  and  skill  directly  affects  the  final  count.  The  average  of  a 
limited  number  of  samples  was  considered  as  the  “best  estimate”  of  the  actual  fog  number-density. 
The  second  major  limitation  of  this  method  centered  around  the  very  poor  contrast  of  these  fog 
particle  cratei  j>  in  the  magnesium  oxide  film.  Figure  43  shows  a photomicroscopic  picture  made  of 
a given  impact  area.  Fatigue  after  several  hours  of  close  study  was  undraibtedly  responsible  for 
some  error.  The  third  area  for  measurement  error  involves  the  collection  efficiency  of  the  impactor. 

No  published  data  were  found  which  attempted  to  answer  this  question.  The  collection  efficiency, 
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Figure  43.  Plioloniicrogruplis  of  10  //  (upper)  and  7 p (lower)  craters  left  by 
ice  fog  particles  on  a magnesinw  oxide  slide. 


especially  lor  the  smallest  particles,  must  certainly  be  poor  since  on  several  slides  there  were 
definitely  traces  of  craters  smaller  than  about  4 p diameter  wliicli  were  of  too  poor  a contrast  to 
count.  Attempts  to  improve  their  readability  by  depositing  thicker  and  thinner  magnesium  oxide 
coats  did  not  help,  nor  did  efforts  at  staining  the  slide.  The  meteorology  group  at  the  U.S,  Army 
Cold  Regions  Research  and  Engineering  Laboratory  uses  a single-stage  impactor  which  samples 
10  cm'  of  volume  and  causes  the  aerosol  particles  to  impact  upon  a small  area  on  a glass  slide 
containing  a drop  of  silicon  oil.  A photomicroscopic  picture  thus  "sees”  the  entire  particle 
collection,  and  can  later  be  enlarged  for  counting.  However,  this  technique  requires  tha*  the 
ambient  temperauire  he  below  freezing  before  and  during  the  time  the  particles  are  being  photo- 
graphed to  prevent  particle  evaporation.  Tliis  was  clearly  not  possible  for  this  experiment,  thus 
ruling  out  this  technique. 
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Experiment 

no. 


Table  VI.  Measured  partible  diameter  (number/ cni^). 

N 

<?H  7-1Z\l  I2l!)ft  19-2611  >26  fi  total 


T - -47  + 0.5°C 


Laser  Mtinction 
(6328  A)(dB/m) 


1 

_ 

121 

13 

3.73 

1.12 

189 

1.33 

2 

— 

155 

4.15 

0.218 

- 

159 

1.16 

3 

29.50 

164 

8.62 

0.144 

0.088 

202 

1.20 

Avg 

10 

164 

8.6 

1.4 

.406 

183 

1.23 

N.J, 

.5.5 

89.5 

4.7 

,3 

183 

— 

T - 43''' 

to 

167* 

12.2 

179 

1.42 

11 

lf)2» 

6.35 

168 

1.50 

Avg 

165» 

9 

174 

1.46 

”o  N »|i 

94. 8» 

!..2 

T -39  + .5°C 

8 

151.2* 

12.7 

163.9 

1.46 

o 

174.5* 

1.52 

176.0 

1.42 

Avg 

163* 

7.1 

170 

1.44 

‘'o  N T 

96* 

4 

— 

T -35  + .5°c 

4 

85.5 

17.2 

0.58 

1.26 

104.6 

1.37 

5 

_ 

164 

76.4 

3.25 

0.8 

244.5 

1.33 

() 

117 

48 

.58 

1.2 

166.8 

1.14 

7 

10.7 

74 

10.2 

4.6 

.59 

100.1 

1.04 

Avg 

2.68 

no 

38 

2.25 

.96 

154.2 

1.24 

N-r 

1.68 

71.3 

24.7 

1.46 

.624 

T = -30°  C 

18 

96.2 

15 

4 

133.2 

1 45 

6-12g 

12-20/i 

’20g 

216 

37 

32 

32 

317 

1..50 

•rvl2/i 


Tlie  number  densities  repotted  in  Table  VI  are  therefore  lower  limits,  particularly  regaroin^t  the 
smaller  sizes.  Table  VII  is  a distillation  of  the  raw  numbers  appearing  in  Table  VI.  Particle 
diameters  reported  here  were  made  by  comparison  of  the  crater  sizes  witli  a microscopic  scale  (in 
tire  microscope  eyepiece)  calibrated  from  a stage  micrometer. 

Several  conclusions  can  be  drawn  from  this  work.  From  these  tables  it  is  clear  that  tlie  large- 
size  ice  fog  particles  are  associated  with  the  warmer  fogs;  the  majority  (95%  or  better)  of  the 
particles  associated  witli  fogs  39"C  or  lower  have  diameters  less  than  about  12  /i.  For  the  same 
laser  extinct ioti.  the  total  number  densities  ol  ‘liese  fogs  at  all  temperatures  are  comparable.  The 
lowest  total  number  density  (100/cm''^)  and  the  highest  number  density  (317/cni^)  are  for  the  warmer 
reported  fogs,  and  tlie  lov/er  temperature  fogs  had  intermediate  number  densities.  j 


5 

t 
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Table  VII.  Gross  statistics. 

Temp 

N/cm^/dB  'm 

Particle 

size  disiribinion  (%) 

(°C) 

6328^ 

5-12ii 

12-19fi 

l!)p 

-30  ^ .5 

149 

217 

34 

38 

8 

-33  ^ .3 

1!?4 

1.34 

73 

23 

2 

-39  ♦ .3 

118 

170 

96 

4 

- 

^ .3 

-43 

- 1.0 

119 

174 

94.8 

5.2 

- 

-47  + 1.0 

149 

183 

93 

4.7 

.3 

Perhaps  the  most  significant  finding  of  all  is  the  reproducibility  of  these  fogs.  For  fogs  -39“C 
and  colder,  the  measured  number  densities  repeated  with  excellent  agreement  from  experiment  to 
experiment.  What  measurement  difficulties  there  were  arose  witli  the  warmer  fogs.  In  the  four 
experimental  measurements  made  on  the  -35°^  fog  the  number  densities  varied  by  as  much  as  a 
factor  of  2.5.  The  characters  of  the  -30°  and  -35°C  fogs  were  much  different  tlan  those  of  the 
colder  fogs.  These  fogs  would  dissipate  more  rapidly  and  generally  rerjuired  tlie  admission  of  more 
steam  fog  for  formation  and  mixing.  While  measured  amounts  of  steam  fog  at  the  colder  temperatures 
would  generally  give  a reproducible  effect,  such  was  not  always  the  case  here.  These  effects  were 
even  more  pronounced  at  -30°C,  but  due  to  the  critical  nature,  of  ice  fog  formation  in  this  temper- 
ature interval,  no  great  efforts  were  made  to  describe  fogs  at  -30°C.  However,  electromagnetic 
measurements  were  made  at  this  temperature  at  the  same  laser  extinction  level  as  were  the  lower 
temperature  measurements. 

These  measured  number  densities  were  used  as  a basis  for  calculating  the  laser  extinction  to 
compare  with  the  measured  extinction.  The  calculations  were  made  by  assuming  the  ice  fog 
particles  were  spherical  and  applying  Mie's  theory  to  calculating  the  absorption  and  scattering  cross 
sections  of  each  particle  size.  This  procedure,  as  well  as  the  optical  depth  calculation,  was  out- 
lined in  Chapter  IV.  Figure  44  shows  five  number-density  distributions  used  in  these  calculations. 
Curves  A and  B are  more  nearly  representative  of  the  actual  impactor  measurements  made  of  these 
fogs,  while  C,  D and  E weight  the  small  particle  sizes  quite  heavily.  Table  Vlll  lists  the  calculated 
and  measured  optical  properties  of  tliese  fogs.  Table  IX  tabulates  the  extinction  cross  sections  for 
ice  spheres  from  1 to  35  //  diameter  for  6328  A laser  light. 

These  calculations  indicate  several  interesting  things.  First,  the  small  particles  are  more 
effective  scatterers  than  the  larger  particles  (not  surprisingly).  However,  a rjlatively  small  number 
of  larger  particles  can  contribute  significantly  to  the  total  extinction  (distribution  B,  Table  VII). 
Secoiia':',  there  is  a factor  of  between  3 and  5 difference  in  the  measured  and  calculated  attenuation. 
Recently  other  investigators  have  reported  similar  results  (Munis  and  Delaney  1972,  J.  Albertine  of 
the  Jolins  Hopkins  App’.jq  Physics  LalAiratoty,  rxivate  cotutiiunlcallctt).  This  diecc^paucy  is  fc%!ir4; 
studied,  but  as  yet  there  is  no  proven  explanation. 

The  principal  point  here  is  the  obvious  difference  between  the  measured  and  calculated  optical 
attenuation.  Three  things  come  to  mind  as  possible  explanations:  1)  the  refractive  index  for  ice  at 
this  w.ivelength  is  inaccurate,  2)  the  measured  number-densities  :ue  wrong,  and  3)  Mie  theory  cannot 
be  applied,  since  these  particles  may  be  mostly  geometric  figures  with  flat,  prismatic  faces.  While 
It  is  certainly  uue  ..hat  improved  measurements  of  ilie  complex  inJe*  oT  iefuu.uoli  ol  .r  e al  6328  A 
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Table  VID.  Cplculated  and  measured  6328  A optical  properties  ot  ice  fog. 


Distribution 

Total 
N cm'^ 

Percentage  ol  particles 
ot  diameters  (%) 

6(1  6-12/1  12/1 

r. 

extinction, 

Calculated 
loss,  dB/m 

Measured  loss, 
dB 'm, 

where  appropriate 

A 

183 

0 

94 



.0271 

.239 

1.23 

B 

170 

29.4 

41.2 

29.4 

.0431 

..380 

1.24 

C 

llSFi 

37.8 

02.2 

- 

.0920 

.816 

n 

•180 

04.2 

35.8 

- 

.0290 

.201 

_ 

E 

4308 

87 

13 

- 

.1071 

.931 

- 

Table  IX.  Optical  extinction  cross  section  for  ice  spheres  at  6328  A 


= 1-302 


Total  extinction  o 

•c  p y rp 

cross  section, 


Diam 

^EXT 

o 

/7f*- 

1 

2.56  V 

3.26 

2 

8.32  X 10“*^ 

2.65 

3 

1.50  r 10~^^ 

2.12 

4 

3.02  X 10~" 

2.41 

5 

3.89  X 10~'* 

1.98 

6 

6.79  >'  10' 

2.40 

7 

7.72  X 10~^* 

2.01 

8 

1.21  X 10~^° 

2.40 

9 

1.27  X 10“'” 

2.00 

10 

1.81  X ur'” 

2.30 

11 

1.86  • 10“'” 

1.95 

12 

2.55  X 10“'” 

2.25 

13 

2.04  X 10“'” 

1.99 

14 

3.43  X 10“'” 

2.23 

1.5 

3.. 57  X 10  '” 

2.02 

16 

V 

C 

o 

2.20 

17 

4.60  X 10  '” 

2.03 

18 

5.43  X 10“'” 

2.13 

j 1.9  X 10“®). 


7’otal  extinction  a 

**  P X T 

cross  section, 

Di  am  QeXT  nr~ 


19 

5.80  X 10“'” 

2.05 

20 

6.64  X 10“'” 

2.12 

21 

7.21  X 10“'” 

2.08 

22 

7,98  X 10“'” 

2.10 

23 

8.76  X 10“'” 

2.11 

24 

9.40  X 10“'” 

2.03 

25 

1.04  X 10“” 

2.12 

26 

1.08  V 10“” 

2.04 

27 

1.21  X 10“” 

2.12 

28 

1.25  X 10“” 

2.03 

29 

1.40  X 10“” 

2.12 

30 

1.43  X 10“” 

2.03 

31 

1.60  X 10“” 

2.12 

32 

1.61  X 10“” 

2.01 

33 

1.81  X 10“” 

2.12 

34 

1.83  X 10“” 

2.02 

35 

2.03  X 10“^ 

2.12 

can  be  made,  it  is  not  likely  that  these  new  measurements  would  be  radically  different  from  what 
was  repo- ted  by  Irvine  and  Pollack  (1968)  for  ice  at  1.0  n- 

In  regard  to  the  problem  of  making  accurate  fog  number-density-measurements,  it  is  worthwhile 
to  quote  some  of  the  remarks  of  Rinehart  (1969),  who  made  a study  of  the  methods  of  measurement 
for  fog  drop  size  distributions.  She  studied  capture  ar'^/or  replication  media  which  included  poly- 
vinyl alcohol,  gelatin,  oil  and  magnesium  oxide,  and  ■ acluded  that  fog  drop  size  distributions 
were  method-dependent  and  none  of  the  above  media  was  sensitive  to  or  retained  for  lengthy  periods 
the  impressions  of  small  drops.  Furthermore,  she  noted  that  ‘‘lack  of  correlation  between  number  of 
captured  droplets  and  visibility  was  present  in  even  the  most  careful  determinations  employing 
impaction  methods.” 
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Figure  44.  Assumed  number-density  distribution  used  in  laser  ex 
tinction  computations. 


One  other  sampling  effect  should  be  mentioned.  The  craters,  seen  in  Figure  43,  may  not 
accurately  reflect  the  true  impacting  particle  size.  A solid  spherical  particle  striking  a thin 
magnesium  oxide  film  would  only  crush  the  film  down  to  the  point  of  contact  of  the  particle  with  the 
underlying  glass.  This,  in  effect,  would  leave  a “crater  smaller  than  the  true  size  of  the  particle. 
This  effect,  while  operative  mainly  with  the  larger  particles,  also  indicates  that  the  measured  ice 
volume  is  really  a lower  limit  to  the  true  ice  particle  volume  per  unit  volume  in  the  cold  chamber. 

It  is  clear  from  these  fog  impactor  measurements  that  considerable  uncertainty  exists  about  the 
true  ice  fog  number^density  distribution.  This  is  particularly  true  in  regard  to  the  smaller  particles 
(those  less  than  about  7 g in  diameter)  since  the  impactor  registered  these  particles  poorly  if  at  all. 
There  does  not  exist  at  this  time  a method  or  technique  suitable  for  making  an  accurate  measure- 
ment of  the  fog  number-density  distribution. 

Although  desirable,  it  is  therefore  not  possible  to  refer  the  electromagnetic  measurements  to  a 
fog  particle  number  density.  However,  a suitable  reference  fog  level  is  the  measured  extinction  of 
a 6328  A helium-neon  laser.  It  was  shown  that  comparable  levels  of  laser  extinction  give  comparable 
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Table  X. 


Electromagnetic  measurements,  normalized  to  1 dB/m  at 


6328^  extinction. 


Laser 


T 

(°C) 

Ar 

{kHz) 

N 

(‘xtinction 
(dB  m) 

-16 

228.6 

2.36 

1.46 

-16 

288.3 

2.16 

1.35 

45.6 

229.6 

2.37 

1.46 

-46 

236.9 

2.44 

1.42 

-46 

2.88.1 

2.11 

1.  16 

46.5 

262.  1 

2.71 

1.90 

-A7 

236.0 

2.48 

1.28 

-46.5 

260.0 

2.68 

1.80 

-16.5 

224.5 

2., 82 

1.17 

-42.5 

250.5 

2.59 

1.57 

-4.'i 

260.9 

2.70 

1.55 

-42.5 

340.9 

3.52 

1.55 

-48 

298.1 

3.08 

1.  19 

-48 

203.4 

2.10 

1.13 

42.5 

279.7 

2.89 

1.40 

-42.5 

324.4 

3.35 

1.61 

-18 

322.5 

3.38 

1.57 

-14 

298.0 

8.08 

1.44 

-12 

286.8 

2.90 

1.51 

-12 

587.0 

6.05 

2.5;} 

-42.5 

881.2 

3.42 

1.48 

-48 

261.8 

2.70 

1.72 

-48.5 

537.8 

5.55 

3.00 

-48 

284.9 

2.42 

1.25 

-89 

769.7 

7.94 

2.76 

-89 

772.0 

7.56 

2.92 

-89 

825.4 

8.51 

8.03 

38 

398.3 

4.05 

2.13 

38 

424.4 

4.37 

2.02 

-38 

333.2 

3.48 

1.61 

-37 

328.2 

8.39 

1.35 

-88 

270,5 

2.79 

1.28 

-88 

302.9 

3.13 

1.24 

87 

275.9 

2.84 

1.28 

-88 

333.2 

8.44 

1.59 

-39.5 

331.8 

3.42 

1.31 

-88 

241.9 

2.. 50 

1.22 

-38 

429.0 

4.43 

1.64 

-38 

267.9 

2.76 

1.77 

-38.5 

361.9 

8.73 

1.40 

■88.5 

346.5 

3.58 

1.51 

' dB  'm 

cni~  * 

X/dB/tn 

1,62 

5.50  • 10“'^ 

3.76  X 10 

1 .82 

10.40  V 10 

7.70  / 10  ^ 

1.63 

6.25  > 10 

4.28  ■;  10 

1.72 

9.55  '•  10  ^ 

6.71  .y  10"'^ 

1.65 

8.70  y 10~^ 

5.95  > 10  ^ 

1.48 

4.80  y 10  ' 

2.58  X 10  ^ 

1.90 

5.50  10“^ 

4.80  X 10  ^ 

1.49 

6.40  10  ^ 

8.56  ■.  10  ^ 

1.98 

7.25  ' 10  ^ 

6.20  X 10 

1.6.5 

9.57  > 10"^ 

6.10  X 10“^ 

1.74 

10.40  y 10“^ 

6.70  < 10  ^ 

2.27 

10.30  10“^ 

6.65  ■:  10 

2.59 

9.00  X 10“^ 

7.55  X 10 

1.86 

9.15  V 10  ^ 

8.10  X 10  ^ 

2.06 

1.94  X 10“^ 

8.52  X 10 

2.08 

7.25  X 10'  ^ 

4.50  X ur^ 

2. 12 

18.60  V 10“'^ 

8.65  X 10 

2.14 

5.95  y 10'’^ 

4.14  X ur'^ 

1.96 

M.50  X 10*'^ 

9.60  X 10*'^ 

2.89 

10.90  X 10"’^ 

4.31  X 10 

2.81 

7.84  X 10“’^ 

5.30  X 10*'^ 

1.57 

8.00  X 10“^ 

4.65  X 10"^ 

1.85 

18.90  > 10"'^ 

6.80  X ur^ 

1.94 

7.25  X 10~^ 

5.80  X 10“'^ 

2.87 

25.4  X 10~^ 

9.2  10~’^ 

2.73 

9.0  X 10  ’ 

8.08  X 10 

2.81 

13.1  X 10"^ 

4.32  > 10~'^ 

1.90 

6.95  X 10  ^ 

8.26  X 10”^ 

2.16 

8.70  ' lO"^ 

4.81  X lO’’^ 

2.13 

26.8  X 10“^ 

16.00  V 10"'^ 

2.51 

14.50  X 10“^ 

10.72  X 10“'^ 

2.18 

22.8  X 10”^ 

17.80  X 10“’ 

o r.Q 

4.50  X 10  ^ 

3.63  X 10“’ 

2.22 

8.85  X 10"'^ 

6.91  X 10“’ 

2.16 

6.95  X 10“'^ 

4.87  10“’ 

2.61 

15.2  X 10”'^ 

11.60  X 10“’ 

2.05 

18. 1 X 10"'^ 

14.80  X 10  ’ 

2.70 

.5.22  X 10”'^ 

8. 18  X 10“’ 

1.56 

8.55  y 10~^ 

4.84  X 10“’ 

2.66 

12.8  X lO”’ 

9.15  X -0“’ 

2.37 

12.9  X 10“'^ 

8.55  X 10“’ 
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Table  X (cont’d). 


T 

(■  Cj 

\f 

(kllz) 

N 

Liisor 
fxiinction 
(dl^  III) 

383.8 

(1.03 

3.33 

:K) 

.399.3 

4.  13 

3.07 

ai 

(il3.9 

(1.(13 

3.88 

af) 

383.  (i 

3.98 

3.31 

.•m 

(i08. 1 

0.37 

3.88 

M 

3(13.8 

3.74 

1.8(1 

•a.'') 

(134.4 

(1.34 

3. 33 

;M.r> 

(109.3 

(1.39 

3.93 
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3.  .38 

■M 

331.(1 

3.30 

3.33 

330.3 

3.34 

1.10 

-3-1 
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3,07 
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:i-t 

39  3. 3 
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3.  Ml 

.at) 

370.8 
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300.  1 
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1 .00 
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3.. 30 

33 
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3.79 

1.90 

37.3 
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3.83 

1.37 

31) 

338,9 

3.7(1 

3.07 

-30 

303.9 

3.72 

2.31 

-30 

183.1 

1.89 

3.44 

-30.3 

1.34.8 

1.00 

1.94 

-29 

141.3 

1.40 

1.94 

-39 

108.7 

1.74 

1.9(1 

39 

331.3 

3.39 

1.90 

39 

170.0 

1.70 

1.90 

-30 

383.7 

3 93 

1.0(1 

30 

313.3 

3.30 

1.48 

N (l«  III 

till 

1 

,\  (IB 'm 

3.09 

3.04 

10 

1.03  X 10 

1 .99 

3.44 

10 

3.03  • 10 

3.31 

3.33 

10  *’ 

1.83  X 10"‘’ 

1.73 

0.3(1 

10 
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(1.73 

10 
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3.01 
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10 
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1.04  • 10  '1 
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10 
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1.91  X lo'*’ 

2.13 
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o 

1 

c: 
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1.19 
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. 10“*^ 
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1.33 
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■ 10 
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4.1(1 

■ 10 
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3.03 

70 

• 10  ■*; 

3.80  10 

1.47 

3-98 
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3. 10  X 10 

1.8 

3.92 

X 10^« 

2.3  X 10  ^ 

3.78 

7.93 

-•  10~‘’ 

3.84  X 10 

1.18 

10.1 

■ 10 

4.37  X 10 

0.773 

12.0 

^ lo  '^ 

3.16  X 10^’ 

0.833 

9.(1 

X 10^ 

4.93  X 10'“^’ 

0.734 

10.4 

X 10'“*’ 

3.33  • 10“® 
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3.08 

X lo””*’ 

2.39  X 10~^ 

1 .22 

3.93 

V 
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0.9 

9.33 

• 10"‘* 

4.72  > 10 

1.70 

4.38 

'■  10 
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1.49 

4.33 

X 10 

2.83  X 10 

fof;  impactor  mimber-density  measures  and.  as  will  be  seen  in  the  next  section,  the  eleotroniagnetic 
measurements  correlate  very  well  with  the  laser  extinction  level.  Thus,  the  best  possible  ice  log 
characterization  is  obtained  by  the  impactor  measurements  together  with  the  laser  extinction  level. 

Electromagnetic  measurements  and  applications  to  one-way 
propagation  through  ice  fog 

Sixty-nine  different  electromagnetic  measiirenients  of  the  complex  permittivity  of  a fog  weie 
made  covering  the  temperature  range  -30"  to  -4S"C.  Recorded  simultaneously  with  tlie  cavity  fre- 
queiicy  shift  and  amplitude  change  was  the  laser  extinction,  which  was  put  on  chart  paper  for  con- 
venience in  future  analysis.  Chart  speed  in  all  cases  was  in. /min. 

Since  it  is  obvious  that  an  accurate  measure  of  the  fog  number-density  distribution  is  at  present 
impossible  to  obtain,  it  is  necessarv  to  refer  all  electromagnetic  measurements  to  some  other  refer- 
ence fog  level.  One  such  convenient  level  is  the  6328  A extinction  record.  The  laser  extinction, 
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in  (IB/m,  was  used  to  normal  i/e  the  measured  N-unit  change  ot  the  cavity  refractive  index  and  the 
extinction  coefficient  x^^of  the  radio  wave  in  the  fo^.  That  is,  the  units  for  these  quantities  ate 
N-units/(dB/m  at  6328  ^ and  cnrV(dB/ra  of  6328  A).  The  data  from  these  measurements  are 
presented  in  Table  X,  Figures  45  and  46  show  plots  of  the  normalized  N-unit  change  in  cavity  fre- 
quency and  extinction  coefficients  as  a function  of  temperature,  For  Figure  45,  the  standard  devia- 
tion error  o associated  with  the  spread  in  the  d^ta  points  decreases  with  decreasing  temperature, 
dropping  from  a = 0.642  N-units/(dB/m  at  6328  A)  at  -30°C  to  a = 0,182  N-units/(dB/m  at  6328  A) 
at  -46°C.  The  standard  deviation  error  « was  computed  from  the  relationship 

»■  t 

where  x is  the  mean  valrie  of  y or  the  refractive  index  found  from  averaging  these  experimental 
measurements.  The  mean  values  of  y and  N at  each  temperature  interval  of  interest  were  found  by 
assuming  that  the  measurements  appropriate  to  these  intervals  applied  to  some  mean  temperature 
appropriate  to  that  interval.  The  smaller  standard  deviation  error  of  the  colder  fogs  is  consistent, 
with  the  fact  that  the  colder  ice  fogs  were  much  more  easily  reproduced  than  their  warmer  counter- 
parts. However,  the  refractive  index  of  these  ice  fogs  is  relatively  constant  over  the  temperature 
interval  -30°  fo  -47°C. 

The  plot  of  the  extinction  coefficient,  in  cm“^  (Fig.  46).  shows  an  interesting  rise  associated 
with  the  warmer  ice  fogs.  Since  -30°  fo  -35°C  is  rather  critical  to  ice  fog  formation,  it  is  reason- 
able to  suppose  that  some  percentage  of  the  ice  fog  is  in  the  liquid  state;  either  some  small  per- 
centage of  the  droplets  could  be  liquid,  or  some  ice  crystals  could  have  a liquid  coating.  By  the 
time  that  the  temperature  has  dropped  to  ~39°C  or  lower  the  extinction  coefficient  has  dropped  by  a 
factor  of  4 or  5 to  about  5 x cm“^  (for  1 dB/m  of  laser  extinction).  A part  of  the  error  asso- 
ciated with  the  measure  of  y of  the  colder  ice  fogs  is  the  difficulty  in  accurately  reading  small 
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Figure  46.  Normalized  ice  fog  extinction  coefficient 
measurements  at  97  GHz. 

changes  in  the  cavity  amplitude  response.  F'or  this  reason,  the  standard  deviation  error  bars  are 
larger  for  the  colder  ice  fogs. 

Figures  47-56  are  traces  of  the  oscilloscope  photographs  made  of  the  cavity  response  to  these 
fogs  over  tlie  temperature  interval  -30”  lo  48”C.  Two  pictures  a’  each  temperature  interval  of 
interest  were  chosen  as  representative  samples.  Accon  panying  each  photograph  tracing  is  the 
section  of  the  strip  cliart  with  the  laser  level  at  the  time,  the  photograph  was  taken.  The  original 
photographs  were  not  of  high  enough  quality  and  contrast  to  reproduce  directly,  which  is  why  tracings 
of  these  pictures  are  included  rather  than  reproductions  of  the  original  photographs.  The  high  corre- 
lation between  the  measured  laser  extinction  and  the  electromagnetic  tiieasurements  is  well  demoii- 
st>-ated  by  Figures  45  and  46. 

The  extinction  coefficient  y,  in  cm~*.  and  the  change  in  the  refractive  index  An  were  directly 
measured.  Knowing  these  two  quantities,  the  propagation  constant  fi  in  the  ice  fog  can  be  computed 
from  the  relation 
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Figure  47.  Sample  cavity  response  dtie  to  ice 
foff.  T 3&  C. 
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Figure  4S.  Sample  cavity  response  due  to  i 
foB.  T 30  C. 
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This  expression  can  Ire  simplified  in  this  case,  and  to  a hinh  decree  of  accuracy  is  j;iven 
fi 

The  loss  tangent,  tan  8 t" U' , of  this  medium  can  Ire  readily  determined.  Writing 
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Figure  51.  Sample  cavity  response  due  to 
ice  fog.  T ^ -39°C. 
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Figure  52.  Sample  cavity  response  due  to 
ice  fog,  T = -39°C. 


and  thus 


s .a,-  = 


(n'f-(n"f 

Table  XI  lists  the  appropriate  propagation  parameters  which  characterize  the  propagation  of 
97-GHz  radio  waves  through  ice  fogs  over  the  temperature  interval  -29°  to  -46°C. 

At  the  conclusion  of  each  data  run.  the  fan  was  switched  on  to  dissipate  the  fog.  As  the  fog 
dissipated,  the  cavity  response  curve  was  observed  to  shift  back  to  its  original  unperturbed  state, 
which  is  a strong  argument  against  drift  and  mirror  accumulation  effects. 

With  these  results,  calculation  of  the  effects  of  an  ice  fog  on  97-GHz  radio  wave  propagation 
can  be  made.  In  terms  of  signal  propagation  effects,  it  is  more  appropriate  to  refer  to  a 1-km  path. 
The  attenuation  and  phase  shifts,  appropriate  for  the  temperature  interval  studied,  for  a 1-km  path 

O 

and  an  assumed  1-dB/m  extinction  of  6328  A light,  are  presented  in  Table  Xll.  As  pointed  out 
earlier,  the  electromagnetic  measurements  are  well  correlated  with  the  reference  helium-neon  laser 
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Figure  53.  Sample  cavity  response  due 
to  ice  fog,  T -43^C, 
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Figure  54.  Sample  cavity  response  due 
to  ice  fog,  T -43'^C. 


Table  XI.  Propagation  constan‘.s  of  Ice  fogs  at  97  GHz  for 
1 dB/m  e'Uinction  of  6328  A red  light. 


iV-uni  ts 


Mean 

extinction 

coeff 


nad/m ) 


1.31  + 0.64 
1.96  + 0.46 
2.36  + 0.36 
2.04  + 0.28 
1.69  + 0.18 


3.83  + 1.11  X 10 
2.21  + 0.7  X 10 
8.08  1 4.68  X 10 
6.12  + 1.79  X 10 
.5.0  i 1.71  X 10 


1.08 

X 

10 

1.09 

X 

10‘ 

3.98 

X 

10' 

3.01 

X 

10‘ 

2.46 

X 
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Figure  55.  Sample  cavity  response  due  to  ice 
fog.  T = -46. 5° C. 
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Figure  56.  Sample  cavity  response  due  to  ice  tog, 
r = -46.5°C. 


Table  XII.  One-way  radio  wave  propagation  effects 
(for  1 dB/m  extinction  of  6328  A red  light). 


T 

<°C) 

AtlenuaUoa 

(dB/km) 

DUtereatial 
phase  shift 
(fi-  fi^)-  z (rad/ km) 

-29.3 

1.66  + .48 

2.666 

-34.4 

0.96  + .31 

3.989 

-38.3 

0.35  1 0.2 

4.803 

-42.8 

0.27  + ,J8 

4.151 

-46.2 

0.22  + .08 

3.439 
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level  for  extinctionn  of  2 dB/m  or  less.  Therefore,  it  is  appropriate  to  scale  to  propagation 
parameters  from  this  table  according  to  the  measured  laser  level  so  long  as  this  limitation  is  ob- 
served. 

Comparison  of  computed  and  measured  values 

Change  ol  index  o(  refraction  due  to  the  ice  fog.  ThOygeneral  results  obtained  in  the  previous 
section  may  be  compared  with  the  theoretical  prediction  of  the  Van  de  Hulst  expression 

n'  = 1 + ImlS(O)  I. 

The  ice  fog  impactor  measurements  of  Chapter  Vll  emphasized  that  the  number-densities  re- 
ported were  a lower  limit  to  the  ;rue  number-density  since  the  collection  efficiency  of  the  impactor 
was  unknown  and  the  smaller  particles  were  poorly  sampled.  As  indicated  in  Chapter  Vll,  the 
craters  measured  under  the  microscope  were  taken  to  be  indicative  of  the  particle  size  making  them, 
but  this  may  not  be  accurate.  For  example,  a solid  spherical  particle  colliding  with  a thin  magne- 
sium-oxide coated  glass  surface  would  not  necessarily  make  a crater  the  size  of  its  own  diameter, 
but  could  in  fact  make  a much  smaller  area  proportional  to  the  amount  of  magnesium  oxide  crushed 
down  to  the  point  of  contact  of  the  spherical  particle  surface  with  the  glass. 

Both  the  literature  and  private  communication  with  some  authorities  in  this  area  have  indicated 
that  no  method  or  technique  exists  for  making  accurate,  reliable  ice  fog  number-density  measure- 
ments. The  measurements  made  on  these  ice  fogs  indicated  the  existence  of  number-densities  of 
from  100  to  about  300  particles/cm3,  with  the  apparent  diameter  of  most  particles  less  than  about 
15  ^L. 

In  any  case,  it  is  reasonable  to  take  200  particles/cm^  of  20  g diameter,  which  gives  an  N-unit 
refractive  index  change  of 

(n  - 1)  X 10®  = 200X10®)  _ 3 2 gQ  ^ ^ 

8.43  10^  2 


s 0.52  N-units. 


The  expected  change  in  the  index  of  refraction  is  quite  sensitive  to  the  ice  particle  size,  as  pre- 
viously stated.  This  can  best  be  emphasized  by  making  one  more  brief  calculation  for  30-fi-diam 
particles  with  the  same  density.  Here, 

(n  - 1)  X 10®  = 2”<200)(10®)  3 ^ ^ ^ 

8.43  X 10^  ^ 


^ 1.73  N-units. 


These  calculations  assumed  small  (with  respect  to  a wavelength)  spherical  particles  of  Cv.mplex 
dielectric  constant  = 3.08  - / •004-  The  value  of  S(0),  appropriate  to  these  particle  diameters, 
is  3/2  times  the  value  of  bj""  listed  in  Table  V. 

This  is  not  meant  to  minimize  the  effective  contribution  of  the  smaller  particles.  Large  num- 
bers of  the  smaller  particles  are  undoubtedly  present,  and  their  effective  contribution,  which  varies 
directly  with  their  number  density,  is  likely  some  significant  fraction  of  the  overall  N-unit  measure- 
ment. 


,2  COUPI.EX  refractive  INDEX  OF  ICE  FOG  AT  A RADIO  WAVELENGTH  OF  3 MM 

For  ex™, pie.  10,000  parlleles  ol  6 p <li™e'ei  I*'  »'>""* 

index  of  refract ii'ii  of 


f.f)  1)  ■ Itl 


(j  (iJrr)  ' 1 '■  10^  10*^  X ^ y g o - iri  8 X in*’ 

8.4:i  ' lO'’ 


■ X - y 6.2  ■-  10  “ X,  10® 
2 


^ 0.7  N-unil 

„„c,  xe.en,l  u.,ul.s  ..co  N.„„l.  co.,,0  e.»,ly  co„,e  1,0, „ „.e  more  eeslly  Oe.ecled  l.rpe,  parnCex  In 

the  size  range  16-25  /c  ,, 

Radio  yaye  ions  d«e  ,»  .«e  <»«■  The  expreen.o,,  lor  enerpy  loss  due  lo  an  assemblape  ol  s,.,all 

scatterers,  as  given,  is 


okn  ' , 


, 2n 
where  k — 


and  n 


2n^  X- 
" r 


RelS(0)l. 


;r,.”==rrz;-,=-r;r.;;= 

theoretical  loss  of 


> 10'*  • 10*'  X ?.  . 7 04  • 10  **  • 10~"  ' 

X'  ,,  9 

4.5  X 10*^ 

2.95  ' 10'® 

mid  for  vnn  oarticles  30  n i»  diameter,  the  loss  would  be 


(4?7)(2)  ■ 10~  " 10*’ 
4.5  X 10® 


- ^ 8.9  ■ 10  ® X 10 

2 


J.b 


10  8 cm 


r,'-;4T=T, 

These  calculations  were  apaia  n,ade  for  a dicleciric  constaul  ol  3.08  - I «m- 

T1  cse  calculallons  can  1«  con,parod  will,  ihe  predicled  loss  tor  a water  fop.  In  M„s  case, 
f li  - M5.  diid  for  6-/r-diam  water  spheres 
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Therefore  the  ratio 
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water  ^ 7 ^ 

— 250 

(hj“')ic(i  7x10“^^ 

and  tliiis  tlie  expected  loss  for  a water  fog  would  be  on  the  order  of  250  times  that  for  an  ice  fog  of 
the  same  mimber-density.  The  theoretical  and  experimental  measurements  are  in  far  closer  agree- 
ment than  this  (less  than  a factor  of  lO)  for  the  ice  fog,  and,  as  mentioned  in  the  preceding  section, 
this  discrepancy  certainly  reflects  in  part  the  uncertainty  in  the  particle  size/nuniber-deiisity 
distribution.  Of  course,  some  small  fractional  part  of  the  fog  particles  could  still  be  liquid  at  these 
temperatures,  which  could  account  for  this  discrepancy. 

Experimental  and  theoretical  relationship  jf  n'  with  n”.  One  further  bit  of  information  can  be 
extracted  from  these  data.  From  the  Vtui  de  Hulst  expressions  for  the  complex  refractive  index  of 
the  medium, 

n'  1 + Im  |S(0)l 

n"  ^RelSfO)! 
k3 


S(0)  ~^(2n^  l)(a«  . 

n-i 

where  all  factors  are  the  same  as  previously  described  in  Chapter  IV.  Then  the  ra'io 
n"/(n'  - 1)  io“6 

is  independent  of  the  actual  fog  particle  number  and  distribution,  and  only  dependent  upon  the 
complex  dielectric  constant  of  the  individual  particles.  This  can  be  readily  seen  for  this  case 
where  the  wavelength  is  much  larger  than  the  particle  diameter  (Rayleigh  approximation),  where 

"•O 

n-i  ^ 

and  whtre 


' 45 

- first  magnetic  mode  in  the  sphere 


first  electric  mode  in  the  sphere 
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where  p r k^a 


radius  of  dielectric  spliere 

m = complex  refractive  index  of  tlie  siiiull  spheres. 

i dominates  a^.  hence  the  ratio  RelS(0)!/Im Urn S(0)1  is  the  ratio  of  the  real  to  the  inuiK  nary  part 
f j,  and  hence  tlie  sphere  size  dependence  drops  out. 

Ihe  Fabry-Perot  cavity  measures  the  extinction  coefficient  associated  with  the  foe  losses 
S...OC.  or  ,,„s  w„vele„„„.  ,,„o  parMclo  abso,p,i„„  d„,„i„„,es  »ca„e,i;,!  stS  r„ 
all  pidctital  piuposes  extinction  absorption,  and 

\ ^ 2a 

where  „ is  a field  quantity  and  ^ is  the  extinction  measured  for  the  power  loss  due  to  the  ice  foe 
Another  form,  as  given  eailier.  for  the  complex  index  of  refraction  is 


Therefore, 
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and  lor  a measured  solid  ice  complex  dielectric  constant  of  i 3.1  - / .004. 


S(0) 

\m-  f 2/ 


t 0.463  X 10  ^ 1 j 0,411 


COMPLEX  REFRACTIVE  INDEX  OF  ICE  FOG  AT  A RADIO  WAVELENGTH  OF  3 MM 


75 


then 


^_iS(0)] 

Iru  (S(O)l 


1.13  ^ 10 


The  ratio  n /(o'  i)  lo'^  was  based  upon  data  taken  froiii  the  low  temperature  ice  fo^s; 

5 ^ 10  ^ ' is  well  within  the  possible  error  for  the  thre  ? lower  temperature  fo^s.  A y of  5 

10  ^ ciir*  was  used  since  these  lower  temperature  fogs  are  stable  in  that  they  are  considered  be- 
low the  transition  temperature  for  spontaneous  droplet  freezing.  In  computing  tlie  ratio  of 


Re  |S(0)1 
Im  IS(0)t 


the  mean  dielectric  constant  found  earlier  for  slab  ice  was  used.  Tlie  match  between  the  measured 
value  of 


n"/(n'  - 1)  10'  ^ 

and  the  computed  quantity 

Re  IS(0)I 
Im  |S(0)I 

is  fairly  good,  and  these  results  are  not  unexpected  in  light  of  the  results  of  the  preceding  section. 
With  the  error  bars  involved  on  the  measured  quantities,  going  into  the  above  ratios,  it  is  possible 
to  obtain  a near-perfect  agreement,  although  this  is  not  the  important  point  here.  The  point  is  that 
the  Fabry-Peroi  cavity  measurements  which  can  be  related  to  the  complex  dielectric  constant  of 
the  ice  tog  part'.cles  are  in  substantial  agreement  with  previously  made  slab  ice  complex  dielectric 
constant  measu.ements. 


Vni.  SUMMARY  AND  CONCLUSIONS 

Tlie  main  thrust  of  this  work  was  to  make  measurements  of  the  index  of  refraction  and  signal 
loss  associated  with  propagation  of  3-mm  radio  waves  through  an  arctic  ice  log.  This  work  was 
originaP>  'notivated  try  the  planned  operation  of  3-mm  pulse  radars  on  surface  effects  vehicles  to 
be  operaieo  ove:  arctic  terrain.  These  radars  scan  horizontally  at  near  grazing  angles  of  incidence 
and  must  be  capable  of  seeing  objects  through  low  temperature  fogs  soon  enough  for  evasive  action 
to  be  taken.  The  low  temperature  ice  fog  is  a relatively  common  meteorological  occurrence  in  the 
Arctic  and  Antarctic,  and  almost  nothing  was  luiown  of  its  effect  on  3-mm  radio  waves.  These 
measurements  are  thus  a “first  ever”  attempt  to  characterize  the  3-mm  radio  wave  propagation  effects 
of  an  ice  fog. 

It  was  also  of  specific  interest  to  make  measurements  of  the  complex  dielectric  constant  of 
slab  ice  at  97  GHz  to  compare  these  measurements  with  the  high  frequency  limiting  values  of 
typically  used  in  the  Debye  equations  and  found  in  the  literature  for  frequencies  up  to  about  25  GHz. 
The  complex  dielectric  constant  of  pure  ice  was  found  to  be  (3.08  + 0.2)  - 1(3.99  + 1.50)  x 10~*^- 
The  real  part,  3.08,  is  in  excellent  agreement  with  the  limiting  values  reported  in  the  literature  for 
frequencies  up  to  24  GHz,  but  the  imaginary  part  is  slightly  higher  at  97  GHz  than  the  reported 
values  for  lower  frequencies.  Measurements  of  samples  of  the  ice  prepared  from  water  samples 
taken  from  a local  lake  gave  sligtitly  higher  values  for  the  complex  dielectric  constant,  but  the  error 
bars  asbociated  with  these  measurements  are  all  overlapping. 
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The  teeliiiique  used  to  make  tliese  measurements  involved  a quasi-optical  free-spaee  approach 
commonly  described  in  the  literature  and  suc.cessfully  used  at  lower  frequencies.  Prior  to  usinK 
this  technique  directly  on  makinf?  the  ice  measurements  an  extensive  sequence  of  measurements  was 
made  on  materials  of  known  permittivity.  One  such  material  with  a well  behaved  and  known  permit- 
tivity is  polystyrene.  This  calihration  sequence  was  done  primarily  to  estalilish  the  necessary 
erne  and  skill  prior  to  making  the  ice  measurements. 

The  ice  fog  measurements  were  made  with  a F'abry-Perot  cavity,  a device  commonly  used  in 
optical  interferometry.  'I'he  Fahry-Perot  resonator  has  been  adapted  to  millimeter  wave  usage,  and 
ample  precedent  for  its  use  exists  ni  the  literature.  The  Fahry-Perot  resonator  built  for  this  experi- 
ment is  an  opiMi-sided  semi-confocal  resonator  (the  fields  are  established  between  one  spherical 
and  one  flat  mirror)  with  a measured,  loaded  Q of  140,000  at  an  operating  wavelength  of  3,09  nim. 

The  semi-conl'(X!al  imvity  was  selected  for  use  primarily  for  the  convenience  of  its  physical  size 
and  its  ease  of  mechanical  construction.  It  can  be  demonstrated  that  the  change  in  the  index  of 
refraction  of  a medium  is  measured  by  the  change  in  the  cavity  resonant  frequency,  and  the  change 
in  the  amplitude  of  the  cavity  response  is  related  to  the  losses  of  the  medium. 

A method  employing  a ptiase-locKing  technique  was  used  to  lucktlie  klystron  to  a higli  multiple 
of  a stabilized  reference  1.7-GHz  oscillator.  The  difference  between  this  reference  frequency  and 
the  klystron  ftetitteHr;y . the  W.  war  cw.p'ttttd  with  oulpw  o(  a lx  in  a |iui*.i  Irf-kitv 

system,  and  the  dc  error  output  was  used  to  steer  the  klystron  in  such  a manner  that  its  frequency 
deviation  at  07  GHz  followed  the  frequency  deviation  of  the  sweep  generator.  Thus  the  klystron 
trequency  was  locked  and  swept  i M over  the  Fabry-Perut  cavity  ftequeiicy.  Tlie  output  of  the 
cavity  was  sampled  tlirough  a square-law  detector  loosely  coupled  to  the  cavity  and  displayed  as 
the  ordinate  oti  an  oscilloscope.  Marker  lieats,  whose  frequencies  could  be  measured  with  high 
precision,  were  also  displayed  superimposed  on  the  horizontal  oscilloscope  sweep,  and  these  were 
used  to  measure  the  cavity  fiequeticy  shift  between  conditions  of  ice  fog  and  no  fog  in  the  cavity. 
This  technique  of  sweeping  the  cavity  and  photographing  the  response  proved  especially  convenietit 
for  information  readout  and  storage. 

Ice  accumulation  on  the  mirror  surfaces  was  investigated  by  grazing  a laser  beam  across  the 
metal  surface  and  photographing  the  light  scattering  from  the  adhering  ice  particles.  Time  lapse 
photography  showed  that  this  accumulation  occurred  only  during  the  first  few  seconds  after  entry  of 
the  warm  steam  fog.  and  that  the  total  accumulation  was  small  when  compared  with  the  light  scatter- 
ed by  a thin  3-/i  polystyrene  dust  coat.  Furthermore,  it  was  seen  that  if  a thin,  low  loss  dielectric 
material  covered  the  mirrors,  there  was  no  measurable  change  in  cavity  response.  When  the  environ- 
mental chamber  was  shut  down  and  quickly  opened,  a heavy  white  frost  coating  was  quickly  formed 
on  all  interior  metal  parts  from  the  inriishing  warm,  moist  air.  This  had  no  effect  on  the  cavity 
response,  but  just  as  soon  as  this  frost  coat  started  to  liquefj  the  response  changed  radically. 
Additional  confirmation  of  mirror  accumulation  effects  came  when  the  interior  fan  was  restarted  after 
fog  electromagnetic  measurements  hod  been  made.  As  the  fog  was  dissipated,  the  cavity  response 
curve  could  be  observed  to  shift  back  to  its  original  no-fog  response.  This  is  also  indicative  of  a 
iiiiiiinuiii.  amufliit  of  n.ecb4tiical  atrl  tlecKleal  tlrift  effect* . 

An  environmental  chamber  of  about  70  ft®  volume  was  constructed  of  Styrofoam  to  house  the 
ief.iHct,.t;ur.  evaporators  and  ttie  Falrty-J*eroi  cavity.  Waveguide  m lljjjugb  llw  rljitr.W  waJls 
to  connect  with  the  electronics.  Mirror  movement  was  accomplished  through  remote  drive  motors 
suitably  geared  down  so  as  to  move  drive  screws  in  three  mutually  orthogonal  directions  tc  control 
the  spherical  mirror.  'I'he  refrigeration  cumpreBbur,  a 3-hp,  3-pliase  unit,  was  ttsuf reeled  from  an 
old  piece  of  surplus  equipment.  To  get  to  t'  ' low  temperati  res  necessary  for  this  work  it  was 
necessary  to  use  a lower  boiling  point  refrigerant,  R-502,  Chamber  temperatures  as  low  as  -50  C 
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could  be  achieved  and  maintained.  To  minimize  vibration  effects  during  the  electromagnetic 
measurements  it  was  necessary  to  stop  both  the  fan  and  the  compressor.  Temperature  changes  in 
the  chamber  f’uring  this  time  were  rarely  observed  to  exceed  two  degrees. 

Initially  it  was  planned  to  correlate  fog  number-density  measurements  with  the  electromagnetic 
measurements  This  plan  hinged  upon  having  a reliable  method  for  obtaining  these  measurements, 
or  at  least  having  experienced  help  in  obtaining  these  data.  As  it  turned  out,  the  help  never 
materialized,  and  a borrowed  cascade  impactor  was  all  the  instrumentation  available.  This  impactor 
used  magnesium  oxide  coated  slides  upon  which  the  ice  fog  particles  would  “crater.”  leaving  pits 
the  size  of  the  fog  particles  to  be  counted.  Slides  made  of  the  colder  fogs  were  generally  fairly 
consistent;  much  more  difficulty  was  encounter“d  in  trying  to  obtain  any  degree  of  repeatability  with 
the  warmer  fogs.  Chi  many  of  the  slides  taken  wi»,i  this  impactor  there  were  very  faint  traces  of 
particles  less  than  about  4 /i.  These  were  too  faint  to  satisfactorily  count,  and  thus  these  measured 
number  densities  represent  only  some  kind  of  lower  limit  to  the  Uue  fog  number  densities.  During 
the  later  stage  of  this  experiment  a single-stage  silicone  oil  impactor  was  used  to  make  fog  density 
measurements.  This  device  too,  proved  to  be  unsatisfactory  and  inconclusive.  This  type  of 
impactor  should  only  be  used  in  an  environment  at  or  below  the  freezing  temperature.  As  it  turned 
out.  the  collected  particles  started  to  disappear  from  the  oil  base  collecting  medium  as  soon  as  the 
impactor  was  withdrawn  from  the  environmental  chamber.  By  the  time  these  collections  were 
photographed,  an  average  of  20  seconds  had  elapsed,  far  too  long  a time  for  any  accurate  analysis. 
However,  these  oil  base  impactor  measurements  did  show  a very  substantial  number  of  small 
particles,  and  surprisingly  there  appeared  to  be  a small  number  (30  or  so  on  one  particular  photo- 
graph) of  very  large  (>  25  n diameter)  particles  per  cubic  centimeter  present. 

There  is  also  some  question  about  the  accurate  particle  size  replication  of  the  magnesium-oxide 
craters  formed  by  the  impacting  ice  fog  particles  The  "craters”  measured  and  studied  .der  the 
microscope  may  not  have  had  diameters  equivalent  to  the  ice  particle  diameters,  but  in  fact  the 
crater  size  could  have  been  smaller  than  the  actual  particle  sizes.  A solid  spherical  particle,  for 
example,  would  crush  the  magnesium-oxide  coating  down  only  to  the  point  of  contact  with  the  under- 
lying  glass,  which  would  make  a "crater”  much  smaller  than  the  impacting  particle.  This  is  un- 
doubtedly responsible  for  some  of  the  discrepant;y  between  the  theoretical  and  experimental  results 
(as  noted  in  Chapter  VII),  as  the  magnesium-oxide  films  deposited  on  these  slides  tended  to  be 
quite  thin  vhich  enhances  this  effect. 

In  order  to  achieve  some  degree  of  uniformity  from  fog  to  fog.  a helium-neon  la.ser  beam  was 
shot  through  a 2-m  path  in  the  fog.  and  all  electromagnetic  and  impactor  measurements  were  made 
at  or  near  a laser  extinction  level  of  l.c  dB/m.  The  required  fog  density  was  achieved  simply  by 
inputting  and  mixing  the  fog  with  the  fan  until  tlie  required  extinction  level  was  readied.  W.iile 
somewhat  arbitrary,  this  amount  of  laser  txtinction  was  chosen  since  it  was  well  within  the  uynamic 
range  of  the  optical  detector,  and  was.  for  the  total  path,  a total  extinction  level  checked  and  cali- 
brated by  an  optical  filter  of  known  transmission  properties.  In  addition,  this  fog  density  provided 
for  a reasonable  amount  of  cavity  response  (frequency  and  amplitude  change).  Visually,  this  fog 
level  appeared  to  be  at  least  as  dense  as  any  warm  fog  ever  encountered  in  a southern  swamp  on  a 
cool,  calm  evening. 

Calculations  of  the  total  extinction  of  the  6328  A laser  red  light  were  made  and  compared  with 
the  measured  extinction.  These  calculations  used  as  input  data  the  nmasured  icjniber -density  and 
particle  sizes  of  these  ice  logs.  However,  the  calculated  extinction  was  between  and  '4 
(depending  upon  the  particular  number  density  distribution  used)  of  the  actual  measured  extinction. 
This  says  that  either  the  number  distribution  used  in  the  calculations  is  somewhat  wrong,  or  Mie's 
theory  has  been  incorrectly  applied  to  the  ice  fog  particles  (that  is.  the  particles  may  be  tiny 
prisms,  or  geometrical  columns  with  flat  faces).  Thi..  information  was  communicated  to  interested 
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groups  at  Jolins  Hopkins  Applied  Physics  Laboratory  and  to  the  Cold  Regions  Research  and 
Engineering  Laboratory.  This  finding  is  in  substantial  agreement  with  their  results  for  ice  fogs, 
and  has  not  yet  been  resolved. 

One  of  the  major  problems  yet  to  be  solved  is  how  to  accurately  measure  ice  fog  particle  num- 
ber densities.  One  group  at  White  Sands  made  a study  of  impactor  types  and  various  replicating 
media,  and  concluded  that  fog  drop  size  distributions  were  method-dependent.  They  also  concluded 
that  no  method  was  satisfactory  for  detecting  or  capturing  small  droplets  less  than  4 fi  in  diameter, 
and  that  “lack  of  correlation  between  number  of  captured  droplets  and  visibility  was  present  in 
even  the  most  car'^iul  determinations  employing  impactor  methods.  It  appears  that  either  small 
droplets  have  not  been  captured  and  replicated  or  that  factors  other  than  scattering  by  water  droplets 
contribute  heavily  to  tlie  deterioration  of  visibility  in  fog.” 

It  was  therefore  concluded  that  the  electromagnetic  measurements  could  not  be  referenced  to  the 
fog  impactor  measurements.  The  most  reasonable  reference  standard  for  these  electromagnetic 
measurements  was  taken  to  be  measured  extinction  of  a 6328  A helium-neon  laser.  The  laser  extinc- 
tion. supplemented  by  the  impactor  measurements  and  recorded  temperature,  thus  constitute  the  best 
meteorologic^al  de.scription  of  fog  which  could  be  obtained. 

By  being  as  consistent  as  possible  in  generating  these  fogs  over  the  temperature  interval  of 
interest,  the  measured  changes  in  the  index  of  refraction  ran  from  about  v,  to  2'A  N-units  for  i dB/m 
of  laser  extinction.  The  extinction  coefficient  dropped  from  a higli  of  3.8  10  ® cm'*  at  29‘^C  to 

8 X 10  ^ cm'*  at  -38'  C.  and  then  dropped  very  slowly  to  5 » 10"^  at  -46°C.  The  temperature 
interval  -30’  to  -36°C  is  somewhat  critical  to  ice  fog  formation.  A critical  literature  review  indi- 
cated that  natural  ice  fogs  almost  never  occur  at  temperatures  above  -30”C.  and  occur  much  more 
frequently  and  with  increasing  severity  (greater  reduction  of  visibility)  with  the  onset  of  lower 
temperatures. 

Theoretical  calculations  using  the  Van  de  Hulst  expressions  tend  to  confirm  the  experimcnial 
measurements.  Both  the  refractive  index  calculations  and  the  loss  calculations  tend  to  he  low  vbut 
within  a factor  of  10)  in  comparison  with  the  experimental  measurements.  TTie  experimental  measure- 
ments indicate  that  a higher  volume  of  ice  p.  tides  per  unit  volume  exists  in  the  cavity  than  is 
indicated  by  the  impactor  measurements,  which,  for  reasons  given  earlier,  are  thought  to  be  lower 
limits  to  the  true  number  densities.  If  the  fog  were  entirely  liquid,  the  theoretical  calculations 
indicate  that  the  loss,  for  example,  would  be  more  than  200  times  what  it  is  for  an  ice  fog.  It  is 
also  possible  that  some  small  fraction  of  the  fog  particles  are  liquid  or  have  a liquid  coat.  This, 
together  with  low  particle  samples,  could  account  for  the  d’screpancy  between  the  measured  and 
calculated  values  indicated  in  Chapter  VII. 

The  ratio  of  the  experimentally  measured  quantity 


n" 

- 1)  10'® 


can  be  compared  wuh  the  theoretical  ratio  obtained  for  a very  dilute  scattering  medium 

Re  IS(0)l 
Im  IS(0)l 

where  S(0)  is  formed  from  the  series  of  the  complex  amplitude  coefficients  associated  with  the 
electric  and  magnetic  modes  present  in  the  dielectric  spheres.  For  the  Rayleigh  limit  it  can  be 
shown  that  there  is  no  particle  volume  dependency  in  the  above  ratio.  Therefore  the  ratio  Re  IS(0)l/' 
Im  IS(0)l  depends  upon  the  complex  dielectric  constant  of  ice  used  in  the  computation  of  S(0).  When 
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each  of  those  ratios  is  computed,  in  one  ease  from  a knowledge  of  the  complex  dielectric  constant 
ice  and  in  the  other  case  using  the  measured  cavity  extinction  and  index  of  refraction  they  are 
less  than  a factor  of  six  from  each  other. 

Therefore  the  huger  extinction  (aiefficienls  and  their  rapid  drop-off  above  about  -36'  C are  not 
surprising,  sincte  some  small  fraction  of  the  ice  fog  may  still  be  liquid,  either  existing  as  i.solated 
liquid  diops  or  as  a liquid  coating  on  some  ice  crystals. 

Loth  the  refractive  index  measurements  and  the  extinction  coefficients  coirelated  very  well 
with  the  measured  laser  extinction.  While  ,t  is  true  that  the  laser  is  sensitive  to  small  particle 
< oneentratious.  and  the  cavity  is  sensitive  to  larger  particles  (since  its  response  depends  upon 
tie  VO  nnie  of  ice  in  its  fields),  the  particle  distribution  is  apparently  not  changing  radically  over 
the  temi^ratnre  range  -29"  to  -46  C since  there  ,s  both  a uniform  warm  fog  input  at  all  temperatures 
and  good  correlation  with  the  measured  electromagnetic  parameters. 

Since  the  impactor  measurements  of  these  fogs  are  quite  similar  to  the  measurements  of  estab- 
oc  mr^inffoes  '!lv 

g g.  , .x(  ept  in  density.  But  bet^ause  of  the  good  correlation  with  tlie  measured  lielium- 
ueon  laser  extinction,  it  is  appropriate  to  scale  lliese  measured  values  of  y and  N down  to  airive 
at  ji  expected  vaim*  for  97  GHz  radio  wave  extinction  and  phase  sliift  in  ai  ice  fog  of  lesser  den- 
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APP«5NDIX 

Tho  derivation  of  the  jiropagating  E-field  through  a four-boundary  three- 


layer  medium  is  as  follows; 
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y Tl^  L 3 3 
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R egion  4 : 


F - A o f B e 
X 4 4 


H = A Ta. 


y ! 4 


B e 
4 


Region  5: 


E = A .• 
X 5 


-j03^ 


1 r _ 


H - — j A_e 

y Tl.  L 5 


Match  Boundary  Conditions 
at  1-2  boundary,  z = Z = 0 


E ; 

X 


A,  =A^4B^ 


H : 

y 


A - B ' =-i 

LI  1 _j  71. 


I _ 1 f 

o \ ^2j 


at  2-3  boundary,  z = d 


E : 

X 


-r^d  r^d 

A c + B e ^ A,  + B, 
^ ^ 3 3 


H : 

y 


1 ^ "^2*^  ^2'^'!  1 


J Tl, 


Z + d 4 L 


z > Z -f  2d  4 L 


A3-B3 


at  3-4  boundary,  z = d + L 


E ; 

X 


-r^(L)  r (L) 

A3  e + B e = A + B^ 

J 3 44 


H : 

y 


r[ 


-r,(L) 


A 3-'  '■3(0  . 

3 = ■ ®3  ' ] = ii7  L\  - 


SuiasmUtim 


Z + 2d  + L 


ikiiiMiiilidijIUiiili 
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at  4-5  boundary,  z = 2^  + L 


A , e + B , c = A5 
4 4 


-V 


4 •!  1 


1 r 4 4 ■'  1 

— A e - B , e U — A, 
L 4 4 J Tl^  5 


T1  = T1  - Ti 
2 4 g 


Solve  8 equations  in  7 unknowns  simultaneously 


qd  -r^d 

'’■1  " ' -'b  ■ ^4 


-r.ci  -2i’ci 

4 4 

B = A e - A e 
4 5 4 


r ‘^4"^  ■^4'^  ; ^d 

Ke  ^ - ,A  ^ - A ..  ScM.A  A. 


4 


r ■''4''  -Vn 

A e -A  +A,e  = — 

L 4 5 4 Ti 


■ ^_4 
- 11. 


■-  'Vi  n 

! = LC"‘]'^e 


A 4^1a 

4 2 L 'n5  J 5 
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-r  d 


1.:^  2 Tl.  J 5 


A -I’d  T1 


[-c] 


-r  L 


®3  " = - Aj  e 3 


’^3^  -2r,L 

Bj  = A4  e + B4  a . Aj  e 


-^3^  rL  T1 

^3®  -B  e =~rA-Bl 

^ 3 \l  4] 


A '^3^  r "1^3^  -2F  L f^L 

^3  ® ■ [^4  ® ® - A e ^ "Je  " 


A ‘^3^  _ 

® - A - B + A e ^ =-l  1 


4 4 3 


i - S/ 

L 4 4j 


-BjL  Tlj 


r 
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A_  r,d  r L r d . T1  ^ T1  , T) 

. = — e ^ e c*  + — )(  1 + — n 1 ) 

^2  8 ^ V T1  71  yV  T.  y 


V Tl- 


R ® '-  V Tiy^  ti^^v  r,^y 


3)0  3)r.  3) 

8 ‘'V  n^yV  Ti^yv  Ti^y 


A r.d  -r,L  -r  d ti  ti 

— e "■  e e ^ ^1  - -^  )(  1 + -^Y 1 - 

8 ^ . TlVV  71,A  Ti3 


^z~  2^  L Y J L 4 ^ ^ V i\^y 


^ "3^  -Vy,  Vi 

^ 3 ® ( -nJ^  11^- J 


I '^'^1  ^Yi.- 


3)^.  3)1 

4 \y^ 


A^  -r,d  FjL  r,d  y3^^,,  V 


3^-^^'‘e-\'^^V.3)o3)o3) 

8 V T|jA  Tig.'' 


A5  -V  -F3L  V y^, 

*T'  " ' I y. 
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-r^d  .FjL  -r^d 

+ -|-  e 


e 


(i  +!iy,  .V, 


\ =T6 


Aj  Tjd  TjL  r.d 


e e e 


^ tiY'-  11,  a Ti,y'. 


A F-dpL-rd  T1  T1  n n 

5 2 3 4A  1\/'  ?"\/' 


16*"  V*  ■*■  - — jfl  - -- Vl  + — ^ 

X 11,  y . T1  yV,  T|  T]  ^ 


A_  F,d  -F,L  ~r.d  T1  T1  Tl  Ti 


A5  -r^d  r,L  r d^  ti  ti  ,1  ^ 


■^c  ^ ^ ^ ^1  Tl  Tl  Tl 

^0-Ax,.J)C-X)0-^ 


H5) 


*t6  ' ' ' -is^)0 +Xa' ~¥> 


nyv  Hj 


A3  -r^d  -F3L  -r^d 
■"T6  ® e . e 


0-tX-4x-XX‘-S 


Therefore 


A5  ~ I- '^2  ■‘'  ■•■  ■^2)  ^4)  ® 


(F2d  + r^L  + r^d) 
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+ i\  + - Tl^)  e 


+ {\  + \ ){T\^  - - T1^)(T1^  f Tl^)  e 


+ (TI2  + TIjHTI^  - Tl^HTl^  4 - Tl^)  e 


+ (\  - "niXTI^  - J\^m^  4 T]^)(Tl5 
+ i\  ■ T1j)(T1^  - TI^HTI^  - Tl^HTl^  - Tl^)  e 


(r2d4F3L-  r^d) 
(r^d-  r3L4r^d) 
(r^d-r^L-r^d) 
(-r^d  4 r^L  4 r^d) 
(-r^d  4 r^L  - r^d) 


+ - TIjXTl.  4 Tl^HTl^  - TI3HTI3  4 T^)  e 


(-r^d  - r^L  4 r^d) 


(-r  d - r L - r d) 

+ (TI2  - TljHTl^  4 n^XTl^  4 Tl3)(Tl5  - T1^)  e ^ ^ J 

Now,  thickness  of  medium  2 = thickness  of  medium  4,  T|  =71  = 71  . and 

2 4 'g 


r = r . = r 

2 4 g 


1 1 F 

x ' 7TX L*\  ^ ’'3”^  “ 

5 16  71  71,71,,  ^ JggJSg 

g 3 '5 


(2F  d 4 r,L) 
g 3 


MTIg  + 71j)(7l3  4 71g)(71^  . . Tl^)  e 

(2r  d - r.L) 

+ (Tig  + 71i)(7l3  - 71g)(71g  - 713X715  + Tig)  e ® 
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+ \ - ^g)(Tlg  + T|3)(T1^  - e 


(-T3L) 


+ (Tig  - T)j)(Tl3  ■ ^3X^5  + T1  ) e 


(r,L) 


(-2r  d f r L) 

+ (Tig  - TI1XTI3  - n^HTlg  - Tl^Hiq^  - T1  ) e g ^ 


(-r,L) 


+ (^g  - T1j)(n3  " - ^3>(Tl5  + Tig)  ' 


^ (-2r  d - r L) 

+ (Hg  - TliHTlj  + IgllUg  + Djjdlj  . -n^)  e S 3 

which  can  be  simplified  somewha*.  for 

^ ^5  ■ either  side  of  the  dielectric  sandwich) 

Break  this  up  into  four  constants . C,  D.  E and  F defined  as  follows: 


C = 


16T1  TIjH^ 


2r  d 

g 


-2r  d 


+ 31,1(313  - 'lg){Tlg  - ® ] 


D 


•Mg-ljTls  ^ « 


[(\  + ''l|)(Hj-  ’IgX'lg  - lljKTlj  + TIJ 


2r  d 

. g 


-2r  d 

g 


+ (Tig  - Tip(7l3  + T)g)(T]g  + 7l3)(7l5  - Tig)  e « ] 


E = 


^ '•"”3  ^ - ^3><\  - \) 

s 
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‘15 


Thus,  at  last, 


A TjL  -rjL 

— = (C  T E)  e 4 (D  + F|  e 

2r  d -2r  d 

g g 

At  millimeter  wavelengths  and  for  glass  holders,  e » e . and 

-2r  d 

thus  the  terms  involving  c ® in  C and  D may  be  neglected.  Thus, 


A,  r,L  -r,L_  zr  d 

= -2'  ’ ^ 


where  Kj  and  contain  only  the  dominant  terms  of  (C  + E)  and  (D  + F), 


respectively.  Thus, 


r 2 ^3^ 
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W7 


-zr  d . . 

and  the  factor  e ® attenuates  both  the  amplitude  and  phase  variation. 
Thus  it  is  important  that  the  glass  dielectric  holding  material  be  as  thin 
as  possible.  From  this  expression  for  E(5)  it  is  clear  that  the  measured 
ph'.se  variation  would  not  be  a linear  function  of  (water)  thickness,  but 
would  rather  have  some  perv  iic.  damped  variation  about  a straight  line 
best  fitting  the  observed  data  points. 


